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General Reverse Flow and Variational 
Theorems in Lifting-Surface Theory 


A. H. FLAX* 
Cornell Aeronautical Laboratory, Ine. 


SUMMARY 


) By introducing the integral equation and auxiliary conditions 
ladjoint to the integral equation and Kutta condition of lifting- 
surface theory, certain general theorems in lifting-surface theory 
Parederived. It is shown that the adjoint to the integral equation 
ofa lifting surface is the integral equation for the same surface in 
reverse flow and that to satisfy the requirements on such an ad- 
Fjoint, the Kutta condition must also be imposed on the solutions 
of the equation for the surface in reverse flow. The adjoint rela- 
tionship between wings in direct and reverse flow is shown to 
lead to a general class of relationships between the characteris- 
Mics of wings in direct and reverse flow for either subsonic or 
Supersonic speeds. Also, it is shown that the total lift, pitching 
Moment, and rolling moment on a wing of any camber and twist 
may be expressed as a closed-form integral over certain pressure 
Mistributions of the same wing in reverse flow. These distribu- 
tio s correspond to the case of the wing in reverse flow at con- 
Pstant angle of attack, pitching with uniform velocity, and rolling 

ith uniform velocity. The formulas may be regarded as gen- 
ralizations of the Munk two-dimensional airfoil integral relations 
to which they reduce for infinite aspect ratio. Finally, an ad- 
Mint variational principle for lifting surfaces in subsonic flow is ob- 
Mained. In combination with the form of the Rayleigh-Ritz 
Method, this is shown to lead to suitable methods for approxi- 
hate solution of lifting-surface problems. On the basis of the 
Wariational principle, certain apparently arbitrary or intuitive 
Steps in the approximate theories of Prandtl, Weissinger, Reiss- 
Mer, and Lawrence are shown to be logical applications of the 
Same method. 


SYMBOLS 


free-stream velocity 
local perturbation velocity in the free-stream 
direction divided by U 
local perturbation velocity in the free-stream 
direction divided by U’ for wing in reverse 
flow 
Cartesian coordinates in free-stream direction 
Cartesian coordinates perpendicular to x and 
lying in the base plane of the wing 
Received June 8, 1951. 
* Head, Aerodynamic Research Department. 


Cartesian coordinate perpendicular to x and 
y and measuring distance from the base 
plane of the wing 
M = free-stream Mach Number 
p(x, y) lift per unit area = 2pU°%u 
p(x, y) lift per unit area = 2pl/u for wing in reverse 
flow 
w(x, y) local slope of lifting surface in stream direc- 
tion 
local slope of lifting surface in stream direc- 


tion for wing in reverse flow 


w(x, y) 


free-stream mass density 

wing semispan 

wing chord 

distance from wing leading edge 

denotes integration over area of wing 
kernel of lifting surface integral equation 
kernel adjoint to K 

VM? - 1 

wing angle of attack 

wing pitching velocity (rad. per sec.) 


K(x, y: &n) = 
A(x, y: én) = 


dynamic pressure = (1/2)pl” 
wing rolling velocity (rad. per sec. ) 
variation of 


constants 


spanwise lift distribution function 
chordwise lift distribution function 


(1) INTRODUCTION 


ATTENTION HAS BEEN FOCUSED on the 


eens: 

existence in linearized wing theory of some rather 
unexpected relationships between the characteristics of 
wings in reversed flows, or reverse flow theorems as these 
relationships are sometimes called. von Karman! ap 
pears to have been the first to point out such a relation- 
ship. He showed that the supersonic wave drag due to 
thickness of a body would be unaltered by a reversal of 
direction of flight. Hayes* then showed, by other meth- 
ods, that not only would the drag be unaltered but also 
that, for a restricted class of plan forms, the slope of the 
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lift curve would remain the same for such a reversal of 
flow direction. Munk,’ using still another line of rea- 
soning, proved both the drag theorem and the lift 
theorem, the latter for a somewhat less restricted class 
of plan forms than Hayes, although he was not explicit 
in defining all the restrictions. The present author‘ 
showed, for wings with supersonic leading and trailing 
edges, that not only were there relations between lift 
and drag in reversed flows but also that were was a 
more general class of relationships between the sta- 
bility derivatives in reversed flows. Such relationships 
were extended to a somewhat less restricted class of 
plan forms in supersonic flow by Harmon.’ Brown, 
using an extension of Munk’s approach, demonstrated 
that the drag due to thickness, the slope of the lift 
curve, the damping in pitch, and the damping in roll 
were all unchanged by a reversal of flow direction for 
any plan form and either subsonic or supersonic flow 
velocities, as long as linearized wing theory could be em- 
ployed. 

The present paper presents a general form of the re- 
lationships between the characteristics of wings in direct 
and reverse flows, which includes all of those previ- 
ously known, as well as an infinity of others, and is 
applicable, within the limits of linearized lifting surface 
theory, to wings of any plan form in either subsonic or 
supersonic flow. The basis for the present result is the ob- 
servation that the reversalof flow direction hasthe effect 
that the kernelof theintegral equation of linearized lifting 
surface theory is the mathematical adjoint of the kernel 
appropriate to the original direction of flow. In order 
further to satisfy the mathematical requirements on 
adjoint solutions to these equations, it is shown that it is 
necessary that the Kutta condition be satisfied on the 
trailing edges in the reverse flow. Thus, an adjoint 
solution is completely specified as any solution of a 
lifting-surface problem for the same plan form with flow 
direction reversed. Among the interesting results of 
the generalized relationships obtained in this paper 
which are likely to be of practical value are generaliza- 
tions of Munk’s two-dimensional theorems for lift and 
moment to three-dimensional wings. The first of these 
gives the lift on a wing witlf any camber and twist dis- 
tribution (or specified normal boundary velocity result- 
ing from pitching, rolling, structural deformation, etc.) 
in terms of the lifting pressure distribution on a flat plate 
of the same plan form with the flow direction reversed. 
The second theorem of this type gives the pitching mo- 
ment on a wing with any camber and twist distribution 
in terms of the lifting pressure distribution on a pitching 
flat plate of the same plan form with the flow direction 
reversed. In the special case of two-dimensional flow, 
these theorems reduce to precisely the Munk result. 
For high aspect ratio wings, for which lifting-line theory 
is applicable, the adjoint solutions are simply any solu- 
tions of the lifting-line equation itself, and the general- 
ized Munk lift theorem gives the total lift on a wing 
with any twist distribution (due to geometrical twist, 


THE ABPRONAUTICAL 


SCIENCES JUNE, 1952 
structural deformation, flap deflection, etc.) jn terms 
of an integral over the lift distribution on an untwisted 
wing at constant angle of attack. A similar result 1s 
obtained for rolling moment on a lifting surface in 
terms of an integral over the lift distribution of a rolling 
wing in reverse flow. 

As a direct consequence of the adjoint relationship 
between the solutions to the lifting-surface problem 
and the solutions of the lifting-surface problem for the 
same plan form with the flow direction reversed, jt js 
found possible to formulate a variational principle of 
the adjoint type for lifting-surface theory. 
ational principle is of interest because it offers a system. 


This vari- 


atic means of obtaining approximate solutions for lift. 
ing-surface problems. Since exact solutions to the ip. 
tegral equation of lifting-surface theory in subsonic 
flow are rarely obtainable, good approximations are of 
great importance. In particular, the variational 
method may be used to obtain approximate one- 
dimensional integral equations to replace the two- 
dimensional integral equation of lifting-surface theory, 
for very high or very low aspect ratio wings. The 
theories of Prandtl,’ Weissinger,* Lawrence,’ and Reiss- 
ner’ are shown to be formulable on the same basis 
from the variational principle. These theories all in- 
volve multiplication of the two-dimensional integral 
equation by some weighting factor in one variable and 
integration over that variable, after a suitable approxi- 
mation to the kernel has been made. In the Prandtl, 
Weissinger, Lawrence, and Reissner theories, the choice 
of weighting function appears to have been made in- 
tuitively or by recourse to analogy with the limiting 
cases of two-dimensional theory or with very low aspect 
ratio theory (Jones). The variational method not 
only gives the form the weighting function must have 
but indicates how any particular approximation may be 
further improved. The two-dimensional integral equa- 
tion of lifting-surface theory may be approximated in 
this way by any finite number of one-dimensional in- 
tegral equations or by any finite number of algebraic 
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Because of the stationary property of the | y ” 
) 


equations. 
grals in the variational method, the accuracy of 


results for any given number of either one-dimensional 


inte 


integral equations or algebraic equations is considerably 
better than if the same number of equations had been 
obtained by simple point collocations, such as were 
carried out by Blenk"! in lifting-surface theory or by 


Glauert!? in lifting-line theory. In the limiting case of xé < [ 
3 t= 





jifting-line theory, the adjoint variational principle re- 
duces to an ordinary variational principle, which has 
been discussed by Ziller®* and the author.** 

A recent paper by Ursell and Ward’ presents the 
same general reverse flow relationship as developed in 
this paper but not the variational principle. Their 
line of attack is through the general equations for 
linearized flow potential and involves the use of mo- 
mentum and energy integrals over properly selected 
bounding surfaces in the fluid. The present method is 
somewhat simpler, starting from the established in- 
tegral equation of lifting-surface theory and retaining out. The variational principle then results almost as a 
a direct connection with lifting-surface theory through- direct consequence of the reverse flow relationship. 





Fic. 2. 


(2) THe INTEGRAL EQUATION OF LIFTING-SURFACE THEORY AND ITS ADJOINT IN INCOMPRESSIBLE FLOW 


The linearized theory of the lifting surface has been formulated in a variety of ways, but for the present purpose, 
the integral equation formulation is most direct. For concreteness, the particular case of incompressible flow will 
It will then be shown that the same considerations apply to compressible subsonic and super- 


be considered first. 
For incompressible flow, the form of the lifting-surface integral equation derived by 


sonic linearized wing theory. 
Lawrence’ from a result of Reissner’s!‘is most suitable. This is 


1 oO _ 5 u(&, n) V (x — &)? + (y — 9)? 
w(x, y) = / J dé dn = | + * . (1) 
: 2r Oy A y=-*% (x — &) 


where principal values are to be taken when necessary. The subscript A beneath the integrals denotes integration 
over the wing, and w(x, y) is the local angle of attack, which is assumed to be small enough everywhere so that the 
assumptions of linearized theory are valid. By the same token, the lifting surface is assumed to lie in a plane in- 
sofar as calculation of distances between elements or of integrals over the surface are concerned. w(é, n) is the 
jump in the streamwise component of the perturbation velocity across the wing surface made dimensionless with 
respect to the free-stream velocity LU. The lift per unit area is given as 
p(x, vy) = 2pU7u (2) 
The coordinate system is shown in Fig. | and is one in which the direction of positive x and £ is in the direction of 
U, the free-stream velocity. 
Now, if an integral equation adjoint’ to Eq. (1) is desired, it is necessary to interchange x and é and y and 9, re- 


spectively in the kernel. Thus, an adjoint might be 


lof ii(£, n) V(x — 22 + (y — 9)? 
M(x, y) = / ti dé dn : 1+ é = : | (3) 
. 2r Oy A yn &—*<x 


where ®(x, y) is any distribution of local angle of attackt and a(x, y) is the corresponding distribution of lift. Here 
itis to be noted that the y and n under the integral signs have not been interchanged, nor has the differentiation 
with respect to y outside the integral sign been changed to differentiation with respect to 7 in seeming violation of 
the definition of the adjoint equation. The.reason for this is connected with the singular nature of the integral, 
which does not permit the differentiation to be done under the sign. If this were possible, the kernel of the integral 


equation would actually be 
t To keep some of the expressions obtained positive in the sign of direct and reverse flow lifts, w is defined as the negative of the 


usual definition here and in all that follows. This amounts to reversing the usual direction of positive slope for the reverse 


flow 
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/ 
oj 1 V@e-8+0-9') 
< l An 4 
dy ly — | . (x — &) iF (4 


Since the term in braces is an antisymmetrical function of (y — 7), its derivative with respect to y isa Symmetrical 
function of (y — ) and therefore is unaffected by interchange of y and 7. This basic symmetry property of the 
operator (actually the sequence of operations: double integration followed by differentiation) is not affected by 
the singular nature of the integral as long as the functions u(x, y) and “(x, y) are not themselves singular. Ag js 
well known from the theory of lifting surfaces, these functions may become singular on edges, and this problem 
will be considered subsequently. 

Before proceeding with any further discussion, it will be useful to consider a physical interpretation of Eq. (3). 
This can readily be done by considering the lifting-surface integral equation appropriate to a reversal of flow direc. 
tion. Eq. (1) would still be applicable if, as shown in Fig. 2, the direction of positive x and £ were redefined so as to 
coincide with the new direction of the free-stream flow. If it is, instead, desired to retain the definitions of positive 
‘x and £ as shown in Fig. 1, it is only necessary to replace x by —x and by —é. This procedure leads directly to 
Eq. (3). Thus, Eq. (3) may be regarded as the integral equation of the lifting surface in a reverse flow when the 
definition of positive x and é is retained as in the original flow direction. 

It is well known that integral equations of types (1) and (3) do not have unique solutions without the specifica. 
tion of additional conditions on the solutions, u(x, y). In the case of Eq. (1), which is the lifting-surface equation, 
the specification of the Kutta condition at the trailing edge of the wing—namely, that u(x, y) be zero at such an 
edge—is sufficient to make the solution unique. But nothing up to this point has imposed any such condition 
on the solutions of Eq. (3). Therefore, it is not possible to state what solutions for a(x, y) are admissible as solu. 
tions of the adjoint equation. If the admissible solutions of the adjoint integral are actually to be those corre- 
sponding to the lifting surface in reverse flow, it is necessary to show that an admissible solution of the adjoint 
equation must satisfy the Kutta condition on the trailing edges in the reverse flow—i.e., on the edges that were 
leading edges in the original problem. This will be done next. 


To simplify the following discussion, Eqs. (1) and (3) will be written, respectively, as 


w(x, y) = fa u(é, n)K(x, vy: &, n) dé dn 5 
A 

@(x, y) = J a(é, n) K(x, y: & 9) dé dn ( 
A 


+ . ° . ° . . - . . 
where the symbol ‘f° stands for the operations of double integration over the area, followed by partial differen- 
tiation with respect to y. Now, in order that the solutions of the adjoint equation satisfy the usual mathematical 
requirements,'° it is necessary that 


f u(x, y) J u(é, n) K(x, y: &, n) dé dn dx dy = f u(x, y) / u(t, n)K(x, y: &, 9) dé dn dx dy 7) 
A A A JA 


where the integral symbol, ~, denotes double integration over the area of the wing for any admissible func- 
tions u(x, y) and a(x, y). It can be seen that the expression on the right-hand side is essentially that on the 
left, with the orders of integration over &, 7 and x, yinterchanged. The operation of differentiation with respect 
to y, which is interposed between the integrations, causes no difficulties because of the symmetry of the entire 
operation in y and n, as discussed previously. In essence, ‘/ K and *f K are identical with regard to y and y. This 
is brought out more clearly in Appendix (II), where, as a result of some manipulation, the symmetrical form in 
y and 7 of the expressions in Eq. (7) is made more apparent. Thus, if the interchanges in order of integration 
are carried out on the right-hand side of Eq. (7), the expression becomes 


Z u(&,n) J u(x,y) K(x, y: &,n) dx dy dé dy s 
A A 


*{ now denoting integration with respect to x and y followed by differentiation with respect to 7. Then it may 
be seen that expression (8) is identical to the expression on the left-hand side of Eq. (7), since x, y, &, and 7 are 
merely integration variables and do not appear in the results. 

It remains to consider under what conditions an interchange of orders of integration is permissible. If the fune- 
tion f(x,y,,n) = a(x,y)u(é,n) is everywhere bounded, continuous, and differentiable, including the vicinity of 
points where x = £ and y = 7, then such an interchange is justifiable, even in spite of the fact that the integrals 
are only definable in terms of their Cauchy principal value (see, for instance, Schwartz"). Now, the function 
u(x, y), representing the local lift distribution over a lifting surface, will satisfy this condition everywhere except 
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possibly near the edges. Near an edge, u must either go to zero or become infinite as 1/+/d, where d is the distance 
from the edge, if finite values of the local angle of attack are to be obtained from Eq. (1)."" The Kutta condition 
specifies that u(x, y) go to zero on trailing edges; on leading edges, it goes to infinity as 1/+/d, where d is the dis- 
tance from the leading edge. Therefore, the conditions for interchange of order of integration as in Eq. (8) above 
will be satisfied if it is required that a(x, y) go to zero on the leading edges of the original wing—i.e., the trailing 
edges of the wing in reverse flow. But this is exactly the requirement of a Kutta condition on the adjoint solu- 
tions for the wing in reverse flow. Thus, a sufficient condition for the specification of an adjoint solution a(x, y) 
is that it satisfy the integral Eq. (3) for a lifting surface in reverse flow and that it satisfy the Kutta condition at 
the trailing edges in reverse flow. 

In order that the adjoint solution a(x, y) corresponding to a given @(x, y) be unique, it is further required that 
it be shown that the satisfaction of the Kutta condition is not only sufficient to specify the adjoint but also is 
necessary. This is done in Appendix (I), where it is shown that the interchange of orders of integration required 
in order to prove Eq. (7) is not permissible.if both u(x, y) and a(x, y) become infinite as 1/+/d at any edge. The 
difficulty arises from the singular behavior of the integrals (7) near such an edge. Ifu ~ C/~/dand ii ~ C/Va 
in the vicinity of the edge, then the leading-edge singularity contributes a term proportional to CC to the integral 
for one order of integration and proportional to —xCC for the reverse order. The form of these expressions and 
the fact that they arise from edge singularities suggest that the difficulty arises from the same source as requires the 
insertion of the leading-edge suction into calculations of drag from lifting-surface theory. This is indeed the case; 
in fact, if @ is taken to be u, the expressions (7) become the drag integrals of lifting-surface theory, and the term 
under consideration becomes 1C°, which is exactly the term annulled by leading-edge suction." 


(3) THe INTEGRAL EQUATION AND Irs ADJOINT FOR COMPRESSIBLE FLows 


For subsonic compressible flows, it is apparent that the application of the three-dimensional Prandtl-Glauert 
transformation ® to both the integral equation and its adjoint changes both in exactly the same way, so that all the 
results previously derived for incompressible flow are applicable without further discussion. 

For supersonic flows, some further investigation is necessary, since all points on a lifting surface do not mutually 
interact because of the propagation of disturbances from a point only inside the downstream Mach cone from the 
point. In turn, a given point is influenced only by disturbances emanating from points within its upstream Mach 
cone. Using the system of axes shown in Fig. 1, with the addition of a z axis at right angles to the wing surface 


which forms the plane z = 0, the expression for w is® 


a it (é, n) dé d th (9) 
w= hm Us, N) dE ayn 9 9 9 » 2 2)t'/2 ‘ 
: 0 F O02 ( [(y — n)? + 27] { (x — t)?— B[(y — n)? + 2 } 


where 8 = +/ \/? — 1, and the area of integration c consists of points on the wing in the Mach forecone of the point 
tn(x <£). This is the region in which the radical in the denominator of the integrand remains real. It may be 
seen that the kernel of this integral equation has the same basic symmetry properties as the subsonic kernel of Eq. 
(1)—namely, that it is symmetrical in y and 7, and antisymmetrical in x and —since the operation of partial dif- 
ferentiation with respect to z has no effect on these properties. 

The adjoint integral equation is _— obtained by interchanging x and é and y and 9, respectively. This gives 


u(t, n) dt, dn(é — x)z 
w= ta > = : arr. (10) 
Tv -—>0 Oz ‘(iy aaa n)? + 27] { (x — ¢)? _ B*l(y = n)- + 2"); : 


The integral is now to be evaluated over the area c’(é > x) in accordance with the rule that the roles of x and & are 
to be interchanged. It is seen that, as in the subsonic case, the adjoint integral Eq. (10) is simply the lifting- 
surface equation corresponding to a reversal of flow direction (Fig. 3). 


Before proceeding, it will again be convenient to abbreviate the notation by writing 


w= J (£, n) K(x, v: & ) dé dn (11) 


= fu n)K(x, y: & 9) dé dy (12) 


where now ** {© denotes double integration, followed by partial differentiation with respect to z and then by a 
limiting process as 2 goes to zero. Then, again, to meet the requirements on the adjoint solutions a(x, y), 


f u(x, y) / u(t, n) K(x, y: &, 9) dé dn dx dy = Z u(x, y) I a(t, n)K(x, y: &, ) d& dn dx dy (13) 
4 F A c’ 
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This will again be true if an interchange of orders of integration in the repeated integral on the right-hand side ‘. 
permissible, since then the expressions on the two sides will be identical, x and ¢ being integration indices only, At 
first glance, the possibility of such an interchange may seem surprising, since the integral Eqs. (11) and (12) involve 
integrations over different areas (c and c’). Referring to Eq. (13), however, it may be seen that the interchange 
of order of integration corresponds physically merely to calculating the product of #(£, ») and the downwash pro- 
duced at each point (é, 7) on the area c’ by the doublets of strength u at a given point (x, y), integrating over é and 
n, and then integrating over all x and y, rather than calculating the downwash produced at point (x, y) by aj] 
doublets of strength u(é, 7) in area c and then multiplying by a(x, y) and integrating over all x and y. Clearly 
these two methods of integration should be equivalent as long as the function f(x, y, —, 7) = u(x, y)a(é, 7) is bounded, 
continuous, and differentiable. If there are subsonic leading edges, u becomes infinite as 1/+/d as in the purely 
subsonic case and again # must go to zero at such points for the interchange of order of integration to be per- 
missible. Thus, at the subsonic trailing edges of the reversed flow (i.e., subsonic leading edges of the original 
flow), the Kutta condition must be satisfied. This condition is again just sufficient to completely specify the ad. 
joint or reverse flow solutions 7. 


(4) THE GENERAL RELATIONSHIP BETWEEN REVERSE FLOWS 


The adjoint relationship between reverse flows leads immediately to a general relationship between certain signifi- 
cant integrals in such flows. If Eq. (5) is first multiplied by #(x, y) and integrated over the wing area, it be- 
comes 


/ w(x, y)u(x, y) dx dy = J u(x, y) J u(f, n)A(x, y: & 9) dé dn dx dy (14) 
JA A A 


Similarly, multiplying Eq. (6) by u(x, y) and integrating over the wing area gives 


J @(x, y)u(x, y) dx dy = I 
A A 


By relationship (7), the right-hand sides of these two equations are equal, so that 


/ w(x, y)u(x, y) dx dy = y) W(x, y)u(x, y) dx dy (16) 
JA A 


By using Eqs. (11) and (12) in place of Eqs. (5) and (6), Eq. (16) may be seen to be true for supersonic flow as 


a > 


* 
u(x, y) J u(é, n)K(x, y: & n) dé dn dx dy (15 
A 


well. Eq. (16) is the general reverse flow relationship that has been obtained by different means and in a slightly 
different form by Ursell and Ward."* 


(5) Sprecrric RELATIONSHIPS BETWEEN CHARACTERISTICS OF WINGS IN REVERSE FLOW 





(A) Lifts 


Ly7 


Letw = aand® =a. Then, 


* 


a J u.dx dy = a | u dx dy (17) 
A A 


or the lifts of a wing in direct and reverse flow are equal 
at equal a. 


(B) Pitching Moment and Lift Due to Pitch 


Let w = aand ® = gx/U (wing pitching about x = 0). 


Then, 
a . . 
a idx dy = LV ux dx dy (18) 
A A 


Thus, for unit a and q/U, the integral of the lift, 7, on a 








nb 


wing pitching about + = 0 is equal to the integral of the 





moment of the lift, wx, on a wing at constant angle of 
attack. Fic. 
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LIFTING 


(\C) Damping in Pitch 


Let w = gx/Uand ® = gx/U. Then, 


‘ | . = 
- ux dx dy = — ux dx dy (19) 
U fa l A : 


Thus, the damping moments on a pitching wing are the 
same for a wing in direct or reverse flow. A similar 
relationship for damping in roll obviously results if 


wand @ are taken as py/U. 


(D) Drag Due to Lift 

If a lifting surface at constant angle of attack, a, 
has only supersonic leading and trailing edges so that 
there is no leading edge suction in either direct or re- 
verse flow, then D = La, and the equality of lifts al- 
ready demonstrated is sufficient to demonstrate the 
equality of drag due to lift in direct and reverse flow. 
It should be noted that side edges parallel to the free- 
stream flow are permitted, and any number of interact- 
ing regions from the tips will not alter the result. If 
subsonic leading edges are present in either the direct 
or reverse flow or in both, the drags due to lift will differ 
by the difference of the leading-edge suctions in the two 


cases. 


E) Yawing Moment Due to Rolling 

Again, if the wing has no subsonic leading edges in 
either direct or reverse flow, the yawing moments due 
to rolling will be equal, since at any point the yawing 
moment contributed will be a times the rolling mo- 
ment, and the equality of the damping in roll for direct 
and reverse flow has already been shown. 


F) Rolling Moment Due to Sideslip or Differential Wing 

Deflection 

In subsonic flow, the effect of sideslip on a wing with 
dihedral and small sweep or of differential wing deflection 
is to produce an increment of angle of attack which is 
positive on one wing and negative on the other. This 
may be represented by w = Aa sgn y, where a is the 
differential angle of attack and sgn y denotes the sign 
of y. 

Let w = pyand ® = Aasgn y. 


Pp [ ay dx dy = Aa J usgnydxdy (20) 
JA A 


For a spanwise symmetrical wing, the integral on the 
right-hand side is simply the sum of the absolute values 
of the lift on the two sides of a rolling wing. Thus, for 
unit p and Aa, the rolling moment due to Aa is simply 
equal to the sum of the absolute values of the lift on a 
rolling wing. For supersonic flow, the Aa must be 
considered as arising from differential wing deflection 
only, since sideslip produces, in general, changes in re- 
gions bounded by the Mach cones from the tips, as well 
as changes in effective leading-edge Mach Number 
which may be important. 


Then 
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(6) GENERALIZATIONS OF THE MUNK Two- 
DIMENSIONAL AIRFOIL THEOREMS 

The general reverse flow relationship (16) may be 
applied to a more general class of problems to give 
formulas for total lift, pitching moment, and rolling 
moment on a wing with any specified boundary condi- 
tions, w = f(x, y), which reduce, for the two-dimensional 
case, to the well-known Munk formulas for lift and 
pitching moment. 


(A) Lift Formula 


Let w = aandw = f(x,y). Then, 


ay 


a J udx dy = a) uf(x, y) dx dy 
A A 
4 < = 
a) udxdy = y) f(x, y) dx dy (21) 
A A @ 


Now the lift, L, is given by 


L= zou f u dx dy 
A 
~ fet. 
L = pl | f(x, y) dx dy 
Qa 
L= f P f(x, y) dx dy (22) 
A a@ 


where is the local lift per unit area on the wing in re- 
Thus, the lift distribu- 


or 


so that 


or 


verse flow at angle of attack a. 
tion on the wing in reverse flow at constant unit angle 
of attack is the influence function for the total lift on 
the wing in direct flow. This formula is particularly 
useful for wings symmetrical chordwise about a span- 
wise axis, since then (x) is simply u(—.x), and no ac- 
tual calculation with the wing in reverse flow need be 
done. Also, in lifting-line theory, every wing is the 
same in direct and reverse flow, leading to a simple 
application of Eq. (22) to flapped and _ twisted 
wings 

In the limiting case of two-dimensional thin air- 


foils,?! 

u/ja = VL — (2x/c)]/[1 + (2x/c)] 
Thus, substituting —.x for x, 

fi/a = —V[1 + (2x/c)]/[1 — (2x/c)] 


and 


“el? ! + (2x/c) 
L = —29l/* f(x) @/ dx 
4 | 2 INT (0x/0) 


where c is the chord and f(x) is the local slope of the air- 
foil. This is exactly Munk’s formula for lift. 
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(B) Pitching Moment Formula 


Let ®@ = g(x/U) and w = f(x,y). Then, 


: J ux dx dy = I uif(x, vy) dx dy (23) 
U Ja A 
M1 
F ux dx dy = ri l f(x, vy) dx dy 
A A q 


0, A is given by 


or 


The moment about x = 


M, = 2pU? / ux dx dy 
A 
" mM . 
M, = 2pU* rf I(x, y) dx dy 
A q 


so that 


or 


(7) THE 
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Thus, /,U/q, where /;, is the lifting pressure distribution 
on the wing in reverse flow pitching with angular ye. 
locity g about the axis x = 0, is the influence function for 
total moment about x = 0 of the wing in direct flow. 
Again, this reduces to Munk’s formula in two-dimen- 
sional flow as follows:*° 


pi/q = pcU V1 — (2x/c)? 


so that 


~ 2x\? 
M, = ct? f f(x) \! = (= ) dx 


which is one form for the moment about the airfojl 


center according to Munk.” 


(C) Rolling Moment Formula 


Let ®@ = py/U and w = f(x,y). Then, 


P 
: uy dx dy = uf(x, vy) dx dy (25 
USA A 
vy 
uy dx dy = ul 
A A p 


The rolling moment is given by 


My, = 2pU” / uy dx dy 
A 

— f#® 

a p 


Here it is seen that p.l’/p plays the role of influence 


coefficient for total rolling moment, where /» is the 
lifting pressure distribution on the wing in reverse flow 


or 


f(x, vy) dx dy 


or 


f(x, y) dx dy (26 


rolling with angular velocity p. 


ADJOINT VARIATIONAL PRINCIPLE 


The results that have been obtained establishing the adjoint relationship between direct and reverse flows over 


a lifting surface [Eq. (7 
flows, 


7)| may be used to formulate a variational principle for such flows. 


Consider, for subsonic 


6] = | fi (mu + wi) dx dy - fi if. Ku dé dn dx ay | (27 
A 


where 6 denotes the first variation in the usual sense of the calculus of variations and wu and @ are restricted to ad- 


missible functions satisfying the Kutta conditions in the direct and reverse flow, respectively. 
limitation, the functions u and 4 are to be varied arbitrarily and independently of each other. 


Subject to this 
Then, 


a= f (wou + win) dx dy — f éu I Ku dé dyn dx dy — fi 4 Kéu dé dn dx dy (28) 


But, by Eq. (7), 
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LIFTING-SURF 


J it J Kéu dé dn dx dy = 
. A 


4 
Substituting Eq. (29) into Eq. (28), 


i] = jf @— / Kit dé dn | bu dx dy + J | w — J Ku dé dn | 6% dx dy 
A JA A 4 


ACE 
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J éu J Ku dé dn dx dy (29) 
A A 


(30) 


If wand # are the correct solutions corresponding to w and @ in direct and reverse flow, respectively, then 


@— I Ki dt dn = 0; w-— J Ku didn = 0; 6] = 0 
A A 
6 f (mu + wi) dx dy — J 7 / Ku dé dn dx dy | = ( (31) 
A A A 


is equivalent to the integral equations of lifting-surface theory for the direct and reverse flow over a given plan 


Thus the expression 


form with local slopes w and @, respectively. 


The value of such a formulation as opposed to the integral equa- 


tions themselves is that it lends itself most readily to approximate solutions by methods analogous to the Ray- 


leigh-Ritz method for self-adjoint problems.*! 


An important feature of solutions by the variational method is 


that the expression in brackets in Eq. (31) is stationary for first-order variations of u and # from the correct values. 


Thus, the errors in the expression itself tend to be minimized even with approximate u and i. 
expression are closely related to lift, moment, and other qualities of physical interest. 


in greater detail in the following section. 


The integrals in the 
This matter will be pursued 


(S) APPROXIMATE METHODS FOR SUBSONIC LIFTING-SURFACE PROBLEMS 


On the basis of the adjoint variational principle, Eq. 
(31), various approximate methods for the solution of 
lifting-surface problems may be formulated. These 
methods are essentially similar to the Rayleigh-Ritz 
method, which is commonly used for self-adjoint prob- 
lems. For most practical purposes, it is usually desir- 
able to simplify the two-dimensional integral Eq. (7) 
of lifting-surface theory in at least two ways: first, by 
reducing the equation to a one-dimensional integral 
equation or a finite set of these; second, by intro- 
ducing more readily integrable approximations to the 
kernel, A, appropriate to the wing shapes and aspect 
ratios being studied. The variational principle pro- 
vides a direct systematic way of accomplishing the first 
of these simplifications; the second simplification is, of 
course, still necessary. The Prandtl,’ Weissinger,* and 
Reissner"’ lifting-line theories for relatively high aspect 
ratio wings will be shown to follow from the variational 
principle plus the high aspect ratio kernel approxima- 
tions, the nature of which have already been clearly 
brought out by Reissner." For low aspect ratio wings, 
the Lawrence theory will be seen to follow from the 
variational principle plus the low aspect ratio approxi- 
mation to the kernel presented by Lawrence.’ In all 
of this work, the two-dimensional integral equation was 
multiplied by some weighting function in one variable 
and integrated over that variable after a suitable ap- 
proximation to the kernel had been made, thus reducing 
the problem to that of solution of a one-dimensional 
integral equation or, in the case of Reissner’s work, to 


two one-dimensional integral equations. For the theo 
ries of this type currently available, the choice of the 
weighting function appears to have been made intui 
tively or by recourse to analogy with the limiting cases 
of two-dimensional theory or of very low aspect ratio 
theory.”” The variational principle leads to all of these 
weighting functions on a unified basis and indicates, 
as well, further approximations to lifting-surface 
theory. 

It should be remarked here that, in principle, the 
two-dimensional integral equation of lifting-surface 
theory can be solved by collocation at a sufficient num- 
ber of points spanwise and chordwise over the wing. 
This, however, usually turns out to be extremely cum- 
bersome, and in few cases is such a procedure used in 
practice. Further, if a small number of points are to 
be used in collocation, it is usually to be expected that 
the variational technique involving a number of free 
parameters equal to the number of collocation points 
would be considerably more accurate. This follows 
from the fact that the variational method involves 
averaging of the satisfaction of the boundary condi 
tions over the wing surface with a weighting function 
such as to make certain integrals closely related to the 
total lift and moments on the airfoil stationary with 
The 
proper weighting functions are, of course, according 
to Eq. (31), the adjoint lift distributions. Some of 


these ideas will become clearer in specific applica- 


respect to first-order errors in the free parameters. 


tions. 
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(A) High Aspect Ratio Theory 


For relatively high aspect ratios, the spanwise distribution is of principal interest, so that a suitable approxima. 
tion for the lift distribution is of the form!" '” 


(c/2) — § i A oc een dO 2) 
u(é, n) = Adi(n) Vic 2) +¢ + Bo2(n) yi - (*) + C¢3(n) (Oo) l — wy, +... (3: 


If the wing is symmetrical about a spanwise axis, reversal of flow direction involves only a replacement of £ by 


—£, so that 


~ . _ |(c/2) + 2\? 8 2¢ 2¢\? 
u(g,) = Adin) -\ (5) = : + Box(n) 1 — (*) — Cos(n)- (“)y! - i ee (3: 


Here, c, the wing chord, may be a function of n, and the axis £ = Ois the axis of symmetry. The constants A, B, 
C, etc., multiplying each term are the free parameters to be established on the basis of the variational principle, 
while the ¢’s are unknown functions to be determined on the basis of the one-dimensional spanwise integral equa- 
tions obtained. The chordwise functions are those originally introduced by Birnbaum"! in two-dimensional wing 


theory; the first, Vv (c/2) — §]/[(c/2) + &], gives the chordwise distribution for constant angle of attack, w = 
g 


constant; the second, V1 — (2€/c)?, gives the chordwise distribution for wat; the third, 2&/c V1 — (2£/c)?, 
gives the chordwise distribution for wa|(2é/c)? — 1/2] (the s-shaped airfoil). As in the Rayleigh-Ritz method, 
6u = (0u/0A)6A + (00/0B)6B (34) 


and since 6A, 6B, 6C, etc., are arbitrary, each variation may be treated as independent. 
To obtain a single one-dimensional integral equation corresponding to lifting-line theory, only the first terms of 
the series Eqs. (32) and (33) are used in conjunction with Eq. (31). Thus, in this case, 
6u = (Ou/OA) 6A 
du = (Ou/OA) 5AS 


Eq. (31) then leads to 


sf (c/2) + x [xa (c/2) -— 
f ’ v — Z ) 
_” YN 2) -—x ’ A a V(c/2) + 


0/2) — x 


7 
dé an| dx dy + 6A f Ply) \ : 2) > & " 


oo (c/2) 
[2 — / KA i(n) Vi. 2) Ms dé in| dx dy = 0 (36) 


Since 6A and 6A are independent and since in this case only the wing in direct flow is of interest, the final result 


See Ste 


tre Str 


desired is obtained by setting the term multiplied by 6A equal to zero, 


ji (y) l(c 2) + “| ' ' sg KAg( (c/2) — 
) wx, Vv aa Z 
A Pty Vic 2)-—x ; , A OW Ni 2) + 


The integration over y may be dropped, since for arbitrary function, ¢), the entire integral will be zero if the con- 


(c/2) +x ri (c/2) — 
wie, y) — | Ks 
INE 2) — =| * ¥) , PAGO VN Cia) a 


is satisfied. The problem now resolves itself into finding a suitable approximation to the kernel, so that the inte- 
gration with respect to and x may be readily carried out. Reissner' has discussed the nature of such approxima- 
tions for the Prandtl and Weissinger lifting-line theories. When the integrations on £ and x have been carried out, 


a one-dimensional integral equation for the spanwise lift remains. 


dé dn | dx dy = 0 (37) 


Sere | Str 


dition 


dé in| dx = 0 (37a) 


Str | Str 


By retaining the first and second terms of the series Eqs. (32) and (33), two simultaneous one-dimensional in- 
tegral equations for the spanwise lift, as originally given by Reissner, are obtained. In this case, 


di = (0u/0A) 6A + (00/0B) 6B (38) 

j I(c/2) —€& 2&\? 

£ =e Ss i) —_ 39) 
u(é, n) = Adi(n) \ (c/2) + + Bd»o(n) \ l % ) (3 


Further, 
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0n/O0OA = V [(c/2) + 
ou/OB = vi ~ (2 


E]/[(c/2) — =] di(n)/ = 
» ( ) 


c)* de(n) \ 


Then, using Eq. (31) and retaining only the equations applying to the wing in direct flow, 


. 


[ (*) j / x[a (c/2) — 
»(y) — Wx, ¥) — \ 4 () 
JA PAY) 'N é JA un Vic 2) + 


Once again, with an approximation for A which per- 
mits the integrations on x and & to be carried out, two 
one-dimensional integral equations result. These may 
be solved simultaneously for A¢,(y) and Bd¢2(y). 
the variational principle, the procedure for generaliza- 


Using 


tion to m simultaneous equations is obvious. 

In lifting-line theory, the restriction of the above pro- 
cedure to wings with a spanwise axis of symmetry is not 
essential, since the chordwise lift distribution is treated 
very approximately. Therefore, the above procedure 
may be carried over to any high aspect ratio wing of 
small sweep. 

Once a one-dimensional integral equation has been 
obtained, it in turn may be solved by collocation or by 
a variational technique. The application of the vari- 
ational technique to lifting-line theory has been carried 
out by Ziller** and by the author” and has been shown 
to be identical to the Lotz method if the chord and twist 
distributions are represented as Fourier series of the 
type used by Lotz. Thus, the Lotz method in Prandtl 
lifting-line theory may be viewed from the standpoint 
of the variational method as the application of Eq. 
(31) along with the Prandtl kernel approximation, and 
with the functions u and # given by 


|(c/2) — x a ee 
u(x, y) = \ A,y'Vs-— y 
i (c/2) +x f : 
+. 


-/9 ‘ 
| _ Mc/2) +x Oe : 
tix, 7) = \ A,y’Vs-—y 
< (c/2) —x 4 j 


where s is the semispan. The success of the Lotz 
method may then be regarded as an indication of the 
power of the variational technique as well. 


B) Low Aspect Ratio Wings 


Low aspect ratio wings have been treated by Law- 
rence.’ His methods may be shown to fit into the 
framework of the variational method and may be gen- 
eralized accordingly. Lawrence deals only with straight 
trailing edges; then, for sweptback leading edges, 
functions w and # having the proper edge conditions 
are, as used by Lawrence, 


[ l(c 2)+ xj ' / K| Aes (c/2) — 
(vy) wx, 7) => é (9) 
{ ane Vic 2) —x { ‘ A in Vic 2) + 


¢ 2¢\2 
g _ ([< / — 
: + Bodo(n) \ l (*) dé dn dx dy = 0 
(41) 
2¢\° { 
; + Bo»(n) \! — (*) | dé dn , dx dy = 0 
(42) 
l 
u(x, vy) = A £')| 
s?— y? (44) 


u(x,y) =AVs? — y? e"(x)| 


where s is the span, which may be a function of x. It 
should be noted that for the wing in direct flow (1) 
there is a 1/\/d singularity at the leading edge, while 
for the wing in reverse flow, (7), the lift goes to zero at 
that same edge as required by the Kutta condition. 
Application of Eq. (31) leads to 


I. s? — y? g’(x) wr y)=— s K X 
SA 1 


A 


g’(x) dé in| dxdy = 0 (45) 
Vs -7” 


By a kernel approximation suitable for low aspect 
ratios, the integrations on y and 7 can be carried out 
leaving a one-dimensional integral equation in the 
chordwise direction. In this case, the integration on 
x may be dropped. 

The foregoing procedure may be generalized to ob- 
tain any number of simultaneous one-dimensional in- 
tegral equations by using the functions 


) Any" 
u(x, y) = = e.(2) 
Vs? — y? 


u(x, V) 


(46) 


\| 
a=) 
4 

| 
7S 


Terms with 2 even correspond to symmetrical spanwise 
loading, while those with m odd correspond to antisym- 
metrical loading. Proceeding by analogy to the ex- 
tremely low aspect ratio case (Jones**), Lawrence’ 
has utilized the term n = | to investigate damping in 
It is perhaps worth re- 


4 


roll of a low aspect ratio wing. 
marking, for the sake of completeness, that the vari 
ational procedure may also be applied to the Jones 
theory, although in this case the integral equation is so 
simplified that exact solution by inversion of the equa 
tion is possible.” 
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In exactly the same manner as illustrated for the 
several cases discussed above, the variational method 
may be employed to construct solutions in subsonic 
lifting-surface theory to any desired degree of accu- 
racy. In particular, sets of simultaneous one-dimen- 
sional equations may be obtained, or alternatively, the 
two-dimensional integral equation may be reduced to 
any number of linear algebraic equations for series 
representations of the solutions involving any number 
of spanwise and chordwise free parameters. The feasi- 
bility of this procedure is limited, however, by the prac- 
tical difficulties involved in performing some of the 
required integrations. For this reason, approxima- 
tions to the kernel of lifting-surface integral equations 
are of the greatest importance. 


(9) CONCLUSIONS 


Through the introduction of the concept of the in- 
tegral equation and auxiliary conditions adjoint to the 
integral equation and Kutta condition of lifting- 
surface theory, a unified approach to many of the prob- 
lems of the theory has been effected. It has been shown 
that the adjoint to the integral equation of a lifting 
surface is the integral equation for the same surface 
in reverse flow, and that to satisfy the requirements on 
an adjoint, the Kutta condition must also be imposed 
on the solutions of the equation for the surface in re- 
verse flow. The adjoint relationship between wings in 
direct and reverse flow has been shown to lead directly 
to a general class of relationships between the charac- 
teristics of wings in direct and reverse flow for either 
subsonic or supersonic speeds. In addition, generaliza- 
tions of the Munk formulas for lift and moment in two- 
dimensional flow to three-dimensional flow have been 
obtained. Finally, an adjoint variational principle for 
lifting surfaces in subsonic flow has been obtained. 
On the basis of this variational principle, certain ap- 
parently arbitrary or intuitive steps in the approxi- 
mate theories of Prandtl, Weissinger, Reissner, and 
Lawrence are shown to be logical applications of the 
same method. 

Further generalizations of this method have been indi- 


cated. 


Appendix (1) 


CONDITIONS FOR INTERCHANGEABILITY OF ORDER OF 
INTEGRATION IN Eg. (7) 


Since the function f(x, y, & 1) = u(x, y)u(é, n) is 
bounded, continuous, and differentiable except near 
leading edges, only the contribution to the integral of 
points near such edges need be considered. Near an 
edge, the flow is essentially two-dimensional, so that 
the integrals corresponding to Eq. (7) for two-dimen- 
sional wings will suffice for this investigation. This 


integral is of the form 


FCAL SCIENCES 


1 ' i , u(é) 
i on J u‘x) J ee Pas dx (I-1) 


The lift distribution in the direct flow, u(), will have Q 
1/+/d singularity with distance d from the leading edge, 
£ = —1, and will go to zero at the trailing edge, ¢ = 
+1." Thus, the function f(x, &) = a(x)-u(t) will 
be regular at all points except possibly at the leading 
edge, where 7(£) may also behave as 1/+/d. The re. 
peated integration, Eq. (I-1), may then be considered 
to be a double integral in the x, & plane if the region 
x= —lto—1+e£&= —1 to —1 + € is excluded, as 
shown in Fig. 3, since the function f(x, £) is then regular 
over the region of integration. In this case, the inter- 
changeability of order for the integral (I-1) has been 
shown by Schwartz.'® It remains then to consider the 
excluded region, or 


—1 -_ i 
1 Te . 1+e u(é) 
i= = u(x) dé dx 
or J! —1 x—€& 


Letx = —1+Aandé = —I1 + uw. 


l Ee © u(p) 
= 1 f u(A) f "ee __ dx (I-3) 


Now, suppose near the leading edge (u = A = 0) 
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(I-2) 


Then, 


u(u) ~ 2C/Vn) 


i =~ 20 Vv ie 


so that 
,=4 ef f 1p dd 
P= or Jo V/r 0 V/ ur “— 


Integrating with respect to u, 
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Expanding the logarithmic term in series, 
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sonic leading edges in the direct flow. But these are 
+ | the subsonic trailing edges in the reverse flow, so that 
this requirement is exactly the Kutta condition for the 

cc 1 it + reverse flow. 
i=5 : > (Qn + 1)2 => .- - The singular behavior of these integrals near the 
leading edge is easily seen to be related to the effect of 
the leading-edge singularity on drag in lifting-surface 
Finally, theory, which must be removed by the introduction of 
I, = rCC (1-6) leading-edge suction calculations into the theory. 
ange a, The drag integral for a two-dimensional lifting surface 
If the order of integration in Eq. (I-5) is interchanged, 


the integral becomes 
‘ ' u(é) 
~~ £4 yg 1 fat J u(x) pi é dé dx (I-8) 
y| , I dy du 2a J-1 —s=—s§ 
T 0 Vyu/o VX A— up 


I,’ =2 

with u(x) ~ 2C/V1 + x near x = —1. Except for 
Clearly, this only has the effect of changing the sign of the region near the leading edge u(x)-u() is regular, 
the integral, since the remainder of the steps go through and the double integral for the remainder of the region 
as above, so that of integration must be zero because of the antisym- 
i! = eff (1-7) metric character of the integral. For the region near 
the leading edge, the previous integration gives, accord- 

Thus, the order of integration may not be inter- ing to Eq. (I-6) with C = sd 
changed unless either C or C is zero at the edge. Ata D= orl? -9 
subsonic leading edge of the wing in direct flow, C is ~ (s-9) 
obviously not zero, so that C must be, if the interchange This is exactly the amount of drag annulled by the 
of order of integration is to be permissible. Therefore, leading edge suction,” so that the drag of a two-dimen- 
the requirement on the solution of the adjoint integral sional wing is zero, which is required by perfect fluid 

equation is that the lift go to zero at edges that are sub-__ theory. 


Appendix (IT) 
SYMMETRIC FORM OF INTEGRALS OVER ADJOINT FUNCTIONS 


The symmetry in y and 7 of the integrals in Eq. (7) is not immediately apparent from the forms used in the 
main body of this paper. It was pointed out that the restrictions on differentiation of an infinite integral under 
the sign caused the apparent lack of symmetry. The symmetry can readily be brought out by a careful applica- 
tion of the rules for differentiation under the sign. The integrals in Eq. (7) are of the form 


9 


. 3 . 1 1 V(x — 8? = 97" 
I, = ii u(x, y) IS u(&, n) | + “ $) = Yy ") | dé dn dx dy (II-1) 
A ov A y>-y °° =—% (x — §) 


where on all edges either u or % is zero or both are zero [Appendix (I)]._ Now consider 


» vA 
ra) l l V (x — &)? »— 7)? 
I, = If u(é, 7) | i e 7 +9 : |ae dn (II-2) 
Oy A yoy yn (x — &) 


Denote the local semispan as b(£); this will be the limit of integration on y. If a small region » = y + € is ex- 
cluded from the integration, the differentiation under the integral sign may be performed, provided due attention 
is given to the variable limits of integration. This corresponds, of course, to recognizing the definition of the 


integral according to its Cauchy principal value. Thus, 


" os sr 1 V(x — &)? + (y — 9)? 
I; = lim / dé if dn u(é, n) | + ; é . 2 + 
+0 < Ov tJ. y-7 y—- (x — &) 


* 1 V(x — &)? »— 9)? 
/ dn ue») | + s é ua, : | (II-3) 


JSJyte 7-—¥y y= 9% (x — &) 


Root . . . . . 
Performing the differentiation under the sign, 
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; 2u(é, y Vix -—e?+e as ] (x — &) 
I; = lim / de 4 : aE - ~— | - 7, dn aL - a 
sae am | (x — &) ~b (y — n)? V (x — £)2? + (y — 9)? 


Substituting this back into Eq. (II-1), 


F 2u(é, y) Vix -—étP+e 
J, = lim I dx dy u(x, y) f dé - ) : [ + ‘] _ 
«> 0 € (x — &) 
(x — &) = l (x — &) 
“ih i+ — dn -| 1+ (II-5 
J->b ey" V (x — &)? + (y — 9)? ye y= V (x — £2? + (vy — nt 


This last equation is clearly symmetrical in y and y, since they only appear in the form (y — )?; the application oj 


—* 


the limiting process can have no effect upon this symmetry. 
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Diffuser Efficiency of Free-Jet Supersonic 
Wind Tunnels at Variable Test 


Chamber Pressure’ 


RUDOLF HERMANNT 
Unwersity of Minnesota 


SUMMARY 


The flow process occurring in the diffuser of a supersonic wind 
tunnel with free jet and its efficiency have been previously in- 
vestigated by the author (see Air Force Technical Report No. 
6334, December, 1950) for the special case, where the pressure in 
the test chamber surrounding the jet is equal to the pressure in 
the exit cross section of the Laval nozzle. 

The present analysis deals with the general case, where the test 
chamber pressure is different from the nozzle exit pressure. By 
using the basic equations of continuity, energy, and momentum 
flow through a control surface, an exact solution for one-dimen- 
sional, nonviscous, steady flow is obtained. 

One of the basic results is that the flow 
tion process in the diffuser is completed 
chamber pressure is above, equal to, or below the nozzle exit pres- 


after the transforma- 
is subsonic if the test 


sure. A supersonic flow is possible, too, if the test chamber 
pressure is smaller than a certain finite limit below the equilib- 
brium 

The results are presented in eight graphs over a Mach Number 
range from 1 to 10, showing velocity and pressure after the trans- 
formation, diffuser efficiency, total pressure ratio, and diffuser 
throat area at variable chamber pressure. The analysis also 
shows, for the first time, the exact mechanism of test chamber 
pressure control for the subsonic solution by means of an adjust- 
able diffuser throat and the quantitative description of a “‘slow 
starting” process of the tunnel. 


NOMENCLATURE 


a = velocity of sound 

f = ratio of nozzle exit cross-section area to diffuser intake 
area [Eq. (03) | 

f’ = combination of f and v |see Eq. (24)] 

F = cross-section area 

K = total pressure ratio 

! = length of free jet 

M* = w/a* = dimensionless velocity 

M = w/a = Mach Number 

pb = static pressure 

be = pressure in test chamber 

po = total pressure 
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po’ = total pressure after losses 

P = p/po = dimensionless pressure 

R = gasconstant 

$s = entropy 

T = absolute temperature 

w = velocity 

y = ratio of specific heats (1.405 for numerical calculations) 

€ = half of expansion angle from nozzle exit to diffuser 
intake (see Fig. 1) 

n = efficiency of pressure recovery |see Eq. (34) | 

v3 = function of /* or M [see Eq. (37) | 

© = dimensionless mass flux [see Eq. (13)] 

p = density 

v = dimensionless test chamber pressure difference [see 


Eq. (07)] 


Subscripts and Superscripts 


L = limit at which Kp, = 1 
max, = maximum 

a = supersonic solution 

B = subsonic solution 

NS = normal shock 

P = transformation process 
0 = stagnation 

l = nozzle exit 

2 = diffuser intake 

° = first throat 

*/’ = second throat 


(1) INTRODUCTION 


. SUPERSONIC WIND TUNNELS with regard 
to their type of test section, we can distinguish 
between those with closed wall test section and those 
with free-jet test section. The flow through the dif- 
fuser, the deceleration or transformation process, and, 
consequently, the pressure recovery are quite different 
in the case of the two tunnel types, as outlined before." ” 
When talking of pressure recovery or diffuser efficiency 
of a supersonic wind tunnel we must, therefore, differ- 
entiate clearly between the diffuser efficiency of a wind 
tunnel with closed test section and with a free jet. The 
term ‘‘diffuser efficiency”’ is used in the former’ * and 
in the present investigation as an abbreviation for 
“total pressure ratio of the diffuser divided by that of 
a normal shock at test section Mach Number.”’ A de- 
tailed reasoning for this definition was given previously 
(reference 2, pages 3-6), and its advantage will be evi- 
dent again during the following analysis. 

One of the basic features of a supersonic free jet is the 
fact that the static pressure surrounding the jet is, in 
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general, not equal to the static pressure in the Laval 
nozzle exit from which the jet emerges. In the case 
of a subsonic free jet, both of these pressures are neces- 
sarily equal. Hence, the supersonic problem has one 
variable more than the subsonic problem. The former 
papers! * calculated flow process and diffuser efficiency 
for the special case of pressure equilibrium. Then the 
pressure surrounding the free jet—i.e., the test chamber 
pressure—is equal to the nozzle exit pressure. For 
practical model testing this case is often preferred. 
The present paper deals with the general case of pres- 
sure nonequilibrium when the test chamber pressure is 
larger or smaller than the nozzle exit pressure. 

Even though pressure nonequilibrium may not be 
used for model testing, the results of this analysis are of 
importance in the design and operation of such tunnels. 
They clarify the flow mechanism of pressure control in 
the test section by means of the adjustable diffuser 
throat, necessary to obtain a desired test chamber pres- 
sure—for instance, pressure equilibrium. They furnish 
numerical data for the needed accuracy of such control 
device. Also from these results, the “‘slow-starting”’ 
process of this tunnel type, so important for all super- 
sonic wind tunnels, is now quantitatively understood. 

This analysis can also be applied directly to the de- 
sign of turbojet or ram-jet altitude test chambers if 
the engine jet is supersonic. In this case, pressure non- 
equilibrium must be used for testing in order to simu- 
late various flight conditions with respect to flight Mach 
Number, altitude, and fuel consumption. 


(2) ESTABLISHMENT OF THE FUNDAMENTAL EQUATIONS 


(A) The Transformation Process 


In establishing the fundamental equations, we refer 
to the procedure given previously" *in detail. The air 
(Fig. 1) comes from the atmosphere or a stilling cham- 
ber with stagnation pressure pp and stagnation tem- 
perature 7); passes the first throat F* of the Laval 
and attains supersonic velocity w,, pressure 
The air 


nozzle; 
pi, and density p; in the exit cross section F}. 
stream emerges from the nozzle exit as a free jet, which 
is surrounded by air at rest in the test chamber with 
the test chamber pressure p,. A certain length / of 
the free jet is required by the model length, its suspen- 
sion to the balance, or for the installation of other test 
devices. One of the main principles of the free jet is 
that it spreads out and increases its cross-section area 
in flow direction because of a turbulent mixing process 
at its surface. 

To obtain any pressure recovery, the free jet has to 
be brought into the diffuser. This requires a ‘‘recap- 
ture’ of the jet by means of a diffuser pickup in the 
intake area F2(F2 > F;) of the diffuser and a transforma- 
tion into a flow between rigid walls. We assume that 
this ‘“‘transformatior process’’ takes place in the intake 
portion of the diffuser with constant cross-section area 
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F;, which corresponds to a relatively long diffuser qe. 
sign with small slope of the walls. 
assume that in this ‘‘transformation zone” a complete 
transfer of all values of state occurs, which results jy 
constant velocity w2, constant pressure pe, and cop. 
stant density p2 all over the cross section Fy, After 
the transformation process is completed, the air flows 
through a convergent-divergent diffuser with adjust. 
able minimum section, through a fixed diffuser, and 
finally into the compressors. 


Furthermore, We 


(B) Assumptions and Fundamental Equations 


We consider the flow as one-dimensional, steady, and 
frictionless, and we use the three theorems of conser. 
vation of mass, axial momentum, and energy for a ree- 
tangular-shaped bounding surface, which extends from 
the nozzle exit to the end of the transformation zone 
and has a constant cross section F, (Fig. 1). The ad. 
vantage of such application of the momentum theorem 
is that the internal processes, such as recapturing and 
transformation, need not be known in detail. Evyi- 
dently, they are irreversible and connected with en- 
tropy increase or total pressure loss (py < 9), which 
will eventually be determined by’the calculation. 

The continuity equation states the equality of mass 
entering and leaving the control volumetT 


Fi piv = Fup, (01) 


The assumption of an airtight test chamber which 
has been made here is of greater importance than might 
be realized. Calculations show that a relatively small 
inleakage of air into the test section can produce quite 
a different flow process. 

The momentum equation states that the increase in 
momentum flow through the control surface equals the 
net force acting on the surface in the same direction, 


Frpxw.” — Fipow,? = Fipy + (F2 — Fi)p. 


i Fopo (05) 


An analysis shows that a neglect of the turbulent 
shearing stresses as forces at the control surface in this 
equation introduces only minor errors. Their main 
effect is the spreading of the jet, which requires that 
the diffuser intake area F, be larger than the nozzle 
exit area F;. This fact, however, is included in Eq. 
(05). 

The energy equation for a closed system states that 
the total energy entering the control volume must equal 
the outgoing total energy. With the aid of the equa- 
tion of state, this can be written 


WwW? y pr W@W y pe 


+ ° + —~- 


T The equations are numbered in agreement with the author's 


manuscript containing the complete calculations. 
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o@-AND RB-SOLUTION IN PRESSURE-VELOCITY * DIAGRAM 


The final states of the a and 8 solution in the dimensionless pressure (7) versus velocity (/*) diagram. Parameters 


are Mach Numbers /, and chamber pressure pv. 
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This equation is valid even under consideration of the 
work of the turbulent shearing stresses, because such 
work represents no exterior work. Since the air in the 
test chamber is in permanent energy exchange with 
the jet, it belongs to the system under consideration. 
The walls of the tunnel, including the test chamber, 
have to be insulated. 


(3) SOLUTION OF THE SYSTEM OF EQUATIONS 


The fundamental Eqs. (01), (05), and (11) are a 
system of three independent equations for the deter- 
mination of the three unknown variables—w», pe, p» of 
the final state. They can be calculated if the three 
variables w, 1, and p; of the initial state, the test cham- 
ber pressure p,, and the geometric boundaries F), F2 are 
given. We now have the general case p, ~ p; with one 
more variable than the former special case p, = /i. 

In order to simplify the set of equations, we introduce 
the area ratio f and the dimensionless test chamber 
pressure difference », 


f=F,/F. <1 (03) 
v = (pe — pr)/Pr (07) 
v = 0 corresponds to pressure equilibrium. v 2 cor- 


responds to test chamber pressure above or below the 
equilibrium. The minimum possible test chamber pres- 
sure is the absolute vacuum p, = Oor vy = —1. Cer- 
tainly, the jet will expand into the vacuum with con- 
siderable increase in cross section, and the one-dimen- 
sional analysis will not hold up to this mathematical 


limit. The maximum test chamber pressure _per- 
mitted will be discussed later [see Section (5B)]. For 
practical wind-tunnel purposes, » = +1 represents a 


reasonable positive limit that is physically possible 
above M = 1.362 [see Eq. (64) ]. 

The momentum Eq. (05) becomes, with Eqs. (01), 
(03), (07), 


f:Wip(W, — We) = po — p,[l + (1 — f)r] (10) 


We introduce the following dimensionless values for 
velocity, pressure, and mass flux: 


w p wp 
= = M*, = P, —. = (13) 
a Po a”‘p 


where P and 0 are well-known functions of 1/* and the 


entropy s {see reference 2, Eq. (5)]. Using 
a**p*/py = y[2/(y + 1)! (12) 
Eqs. (10) and (01) become 
P(M,*, Se) — P(M,*, 51) [1 4 (1 _ f)v| _ 
M,* — M,* 7 
» ¥/\4 1) 
/ as -f-O(M*, 51) (14a) 
7 ( 4 ) 1» 31 
0(M2*, s2)/0(Mi*, 1) = f (14b) 
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In order to eliminate s2, we further introduce the total 
pressure ratio of the transformation process 
Kp — po’ Po (15 
which is connected with the entropy increase by 
(So — 5S) ‘R = —In Kp (16 


The energy Eq. (11) gives, after some conversions and 

with Eqs. (15) and (16), 
P(M*, sc) = Kp-P(M*, s,) (17 

O(V/*, Ss.) = Kp:0(M"*, $1) (18 


We put Eqs. (17) and (18) in Eqs. (14a) and (14b), 
eliminate Ap, omit s;, and obtain 


Kp = f(90(M,*) /0(M.*)] (0 
P P 
/ , (M2") —f1l+(1 ~ 75 ta") = 


») Y/Y l 
~ -f-(M,* — M.*) (93 
. [. = ) i‘ 


Here and in the following, the P and 0 have to be taken 


at the entropy s;. Using the abbreviation 


f = ffi + 1 — foe] (24 
Eq. (23) becomes 
, P 
f’ — (M2*) — — (M*) = 
0 0 


> v/(y-1 
Y (. rt ) -f'(M,* — M2*) (25 
while in Eq. (20) the value f remains. Comparison of 
Eq. (25) with the corresponding equation of the former 
report [reference 2, Eq. (7)] gives the following results 
For v = 0, f’ = f and Eq. (25) becomes identical with 
former Eq. (7), which is a necessary agreement. Be 
yond that, former Eq. (7) becomes identical with Eq. 
(25) if f is replaced by f’. This fact is not trivial but 
allows a reduction of the solution for the velocity at 
y #'0 to that at »y = 0.3 Thus we obtain 


2 ly -—1 | 
7 + 1, vt+i1 I 
2M,* 


2 l — ] l 
Y (ie 2 
—{ 
\ 2M,* 


This equation is identical with the former velocity 
at pressure equilib- 





relation [reference | or 2, Eq. (10) 
rium if f is replaced by /’. However, this reduction is 
valid only for the velocities after the transformation. 
The pressure relation (20) contains f and cannot be re- 
We use subscript a for the positive sign and 
Exclud- 


duced. 
subscript 8 for the negative sign of the root. 
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; i es es oe _ 
ing special cases lor M,* = Land f’ = O, it follows from 


Bq. (29) 


M,,.*-M.,* = 1 (30) 
where 
M..*>1; Mat <1 (31) 
and both are always real and positive. 
Mathematically, our problem has two solutions—one 


supersonic a-solution and one subsonic £-solution. 
Their velocities are interconnected in the same way 
as the velocities before and behind a normal shock. 
These results are the same as those obtained in the 
case of pressure equilibrium. 

Analysis of the total pressure belonging to the a- and 
-solution will show if and under which conditions they 
will be permitted by the second law of thermodynam- 
ics. Eq. (20) furnishes the total pressure ratio of 
the a- and 6-solution 

Kp, = {(0(M,*)/0(M,,*) | 

Kpg = f (0(.M,*) /0(M,4*) | (33) 

The diffuser efficiency or pressure recovery factor is de- 
fined again by 

n = Kp/Kns(M,*) (34) 


where Ays(1/,*) is the total pressure ratio of the normal 


shock at the test section velocity 1/,*. Introducing 
for M* > 1 the function 
3(M*) = O0(M*)/Kys(M"*) (37) 
with the characteristic relation 
3(M*) = 0(1/M*) (38) 


and a finite limiting value for Mach Number infinity 
isee reference | or 2, Eq. (14c) |, the diffuser efficiency of 


the a- and 6-solution results in 


f{[3(M,*) /0(M,,*) | 


ae (39) 


ng = f[8(Mi*)/3(M,,.*) |] (42) 


In order to characterize the diffuser losses it is sufficient 
to use either Ap or n. For the a-solution, Ape is pre- 
ferred because the condition Kp, 1 is significant if 
the a solution is or is not thermodynamically allowed. 
For the 8-solution, y, is preferred because Kp, is, in 
general, smaller than | and for high Mach Numbers is 


inconveniently small. 


(4) DiscussION OF THE FINAL STATE OF THE a- AND 
8-SOLUTION IN THE PRESSURE-VELOCITY DIAGRAM 


General Remarks.—The discussion of the results in 
Sections (4) to (7) is made for the purpose of presenting 
this paper without further mathematical deductions 
simply by description of the graphs, obtained by nu- 
merical evaluation of the formulas above. They are 
calculated with the arbitrary value of the area ratio f = 
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0.64, which represents a good average for supersonic 
wind tunnels with a free-jet length between one and 
one-half and two times the nozzle exit height. It was 
used in the design of the two supersonic wind-tunnel 
test sections (40 by 40 cm.) at Peenemuende and 
Kochel, Germany—now in operation at Naval Ordnance 
Laboratory, White Oak, Md. However, the equations 
allow the calculations with any /-value. 

The final states of the a- and #-solution are repre- 
sented in the diagram of dimensionless pressure (P) 
versus velocity (/*) (Fig. 2). Each point represents 
The initial states lie on the 
The isen- 


one thermodynamic state. 
supersonic part of the upper isentropic line. 
tropic lines below the initial one designate higher en- 
states that are allowed by the second law of 
States above the initial isentropic 
Each final 


tropy 
thermodynamics. 
line are thermodynamically forbidden. 
state is designated by the intersection of two lines, one 
line of Mach Number ./, and one line of test chamber 
The net at the right-hand side (/,* > 1) 
the net at the left-hand side 
Each final a-state 


pressure vp. 
contains the a-solution; 
(\/.* < 1) contains the §-solution. 
point is connected with the corresponding final 8-state 
point by a normal shock process at the; velocity M,.* 
It might be emphasized, however, that ,the transition 
from one initial state to either one of the final states is 
never a normal shock, as this graph demonstrates. The 
diffuser transformation process is a different process of 
a more general nature than a normal shock. 

The régime of the 8-solution lies almost completely 
below the initial isentropic line. The few points at 
low Mach Number and v around +1, which lie above, 
are situated at the same time beyond the positive pres- 
sure limit yvys [see Section (5B)]. Hence, all 8-solu- 
tions within this limit are also permitted by the second 
law of thermodynamics. The régime of the a-solution 
extends to a large degree above the initial isentropic 
line (Kp, = 1) into the forbidden area. Therefore, 
only the a-solutions for low chamber pressure smaller 
than the negative pressure limit vz(./,) are permitted. 
The £-solution at this v; coincides with the lower por- 
tion of the well-known normal shock polar. The line 
for vy = —1 in the a-solution is identical with the isen- 
tropic line for Ap, = f = 0.64. In the 6-solution this 
line is affine to the normal shock polar by the factor /, 


which means nz = /, valid for all M,. 
(5) DrrFuSER EFFICIENCY OF THE 6-SOLUTION 


(A) Variation with Mach Number 


The diffuser efficiency ng of the 8-solution is shown in 
Fig. 3 as a function of the initial Mach Number with 
test chamber pressure v as the parameter. The general 
trend of the curves with different chamber pressure is 
the same as the curve for pressure equilibrium (v = 0), 
which at Mach Number | has the value ng = 0.888 and 


then decreases. All curves approach, for infinite Mach 
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DIFFUSER EFFICIENCY VS INITIAL MACHNUMBER 
Fic. 3. Diffuser efficiency n, of the 8 solution as a function of 


Mach Number J/, with chamber pressure v as the parameter. 
Number, the limiting value nz = f. Positive test cham- 
ber pressure improves the diffuser efficiency; negative 
chamber pressure diminishes the efficiency. At v = 
—1, ng is always equal to f. The influence of test 
chamber pressure on the efficiency decreases con- 
siderably with high Mach Number. 


(B) The Positive Pressure Limit vys of the Normal Shock 


The upper curve (Fig. 3) designated with vygs is the 
positive pressure limit, beyond which the #-solution 
has no physical meaning. It is the test chamber 
pressure at which a normal shock is established in the 
exit cross section F; of the Laval nozzle, followed by a 
free parallel jet with subsonic velocity and pressure 
equilibrium. vys is calculated from the theory of the 
normal shock, with /; equal to the static pressure be- 
fore and p, equal to the static pressure behind the shock 
at Mach Number /;._ This furnishes 


vys = [2y/(y + 1)] (Mi? — 1) (64) 


Obviously, this positive pressure limit is large at high 
Mach Numbers. 

Fig. 3 shows that the curve for ng(vys) starts at Mach 
Number | at the point for ng(v = 0), then increases, 
Using an incom- 
" 


and soon approaches a limiting value. 

pressible approximation for Apg(1/Mya,.*) given in 

the former analysis [reference 2, Eq. (26)] furnishes the 

formula 

; ~Xe -— Dd 

lim ng(yvs) = 1 — (1 — f)?- 
8 . 1)2 

Ms ‘7 + 4) 
This limit is always smaller than unity, and for f = 0.64 
Considering only the normal shock 


(70) 


its value is 0.987. 
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in the exit cross section would furnish y, = 1. The fact 
that 9,(vvs) 1s always below | is due to the transforma. 
tion losses of the following subsonic jet. 

A free-jet wind tunnel can be started “‘slowly”’ (quasj- 
steady) by means of gradually opening the diffuser 
throat [(see Section (7B)]. 
hereby reduced from a high positive pressure pv < y,, 


The chamber pressure js 


gradually down to pressure equilibrium v = 0 (Fig, 7), 
Hence, the diffuser efficiency decreases continuously 
(according to Fig. 3 or 4) down to the value for opera- 


tion. No jump occurs; no swallowing of the shock 
happens. No special reserve of compression ratio of 


the blowers for the starting process is required. In q 
tunnel of this type, the running compression ratio js 
higher than the compression ratio needed during the 


“é 


slow-starting”’ period. 


(C) Variation with Chamber Pressure 


Fig. 4 shows the same ng as in Fig. 3, but it is plotted 
as a function of chamber pressure and with Mach 
Number as the parameter. 
ing test chamber pressure. 
—1 and +1, the change in efficiency is nearly linear, 
Only the beginning of the positive pressure limit curve 
vys(M;) is shown in this diagram. The negative pres- 
sure limit vz(14,), below which the a-solution is per- 
mitted by the second law, is also plotted in this graph 
of the 6-solution. 


ng increases with increas- 
In the v range between 


(D) Linearized Approximation for High Mach Numbers 


An approximation of ng was obtained, which is linear 
in the term 1/M?. The calculations are carried out 
in the neighborhood of V/,,,;.*, linearized in terms of 


AM* = Miar.* — M*. First, Mo,* according to Eq. 


f*064 
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DIFFUSER EFFICIENCY VS. CHAMBER PRESSURE 
Fic. 4. Diffuser efficiency 1g of the 6 solution as a function of 
chamber pressure vy with Mach Number M as the parameter. 
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(29), is expanded as a function of M),9; * in terms of 
AM* for v = 0. Second, the functions of 0(M,*) and 
§(M..*) are expressed by 3(Mmar *) and are expanded 
in terms of AM*. Third, ng is calculated with those 
functions according to Eq. (42) valid for v = 0. Fourth, 


» ¥ Vis introduced by means of the relation between / 


Fifth, conversion is made 


and f’ according to Eq. (24). 
The result of the ex- 


from Mnar.* and AM* to 1/M,?. 
tended calculation is 


(Mh; i; y) = 


eo? 34 


, (97) 
2y My *' 


pit- Ja 4y). 
This equation shows the approach toward the limit- 
ing value nz = f with 1/M? at a given finite chamber 
pressure; the linear variation of ng with the test cham- 
ber pressure at a given Mach Number; and the value 
n3 = fat vy = —1 forall Mach Numbers. The accuracy 
of this equation is very good. For f = 0.64 and vy = 
(the deviation between Eq. (97) and the exact numerical 
values are 0.5 per cent at Mach Number 1.5 and smaller 
than 0.1 per cent from M = 2 up to infinity. Later, 
the partial derivative of Eq. (97) with respect to v is 

needed, 
Ong a + 1 l 


(Mi; f;») = U1 —-f) 


. (98) 
Ov : : 27 M;, 


The influence of the chamber pressure on the diffuser 
efficiency decreases with 1/M,*. For f = 0.64, the de- 
rivative becomes Eq. (98a) with an accuracy of about 
0.5 per cent from Mach Number 2 up, 


On,/Ov = 0.3081/M," (98a) 


(E) Stability of Test Chamber Pressure at the 8-Solution 
Eq. (98) shows that there is always 
Ons/Ov > 0 (99) 


The positive sign, which is not only valid in the linear- 
ized equation but in the total Mach Number range from 
1 to infinity and for all v between »,,;,. = —1 and 
Vmar. = vys, 1S Of great importance. We make a 
“quasi-steady”’ stability consideration as follows: A 
small positive pressure change in the test chamber 
(dv > 0) will cause an increase in efficiency ng and in 
total pressure behind the transformation zone. Hence, 
the mass flow determined by the Laval nozzle now has 
a smaller volume flow in the sonic diffuser throat F,’ 
[see Section (7A)] than before. Therefore, the throat 
is able to take additional mass flow which is drawn out 
of the test chamber, whereby the chamber pressure is 
reduced and the positive pressure disturbance disap- 
pears. A small negative pressure change in the cham- 
ber also vanishes in the corresponding manner. The 
opposite assumption of a negative sign in Eq. (99) 
would result in a completely unstable test chamber pres- 
Sure against disturbances. 
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TOTAL PRESSURE RATIO VS. CHAMBER PRESSURE 
Fic. 5. Total pressure ratio Kp, of the @ solution as a func- 


tion of chamber pressure » with Mach Number J/, as the param- 
eter. 


Obviously, the positive sign of Ong/Ov is the physical 
reason for the stable behavior of the test chamber 
pressure in supersonic free-jet wind tunnels as ob- 
For subsonic free-jet wind 
It is restricted 


served in actual operation. 
tunnels, such reasoning is not valid. 
completely to supersonic tunnels with a diffuser effi- 
ciency dependent on test chamber pressure. It 1s 
known that subsonic free-jet tunnels have a tendency 
for pulsation flow originating at the diffuser inlet, while 
in the supersonic free-jet tunnels at Aachen, Peene- 
muende, and Kochel, the author never observed any 


unstable flow phenomena. 
(6) ToTAL PRESSURE RATIO OF THE a-SOLUTION 


(A) Variation with Mach Number and Chamber Pressure 


In order to describe the transformation losses of the 
a-solution, we use the total pressure ratio Kp,. Fig. 
5 shows Kp, as a function of » with Mach Number as 
the parameter. At a given Mach Number positive 
chamber pressure improves Kp,, and negative chamber 
For v = —1, we have Kp, = 
For a given chamber pres- 


pressure diminishes Kp,. 
f = 0.64 independent of M. 
sure, Kp, decreases with Mach Number. The line 
Kp, = 1, identical with no change in entropy, divides 
the solution into two régimes. The lower régime with 
Kp < 1 means an entropy increase and is permitted 
by the second law of thermodynamics. The upper 
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TABLE | 
i 0 0.1 0.2 0.35 0.5 0.64 0.80 0.99 1 
VL« —1 —(). 834 —0.704 —0.540 —(0.396 —(0.274 —(.147 —0.0070 0 


régime with Ap, > | is connected with entropy decrease, 
and therefore the corresponding flow processes are 
forbidden. Evidently a limiting chamber 
pressure vz belongs to each Mach Number at and below 
which the a-solution is permitted. vy, ranges from 
vy, = —0.536 at Mach Number 1 to vy, = —0.275 at 
Mach Number 10. At pressure equilibrium, v = 0, a 
supersonic solution is not possible at all. This im- 
portant result was obtained formerly (see references 1 
and 2, Sections 4C, 2). 


certain 


(B) Linearized Approximation for High Mach Numbers 


At Mach Number infinity, the exact formula (32) 
for Kp, fails and becomes indeterminate. An expan- 
sion linearized in 1/M? is then helpful and is made in 
the same way as described for ng in Section (5D). The 
final formula gives the limiting value of Kp, for Mach 
Number infinity, 

lim Kp, (f, v) = 
M=o@ 
f 


Toa te oe Ae 120 
p—-@—-pannea-paeo 


The curve for J = © in Fig. 5 is calculated according 
to Eq. (120). The negative pressure limit v,.. for 
Mach Number infinity is now calculated by the re- 
quirement that Eq. (120) is equal to unity, which gives 


1 — f(y — f(y - 1) 
1—f 


Evaluation of Eq. (129) for various f gives the values 
for vz. aS shown in Table 1. 

For conventional wind-tunnel application, v,., has 
a finite, not vanishing, value. 
chamber pressure is needed until the supersonic solu- 
tion is permitted by the second law. 


VLo = (129) 


A definite negative 


(C) Regions of Existence of the a- and 8-Solution 


The total pressure ratio of both solutions Kp, , are 
The curves 
for Kp, are essentially the same as in Fig. 5. The line 
Kp, g — 1 separates the upper forbidden region from 
the lower region permitted by the second law. Also 
included are the boundaries for vys(14,), above which 
our system of fundamental equations becomes invalid. 
It is seen that for the a-solution the function vz(M)) is 
the effective upper v-limit, while for the 6-family the 
upper v-limit is given by the function vys(M,). The 
following results are obtained: 

(1) The £-solution is permitted for all negative 
chamber pressures, for pressure equilibrium, and for 
positive chamber pressure up to vys(W), where the 
normal shock is established in the nozzle exit. 


shown in Fig. 6 with logarithmic ordinate. 


(2) The a-solution is permitted for low negative 
chamber pressure with the upper limit »,(M,) < 9 
The a-solution is forbidden for higher negative pres- 
sure above v,(./;), for pressure equilibrium, and for aj] 
positive chamber pressures. 

(3) Thus, in the negative chamber pressure range 
at and below the limit »;(1/,), both the a- and the 3. 
solution are permitted. Which one will be established 
actually cannot be decided by this analysis alone. |t 
requires a study of the continuation of the a-solution 
through the diffuser under consideration of friction. 
The experimental data [see Section (7B)] can be 
understood qualitatively. 

Fig. 6 may suggest the opinion that it is desirable 
to establish the a-solution, if possible, because of its 
much higher total pressure ratio. However, the flow 
through the total diffuser must be considered [see 
Section (7A)] up to the region, where the supersonic 
velocity is finally converted into low subsonic velocity, 
Kpg is small but represents practically all losses, while 
Kp, is higher but does not include the additional large 
losses downstream of the transformation zone with the 
Mach Number /,, > M,. 


f=0.64 Me | S 
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Fic. 6. Total pressure ratio Kp, of the @ solution as a fune- 


tion of chamber pressure and that of the 8 solution Apg as 4 


function of chamber pressure v with Mach Number 1, as the 
parameter. Regions of their existence. 
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DIFFUSER THROAT CROSS SECTION OF THE 


(7) 
§-SOLUTION AND TEST CHAMBER PRESSURE CONTROL 


(A) Conditions of the Second Throat 


In order to calculate the diffuser throat and discuss 
pressure control by means of the adjustable diffuser 
we have to follow the flow further downstream of the 
transformation zone through the diffuser. In case of the 
g-solution, the subsonic flow in the entrance region of 
the diffuser will be accelerated within itsconvergent part, 
connected with small boundary layers and negligible 
losses by friction. If the pressure downstream of the 
diffuser is low enough, the flow attains sonic velocity in 
the minimum cross section and supersonic velocity in the 
divergent part. (Subsonic velocity is finally obtained 
by anormal shock whose location can be adjusted by the 
compression ratio of the blowers close downstream of the 
minimum section, hence encountering negligible losses.) 
The diffuser acts as an accelerating conventional Laval 
nozzle, not as a decelerating reverse nozzle as in the 
case of a closed test section. The mass flow passing 
the second throat with sonic velocity must be equal to 
the mass flow entering through the Laval nozzle, when 
steady state is present, furnishing the equation 

wd 3(M,*) 


“_ (M,*, f, v) = 


= (105) 
F, : ng(M,*, f, v) 


If M,* and f are given, a certain chamber pressure 
vy causes a certain efficiency | Eq. (42)] and thus requires 
in a steady motion a certain throat area according to 
Eq. (105). If now the diffuser throat is made smaller 
than given by Eq. (105), for instance, by voluntary 
throat control, the mass flow coming from the Laval 
nozzle is not able to flow completely through the throat. 
The difference in mass flow must accumulate in the 
test chamber, raising the test chamber pressure, which 
increases the diffuser efficiency. v and ng rise until 
Eq. (105) is fulfilled. When the diffuser throat is made 
larger than given by Eq. (105), the opposite is valid. 

In case of the a-solution, the convergent diffuser part 
causes a deceleration of the supersonic flow with pres- 
This will produce strong increase in bound- 
oblique 


sure rise. 
separation or 
with considerable 


ary-layer thickness, maybe 
shocks, all processes connected 
losses. Evidently, one-dimensional frictionless analy- 
sis is not applicable for the continuation of the a-solu- 
tion through the diffuser. According to present experi- 
ence, the velocity can be decelerated only partly, and 
no sonic velocity is attained in the minimum area. 
Thus it does not act as a second throat, and chamber 
pressure control by means of the adjustable diffuser is 


not possible. 
‘B) Variation with Mach Number and Chamber Pressure; 
Operation of Pressure Control 


Fig. 7 shows the second throat area of the 6-solution 
, *I DD . : . y . 
F;*'/F, as a function of Mach Number with the chan- 
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as f=0.64 
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SECOND THROAT AREA VS. MACH NUMBER 
B-SOLUTION 
Fic. 7. Relative second throat Fg*’/F;, of the 8 solution as a 
function of Mach Number MM, with the chamber pressure v as the 


parameter. 
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values, linearized approximation, and simplified linearized ap- 
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ber pressure v as the parameter. It also gives the com- 

plete curve for pressure equilibrium through the sub- 
sonic region down to Mach Number zero (taken from 
references | and 2) and the curve for vys. The parts 
of the v curves below curve vys have no physical mean- 
ing. 

The test chamber pressure control operates the fol- 
lowing way: Assuming a fixed Mach Number we may 
start with a diffuser throat belonging to the positive 
pressure limit vys, which produces a normal shock in the 
nozzle exit. Now we open up the diffuser throat gradu- 
ally and pass one p-line after the other, designated with 
decreasing but still positive chamber pressure.  Fi- 
nally, we reach pressure equilibrium at the line v = 0. 
Further enlargement of the diffuser throat produces 
decreasing negative chamber pressure until the nega- 
tive pressure limit vz, is reached. In this whole range, 
only the §-solution exists and the chamber pressure con- 
trol works perfectly. This control is important for 
actual model testing and has been used for about 20 
years in supersonic wind tunnel practice, but it never 
before was understood. At v, we have the alternative 
[see Section (6C)] of the -solution with continued 
chamber pressure control or of the a-solution, where 
chamber pressure control by the adjustable diffuser 
According to the author's 
former experience in operation of free-jet supersonic 


throat no longer exists. 


wind tunnels, the latter will probably occur, at least at 
Mach Numbers above the order of 2 to 3. 


(C) Chamber Pressure Control Term—Linearized 


Approximation 

For some purposes of design and operation of wind 
tunnels, it is valuable to have quantitative data of 
the pressure control which are more explicit than those 
given in Fig. 7. 
throat area with the chamber pressure, called the 


Therefore, the rate of change of the 


“chamber pressure control term,’’ was calculated with 
the aid of Eq. (105) and is represented in Fig. 8 for f = 
0.64, designated as “exact values.’’ It is sufficient to 
present the control term for v = O only because the 
influence of v is extremely small as long as we stay 
within vy = +1 for wind-tunnel purposes. Obviously, 
the sensitivity of pressure control increases consider- 
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ably with high Mach Numbers. A linearized approxi- 
mation is obtained by using ng (97) and On, / dv (98) 


l -_ 7 
L-fvt! 7 l-fy+1 1 
f..2y Me f 24 a (i 


Fig. 8 shows Eq. (114) as ‘‘Approximation.”’ If we set 
the bracket equal to unity, we obtain the ‘simplified 
approximation.’’ The first equation is exact within 
3 per cent and better from Mach Number | to 10. The 
simplified equation is exact within 2 per cent and better 
from Mach Number 4.5. Eq. (114) allows fast calcula. 
tions of the chamber pressure control term for any ; 
value. 


CONCLUSIONS 


The flow process in the diffuser and its efficiency are 
investigated for the general case, when the pressure in 
the test chamber surrounding the free jet is different 
from the pressure in the exit cross section of the Laval 
nozzle. An exact solution for one-dimensional non- 
viscous steady flow is obtained. Velocity and pressure 
after the transformation process, total pressure ratio, 
diffuser efficiency, and diffuser throat area are calcu- 
lated for variable chamber pressure. The exact mech- 
anism of test chamber pressure control by means of 
the adjustable diffuser throat is shown. Other results 
are the proof of the stability of the test chamber pres 
sure, the existence of two different types of flow at low 
negative chamber pressure, and the quantitative de- 
scription of a ‘‘slow-starting’’ process of the tunnel. 
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SUMMARY 


Using the ballistic pendulum air-drag technique with pendu- 
lums of several different lengths, the air-drag coefficient, Cp, has 
been determined for cal. .50 spheres in the region of Mach 
Numbers 1.4 to 5.4. A conventionally suspended pendulum with 
a period of 3.2 sec. was used to obtain most of the data; the values 
of Cp so obtained were consistently about 6 per cent low. It was 
suspected that the time of momentum transfer covered too great 
a fraction of the pendulum period. For this reason, pendulums 
having seismographic suspensions with periods from 17.5 to 30 
sec. were installed. After correction for end effects, the values 
of Cp then obtained agreed closely with the results of other meth- 
ods. The influence of baffle-hole size and correction for end ef- 


fects are discussed. 








SYMBOLS 


a = air-drag coefficient uncorrected for end effect 
p = air density 
o = 2x/P 
A = area of projectile normal to trajectory 
Cp = air-drag coefficient corrected for end effect 
D, = base length between first and second screen 
D, = base length between second and third screen 
E = (1/2)mV? kinetic energy of projectile (no spin) 
AL = apparent additional length of pendulum, representing 
end effect 
L = physical length of the pendulum 
m = mass of projectile 
P= period of the pendulum 
S = displacement of projectile in direction of trajectory 
7, = transit time over D, 
7, = transit time over D, 
V = projectile velocity 
Ve = velocity of projectile at pendulum exit 
Vi = projectile velocity loss to pendulum over path L + AL 
V; = velocity over first base length 
6 velocity over second base length 


X = pAS/m 
X’ = pAL/m 
X”" = pA(L+AL)/m 


INTRODUCTION 
Origin of the Ballistic Pendulum Air-Drag Technique 
| (aio TO STUDIES OF ARMOR for resistance to 
fragment penetration, questions arose during the 
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summer and fall of 1944 as to the reliability of estimates 
of remaining velocity of fragments. The desire to check 
fragment drag coefficients in current use coincided with 
completion of a new high-speed spark chronograph, 
which gave sufficient precision to permit the measure- 
ments being made on an improvised 20-ft. range in a 
basement laboratory The results of these 
chronograph drag measurements were satisfying with 
This brief 


room. 


respect to the immediate issues in question. 
activity in ballistic aerodynamics preceded and helped 
stimulate the conception of making laboratory meas- 
urements of air drag with a ballistic pendulum. 


In 1943, personnel from the Naval Ordnance Labo- 
ratory sought assistance from the Naval Research Lab- 
oratory in their studies of mechanical fuze develop- 
ment. Dr.G.R. Irwin, Superintendent of the Mechan- 
ics Division and Acting Head of the Ballistics Branch 
of the Naval Research Laboratory, proposed that the 
momentum loss in plates be precisely determined using 
a special ballistic pendulum technique. This tech- 
nique consisted of mounting a plate in a_ ballistic 
pendulum and determining the momentum lost by the 
projectile to the plate by computing, from the pendulum 
deflection, the momentum gained by the plate. The 
success of this method led Dr. Irwin to believe that 
momentum lost by a projectile to air could be meas- 
ured in a similar manner. 


In March, 1945, at the Bureau of Ordnance, Dr. 
Irwin discussed with Lieutenants Hunt and Hoffman 
the possibility of measuring air drag with a ballistic 
pendulum. It was agreed that the method was prom- 
ising, and active experimental work was initiated in 
August, 1945, using °/;s-in. steel spheres in a cal. .30 
In August, 1947, pendulums were 
All the 
At- 


smooth bore gun. 
installed for work with cal. .50 projectiles. 
cal. .50 work has been concerned with spheres. 

tention is directed to the advantages of this system: 


(1) The pendulum method measures momentum 
rather than velocity loss and, for equal precision, uti- 
lizes a much shorter trajectory than any of the known 
chronograph techniques for determining resistance of 
gases to penetration by projectiles. 

(2) About 1 per cent precision is obtainable readily 
enough so that the experiment can remain simple and 
convenient. 
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Fic. 1. Conventional cradle pendulum support. 








CEILING SUSPENSION POINTS 


8 = ANGULAR DISPLACEMENT OF CEILING 
SUSPENSION POINTS FROM VERTICAL 











CONE POINT PIVOTS 

Fic. 2. Seismographic-type pendulum showing the angle @ be 
tween the line from ceiling suspension points to cone pivots and 
the vertical. The period p = V//g sin 6, where / is the length of 
the arms. 
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LONGITUDINAL SECTION OF PENDULUM CELOTEX BAFFLE 


Fic. 3. Arrangement of baffle discs. 


(8) The transonic and subsonic ranges are not too 
low for use of the pendulum air-drag technique. In 
contrast to chronograph methods, accuracy will improve 
with increasing velocity of the test missile. 

(4) Missiles may be of any shape that can be fired 
in proper orientation through a 100-cal.* path in air. 

(5) The pendulum method is adaptable for study 
of such features as head shape, the effect of yaw upon 


* That is, 100 projectile diameters 
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drag, scale effect, and the effects of various physical 
conditions of the resisting medium. Adaptations of 
the method to the study of lift are possible but will re. 
quire a considerable amount of experimental develop. 
ment. 


EXPERIMENTAL PROCEDURE 


Measurement Technique.—When a missile is fired 
through a hollow, partly closed pendulum, the pendy- 
lum acts as a trap for the momentum imparted to the 
air. The swing of the pendulum increases in a known 
way with its momentum gain. In a way not fully under- 
stood (since details of momentum transfer from air to 
pendulum are still obscure), the pendulum deflection 
There 
are three principal trajectory sections: (1) the path of 
the missile which creates pressure differential against 


increases with momentum loss of the missile. 


the entrance face of the pendulum; (2) the interior 
and (3) the path that may cause transfer oj 


To meas- 


path; 
momentum through the exit-face opening. 
ure the air-drag coefficient successfully, it is desirable 
to design the experiment so that sections (1) and (3) are 
independent of pendulum length, so that the wall or 
choking effects along section (2) are negligible, and so 
that the air motion inside the pendulum becomes ran- 
dom in a time less than one-thirtieth of a quarter period 
of the pendulum. When these requirements are satis- 
fied, firings at the same measured velocity with several 
pendulum lengths will serve to determine the end cor- 
rection and will, therefore, permit absolute determi- 
nation of the air-drag coefficient. 

The pendulums were constructed so that the en- 
trance and exit faces were similar. All firings were 
made with spheres; choice of this shape eliminated 
the yaw problem and reduced expense, since commer- 
cial ball bearings were used. 

Equipment.—Outside the pendulum room: A cal. .50 
sphere was given a horizontal velocity by firing through 
The amount of gun powder used 
provided the desired velocities. After leaving the gun, 
the sphere passed through a steel cylindrical silencer 3 
ft. in diameter and 10 in. long. Attached to the si- 
lencer on the side facing the pendulum was a blast de- 
flector 1'/, ft. long with 14 square aluminum plates 2 in. 
apart; each plate measured 1 by | ft. and had a 1-in. 
diameter hole for the passage of the sphere. 


a cal. .50 gunbarrel. 


Inside the pendulum room: Pendulum bodies, whether 
of the conventional Fig. 1 or seismographic-type Fig. 2 
suspension, were identical in construction. Aluminum, 
0.032 in. thick, was formed into a cylinder 10 in. in 
diameter and of appropriate length. Art board, '/s in., 
was cut into discs to fit within the pendulum; each 
disc had at the center a hole 4 or 5 cal. in diameter. 
Starting at one end, the discs were placed in the cylinder 
spaced by aluminum frameworks, which served to 
separate the discs by 4 in. These discs act as bafiles 
as shown in Fig. 3. 
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FREE-FLIGHT AIR-DRAG 


Qn the outer surface of the cylinder, an aluminum 
projection 1/,by 1 by 1 in. was fastened so that a normal 
off the 1- by 1-in. face was parallel to the direction of 
pendulum travel. This projection, called a ‘‘pusher,”’ 
served to displace a free pencil lead which rested hori- 
zontally in two V supports on a stationary block. A 
traveling microscope was mounted on a solid stationary 
support in such a manner that the distance between the 
pusher at the pendulum rest position and the end of 
the pencil lead after displacement could be determined. 
The sphere passed through the pendulum without direct 
contact, and the air disturbed by the projectile and its 
shock wave caused the pendulum displacement. 

After the rest position of the pendulum is deter- 
mined, the pendulum must be held there until the in- 
stant before firing. Restraint was effected by a damper 
actuated by a solenoid and controlled by push buttons 
either from the pendulum room or from the firing point. 

The period of the pendulum was determined with a 
stop watch. In the formula F = — Mw*d (where w = 
2x/P and P = the period of the pendulum), the force, 
F, and the displacement, d, were directly measured by 
holding back the pendulum by a weight over a low- 
friction pulley. Thus the mass, M, of the pendulum 
was determined. To support the 176-in. pendulum 
two cradles were used, and three cradles were needed 
for the 273-in. pendulum. 

The velocity of the projectile at the center of the 
pendulum was measured by a 1.6-mce. electronic counter, 
which is operated from pulses fed to it by two photo- 
electric screens placed about 2 ft. from each end of the 
pendulum. The exit velocity is then approximately 
V, — (1/2) Vz, where V, is the velocity at the center of 
the pendulum, computed from the counter observation ; 
l’, is the velocity lost by the projectile through the 
length of the pendulum, computed from the observed 


deflection. 


THEORY 


Drag Calculations 


Let us define the air-drag coefficient, Cp, by the re- 


lationship 
—dk/dX = CpnkE (1) 
whence 
E => Eo exp (—CpX) (2) 
and 
V = Vo exp [(—CpX)/2] (3) 
where 
E = (1/2)mV?, kinetic energy of the projectile (no 
spin) 
V = projectile velocity 
m = mass of projectile 
X = pAS/m 


MEASUREMENT 
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p = air density 
A = area of projectile normal to trajectory 
S = displacement of projectile in direction of tra- 
jectory 
If S = L, 
so AL 
x’ = (4) 
m 
andif S= L + AL, 
X" = pA(L + AL)/m (5) 


where L = physical length of the pendulum and AL = 
apparent additional length of pendulum, representing 
the end effect. 

For small increments in X and S, the relation may be 


written 
J h (CpX ) ( aX ") . 
1+ 7 = exp : = exp : (6) 
BK - - 
where 
V;, = projectile velocity loss to pendulum over path 
L+AL 
Ve = velocity of projectile at pendulum exit 
a = air-drag coefficient uncorrected for end effect 
Cp = air-drag coefficient corrected for end effect 


Eq. (6) is convenient to use since I’, and I’, may be 
readily computed, X’ is known, and aX’/2 can be ob- 
tained from a table of the exponential function when 
the value of 1 + V,/V,is known. The relation 


Cp/a = L/(L + AL) (7) 


follows directly from Eqs. (4), (5), and (6). 
Method of Determining Cp 


If ais plotted against the reciprocal of the pendulum 
length, one can extrapolate to Cp at 1/L = 0. Fig. 4 
is plotted in this manner. It is interesting to note that, 
if this curve is extrapolated so as to intercept the hori- 
zontal axis, this intercept will be the reciprocal of the 
end correction, AL. 

The equation 


Cp = (Lia, om Leas) (L, — Le) (S) 


where L, and Le are any two pendulum lengths and 
a, and a: are the corresponding values of a, may be 
more convenient for the determination of Cp than the 
graphical method. The validity of this equation may 
be shown at once by plotting a vs. L for several values 
of L. The slope of the resulting curve is Cp, the hori- 
zontal intercept is AL, and the vertical intercept is 
aAL. 


AIR DRAG OF CAL. .50 STEEL SPHERES BY A 


CONVENTIONALLY SUSPENDED BALLISTIC PENDULUM 


The values of a were investigated to determine the 
effect of baffle-hole size and pendulum length. As was 
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against the reciprocal of the pendulum length. Data for 4 and 5 
cal. baffle hole diameters is plotted. Graph (b) is the same except 
that the a values were obtained at Mach = 2.0. For comparison, 
the values read from graph (a) and (b) at 1/L = 0 are shown on 
Fig. 5d and are indicated at 4A and 4B. The values were ob- 
tained using a conventionally suspended pendulum. 


expected, the larger the baffle hole, the smaller the a 
for the same Mach Number. Curves of a vs. Mach 
Number for the 5-cal. baffle-hole diameter were made 
for pendulum lengths of 80, 176.6, and 273.3 in., as 
shown in Fig. 5.* Data were taken for a 4-cal. baffle- 
hole diameter at Mach Numbers 2 and 3, using S0- and 
176.6-in. pendulums. Values of Cp, obtained by plot- 
ting the values for the 4- and 5-cal. baffle-hole diameters 
against 1/Z and extrapolating to 1/L = 0, were the 
same (Fig. 4) and shown in Fig. 5. 


AIR DRAG OF CAL. .50 ALUMINUM SPHERES BY A 
VELOCITY- DIFFERENCE METHOD 


To have data with which to compare results of the 
ballistic pendulum air-drag technique, Cp was deter- 
mined by the velocity-difference method. Cal. .50 
aluminum spheres were used, since their low density 
made their velocity loss over a given distance almost 
three times that of the steel spheres. Screens were 
set up 25 ft. apart to give two nearly equal base 
lengths. 

In the calculation of Cp, no consideration has been 
given to the difference in base lengths between screens. 
As shown in the Appendix, if the experiment is re- 

* A, B, and C, respectively. In Fig. 5d these curves are com- 


pared with the values of Cp obtained by Charters and Thomas 
and the values of Cp from Fig. 4A and 4B. 
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peated with the direction of firing reversed, the ayer. 
age Cp calculated for the two directions will be identical 
with the average Cp calculated with base-length dif. 
ferences taken into: consideration. That is, averaged 
Cy values obtained by this method, ignoring possible 
errors in the location of the center screen, will not be in 
error.{ The values of Cp obtained by this method 
are quoted in the report and are averages of this kind. 
The velocity of the projectile was measured over each 
base length, and values of Cp were computed using the 
equation 

Vi/ V2 = exp (CpX/2) = T2/T, (9 
where |’, and V2 are the velocities of the projectile over 
the first base length D, and the second base length D.,, 
respectively, and where 7 and 7» are the correspond- 


+ Actually, this statement is correct only insofar as the com- 
putation method outlined earlier in the report is used. A more 
exact treatment shows that this averaged Cp value will be cor 
rect within second-order terms in (D,; — Dz2)/(D,; + Dz). See 


Appendix for further discussion. 





















































(a) 
a T T 
' } | | r 
L.Or— co’ PENDULUR TT t 4 
i=] | } | 
| | | 
0.9 } 1 i | | 1 1 
(b) 
it , T T T 
Or 
i] 
os} 
(c) 
it : T T T T T T 
Lo} —4 hee A OO et ae cee 
& | | 
09 1 
(d) 
'\("g0" PENDULUM 
762 
1.0 PENDULUM —4——— - 
SC saw’ . 
2733" PENDULUM Das . 
09) on + + + + + + + —{ faa] _ + 
4A | | 
| 12 1@ 16 18 2 22 24 26 286 3 32 34 36 38 
MACH NUMBER 
Fic. 5. a@ vs. Mach Number for cal. .50 steel spheres, 5-cal 
2'/2-in.) baffle hole, conventional pendulum suspension. x 


represents values obtained by Charters and Thomas in (d 
4A and 4B are values of Cp obtained from Fig. 4a and 4b. 
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Cp from N.R.L. free-flight velocity difference method 
compared with Charters and Thomas values. 
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FREE-FLIGHT AIR-DRAG 


ing transit times. Again, time was measured by Potter 
1,6-me. electronic counters. Photocell light screens 
were used for triggering. Data obtained in this man- 
ner are plotted in Fig. 6, together with values obtained 
by A. C. Charters and R. N. Thomas* for °/,6-in. steel 
spheres. It will be noted that the values of Cp ob- 
tained by these two methods agree closely. 


DISCREPANCY IN Cp VALUES OBTAINED BY DIFFERENT 
METHODS 


Values of Cp obtained by extrapolation after plotting 
avs. 1/L are given in Fig. 5d, together with values ob- 
tained by Charters and Thomas, as well as curves ob- 
tained by the ballistic pendulum method at N.R.L. It 
will be noted that the values of Cp corrected in this 
way are about 6 per cent lower than those reported by 


Charters and Thomas. 


Probable Cause of Lower Cp Values 

Since the period of the conventionally suspended 
pendulum was 3.2 sec., it seemed probable that the time 
of momentum transfer was long enough to affect the 
deflection reading. An analog of the pendulum was 
set up using a ballistic galvanometer, a charged con- 
denser, and a resistor. The RC product of the con- 
denser and the resistor, plus the resistance of the 
galvanometer, was used as a measure of the time re- 
quired for the charge to be transferred through the 
coil of the galvanometer. The capacity of the con- 
denser and the voltage applied were adjusted to ob- 
tain a nearly full-scale reading with no external resist- 
ance, 

The period of the galvanometer was 8.6. sec. 
Fig. 7 shows the per cent of maximum deflection 
plotted against the time expressed as the per cent of the 
total period of the galvanometer. 


From Fig. 7, 94 per cent of the maximum deflection 
corresponds to having an RC product that is 4 per cent 
of the galvanometer period. Analogously, with the 
period of the conventionally suspended pendulum equal 
to 3.2 sec., the momentum reading would be 94 per cent 
of maximum if the time of momentum transfer were 
4 per cent of 3.2 sec. or 0.128 sec. If a 99.5 per cent 
momentum reading is desired, the time of momentum 
transfer, judging from Fig. 7, should be about 0.7 per 
cent of the period of the pendulum or about one-sixth 
of the time that gave 94 per cent of the maximum de- 
flection. 


Since increasing the period of the pendulum would 
have the same effect as decreasing the time it 
takes for the momentum transfer, increasing the 
pendulum period by a factor of 6 should give values of 
Cy which agree closely with established ones. 


* Aberdeen Proving Ground, Ballistic Research Laboratory 
Report 514; Charters, A. C., and Thomas, R. N., JOURNAL OF THE 
AERONAUTICAL SCIENCES, Vol. 12, p. 468, October, 1945 
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UsE OF SEISMOGRAPHIC-TYPE SUSPENSIONS TO OBTAIN 
LONGER PERIOD PENDULUMS 


Seismographic-type suspensions? (Fig. 2) were tried 
using pendulums with periods of 17 to 30 sec. The 
values of Cp thus obtained agree closely with those ob- 
tained by Charters and Thomas and with values ob- 
tained at N.R.L. by taking velocity differences. For 
comparison Fig. 8 shows these Cp values, as well as 
those obtained by Charters and Thomas. Good agree- 
ment will be noted, except at Mach Number = 1.6 
where the pendulum values are about 2.5 per cent low. 


CONCLUSION 


The ballistic pendulum method of determining the 
air-drag coefficient is particularly applicable to ex- 
tremely high Mach Numbers, provided the arrange- 
ments include suitably reduced baffle-hole sizes and a 


sufficient reduction of muzzle blast. Since the base 


+ Use of seismographic-type suspensions for this work was sug- 
gested to the writers and given initial trials at N.O.L. by Dr. 


George Shue. 
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TIME-RECORDING STATIONS 


INITIAL DIRECTION 





OF FIRING TRAVERSE 
times (Oty Mts 
| 
= _ple — 
Fic.9 The experimental arrangement. 


length over which the projectile velocity is determined 
is large and since the deflection of the pendulum in- 
creases with projectile velocity, no decrease of precision 
in the determination of projectile velocity or momen- 
tum loss is noted at high Mach Numbers. Because 
the total length of trajectory need be only 30 ft., some- 
what higher Mach Numbers can be maintained than 
are possible with the customary spark chronograph 
ranges. 

This method of determining the air-drag coefficient 
can be used with missiles of any shape which can be 
projected through the pendulum, provided the yaw and 
drift, if present, can be controlled to prevent contact 
between the missile and the pendulum. The conven- 
ience of this method makes it adaptable for use in 
universities either for the determination of the air-drag 
coefficient or for instructing students in aerodynamics. 


APPENDIX 


By Robert E. Roberson 


Referring to Fig. 9 for experimental arrangement, 
let D,; and D. denote distances between screens, defin- 
ing D = (1/2) (D, + Dz), and let At, and Afs be the 
transit times between screens of a projectile fired in the 
direction indicated by the arrow. It is desired to find 
the error made by calculating Cp as if Dj = D. and to 
show how it is reduced by reversing the direction of 
firing and using the average value of Cp calculated for 
the two directions. 

Let us represent by Cp the value of the drag co- 
efficient calculated recognizing the difference between 
D, and D, and by Cp* the approximation to this value 
calculated upon the assumption that D; = D.. Postu- 
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late that (1) the value Cy is calculated from the equa- 
tiont 

2m D, Ate 


Cp = 1 
o ghd” Dal 


(10) 


and (2) that, except for the errors in the approximation 
Cp* = Cp, the errors of the experiment are the same for 
the two directions of firing. 

The first of these postulates makes it unnecessary to 
examine in this Appendix the rationale of the calcula- 
tion procedure. The second allows the errors in the 
location of the center screen to be considered aside from 
the various other experimental errors. 

An additional notation, let subscripts f and 7 repre- 
sent forward and reversed directions of firing and a bar 
above any quantity represent its average value for the 


fand,r experiments. 


From the definition of Cp and Eq. (10), there result 


immediately 

C 2m Dy, Ate s 
(Cp)r n (11) 

- pA D De Aty 

. 2m Ate 
(Cp*), = In (12) 

pA D Af; 

: ; 2m dD, 

(Cp); — (Cn*), = >AD n D. (13) 


the roles of D,; and D, 
Therefore, by analogy 


Clearly, in the reversed firing, 
are effectively interchanged. 
with Eq. (13), 

2m Ds 


(Cp), — (Cp*), = n 
D D ) pAD D, 


(14) 


Taking one-half the sum of Eqs. (13) and (14), one has 


Cp sed Cp" ( 15) 
That is, if the Cp value of the experiment is the average 
of the Cp values obtained by firing through the screens 
in two directions, it contains no error by virtue of being 
calculated on the assumption that D; = Dz. 


argument, con- 


A somewhat more sophisticated 
structed upon a less restrictive base than postulate 
(1), shows that actually the error made in Cp by con- 
sidering D, = Dz is of second order in (D,; — Dz2)/- 


(D, + Ds). 


+ The listing of the more basic assumptions that lead to this 
formula is beyond the scope of this Appendix. 
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Impingement of Spherical Water Droplets 
on a Wedge at Supersonic Speeds in Air 


MYRON TRIBUS* ann ARMAND GUIBERT? 
Unwersity of Michigan and Unwersity of California at Los Angeles 


SUMMARY 
rhe water-drop interception characteristics of a wedge in a 
supersonic flow field are investigated analytically for the case of 
flows in which the shock wave is not detached. The results indi- 
cate that, contrary to the subsonic flow cases considered previ- 
ously, there is no limiting position on the wedge beyond which the 
droplets cease to strike, but rather icing may occur along the 


entire wedge surface. 


INTRODUCTION 


Peconic and experimentally determined icing 
characteristics for various objects have been pre- 
viously reported in the literature.'~* All of the pre- 
vious reports have been concerned with subsonic flow 
at velocities where the compressibility of the fluid is not 
important. Because the present trend in aircraft de- 
sign is toward supersonic flight speeds, it is of interest 
to consider the icing characteristics that will occur at 
these higher velocities. At these very high velocities, 
small changes in wing profile, for example, are expected 
to produce much greater changes in flying character- 
istics than occur at subsonic speeds. 

Although at very high flight speeds it is certain that 
the effects of aerodynamic heating are sufficient to pre- 
vent icing, there still exists a combination of low super- 
sonic speeds and low ambient temperatures at which 
icing may occur. This icing may occur even though the 
stagnation temperature in dry air is above freezing, be- 
cause the presence of the moisture on the wing will cause 
a depression in temperature (similar to wet bulb de- 
pression) with the result that the wing may fall below 
the freezing temperature. This effect has been dis- 
cussed in a previous report.* 

Even when no ice forms, it is of interest to compute 
the rate at which water droplets strike, since the results 
may be of value in erosion studies. 

As an initial contribution to the solution of the over- 
all problem, there is presented an analysis of the water 
interception characteristics of a wedge in a supersonic 
flow field. 


NOMENCLATURE 


a = acoustic velocity ({t. per sec. ) 
a = acceleration vector (ft. per sec.?) 
b = distance from centerline of wedge (ft.) 


Received April 16, 1951. 
* Director of Icing Research. 
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Cp = drag coefficient (dimensionless 
F = force vector (Ibs 
m = mass (slugs 
MJ = Mach modulus (dimensionless 
P = pressure (Ibs. per sq.ft. 
r = radius of droplet (ft.) 
; . U\rpa 
R = Reynolds modulus (dimensionless); R = 2—— 
“ 
t = time (seconds) 
7 = absolute temperature (°R 
U = relative velocity vector (ft. per sec. 
\ = velocity vector (ft. per sec 
x = distance (ft. 


a = wedge half angle (deg 
= “range’’ of a drop in still air ({t.) 


“a = viscosity of air (Ibs. sec. /ft.*) 
p = density (slugs per cu.ft 
w = shock wave angle (deg 


Subscripts 


a = air 

w = water 

0 = initial value 

l = before shock wave 
2 = after shock wave 


Superscripts 
Dots = differentiation with respect to time 


ANALYSIS 


The water interception characteristics of a wedge at 
supersonic velocities may be computed in a rather 
simple fashion subject to the following postulates: 

(a) The water drops are spherical. 

(b) The airflow about the wedge corresponds to two- 
dimensional frictionless flow of a perfect gas with con- 
stant heat capacities. 

(c) The drag force exerted by the air on the drop is 
that experienced by a sphere in a real fluid. 

(d) The shock wave is not detached. 

(e) Nocondensation shock occurs. 

(f) No heat transfers or mass transfers occur to or 
from the drop. 

These postulates have been used in the idealized 
treatment of icing at subsonic speeds with good agree- 
ment between measured and predicted behavior. Of 
the above list, postulate (e) most likely will cause devi- 
ations between predicted and observed performance. 

For the case of a wedge, the supersonic flow field 
(Fig. 1) is simpler than the subsonic flow field. If the 
flow ahead of the shock wave is uniform and parallel to 
the chord line of the wedge, after the shock wave it will 
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m = (4/3)r r°py (2 
a= dU dt (.3 
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again be uniform but parallel to the wedge surface. 
Because disturbances cannot be felt upstream, the flow 
over a finite wedge is initially the same as that over a 
wedge that extends an infinite distance downstream. 
Fig. | reveals the interesting fact that all the droplets 
describe identical trajectories until they either strike 
the wedge surface or flow with the air stream parallel to 
it. The drops arrive at the shock wave with the initial 
Relative to the airflow after the shock 


velocity V4. 
Therefore, asa 


wave, they have a velocity ( V, — 0). 
consequence of the fact that the flow field is uniform 
before and after the shock, the trajectory of a drop rela- 
tive to still air may be determined, and the trajectory 
of the drop relative to the wedge may be obtained by 
adding the air velocity Vs. 

The ‘“‘range’’ of a sphere projected into still air has 
been computed by Langmuir.* Because some of the 
intermediate data indicated but not presented in his 
analysis are needed here, the derivation and computa- 
tion are presented in the following. 

If a drop of radius r and density p,, be projected into 
still air of density p, and viscosity u, then, according to 
Newton's Law, there will be a deceleration because of 
the drag—.e., 

F = ma (1) 


where m, the mass of the drop and 4, the acceleration are 


given by 


The drag force Fis given by 
F = —(1 2) pa U 2 oe r(U U (4 
Cp is the experimentally determined drag coefficient for 
asphere. Therefore, 
dU 3 paCp nat . 
= _— = = \») 
dt Spy l 
The drag coetticient is a function of the Reynolds modu- 


lus for the droplet R. Replacing the velocity by the 


Reynolds modulus, one obtains the scalar equation 


MY) 
S\px r 2rpr 


Finally, then, 


x dR 3 w(t — to) J 
= — (4) 

Ro Cp(R)R? 6 Tt, 

Since Stokes’ Law (i.e., Cp = 24/R) is valid for small 


values of R, if the upper limit of integration approaches 
zero, t becomes infinite; the drop never comes to rest. 
As pointed out by Langmuir, however, the drop travels 
only a finite distance 

Using the reported experimental data for drag co- 
efficients for spheres, Eq. (7) may be integrated. Curve 
A of Fig. 2 shows the resultant variation in Reynolds 
modulus, R, with time when a drop is projected into 
still air with an initial value for R equal to 1,800. 
Curve A of Fig. 2 may be used for other initial values 
of Reynolds modulus by shifting the origin of the time 
scale—i.e., subtracting the time required for R to change 
from 1,800 to the given initial value. 

To determine the range of a drop projected into still 
air, Eq. (8) must be integrated 


U| = dx/dt = (u/2r p,)R (8 


oe *(pt/ pr?) 
(59) 1m)» 
Pu 2r ft (ulo/ p,¥2) Py” 


The integral on the right-hand side of Eq. (9) is the 
area under the curve of R plotted versus (uf/p,4r° 
within the limits ¢ and f determined from the first 
The result of the second integration is 


hence, 


integration. 
shown in Fig. 2 as Curve B. 

The results of the foregoing calculations can be used 
to calculate the water impingement on a wedge at 
supersonic speeds as follows: 

(a) The flow conditions before and after the wedge 
may be found from Reference 7. 

(b) The initial Reynolds modulus of the droplet, 
with respect to the air flow after the shock wave, is com- 
puted using the magnitude of the difference in vector 
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velocities of the air before and after the shock wave as 
the relative velocity. 

(c) With this value of RX, Curve A of Fig. 2 may be 
used to determine the dimensionless time interval for 
the droplet to come to the initial Reynolds modulus 
from R = 1,800. 

(d) Corresponding to this dimensionless time u(t — 
ty) 4r°p,, a value of dimensionless distance (x — Xo) pq 
2rp, may be determined from Curve B of Fig. 2. 


(e) The difference between this value of (v — Xo) p, 
2rp,, and its maximum value (8.63) corresponding to 
R, = 1,800 is the maximum range, A,,~:,, the drop at- 


tains with respect to the air. 

({) Superposing a line of length \,,,,. in the direc- 
tion of the initial relative velocity on the sketch of the 
wedge flow field (Fig. 5), it is seen that there is a critical 
starting position, 0,,,,, such that the drops just reach 
the wedge surface. Drops farther from the chord line 
than 3,,,,, do not reach the wedge. 

(g) The position at which a particular drop strikes 
the wedge surface may be obtained by adding, vecto- 
rially, the distance traveled by the drop, with respect 
to the air, to the distance the air moves in the same time. 


A sample calculation follows. 


PROBLEM STATEMENT 


Consider a two-dimensional wedge of half angle 10° 
traveling at a Mach Number of 1.767 at an altitude of 
15,000 ft. and a free air temperature of —25°F. The 
droplets are 20 microns in diameter (6.56 XK 107° ft.) 
and the liquid water content of the cloud is 0.2 Gm. per 
cu.m. Where do the drops strike the wedge? 

From reference 7 it is found that (Fig. 3) MWe = 1.417, 
P:/P, = 1.655, pa2/pai = 1.428, a, = 1,021 ft. per sec, 
and the shock wave is at an angle of 45°. 

From these data the following are calculated: 


T> T, = (P. P)) (Pa Paz) = 1.699 1.428 = 1.156 

V. Me |T, 1.417 , 

= = 2 = 36 = 0.862 

r, MNT, 1.767 ¥ 1-196 

V, = Mya, = 1.767 X 1,021 = 1,802 ft. per sec. 

> = 0.862 X 1,802 = 1,555 ft. per sec. 

P, = 16.SSin. Hg, pai = 0.00160 slugs per cu-ft., 
T, = 435°R. 

P» = 27.9 in. Hg, paz = 0.00228 slugs per cu.ft., 
To = 504°R. 

The relative velocity of the droplet ( V, — PV.) may be 


computed from the vector diagram of Fig. 4. 


U) = V1,5552 + 1,802? — 2 X 1,802 X 1,555 cos 10° 


U) = 379 ft. per sec. 


The initial Reynolds modulus of the droplet is, there- 


fc re, 
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M,=1.767 
P,=16.88"Hg 
IV,1=1802 ft/sec 


M,=1.417 
P,=27.9"Hg 
IV1=1555 ft/sec 











f,=O.00160 slugs/ft® —P4= 0.00228 slugs/ft® 
T= 435°R T,=504°R J 
=! , = = 
‘ ™ - 
Vv, 
V,= 1802 
U=(Vi-4) 
—~ 3471599 
Fic. 4. 
U| 2rp, 379 XK 6.56 KX 10 & 0.00228 
RR = = : : = 154 
u 0.368 K 10- 


The dimensionless time required for the drop to come 
to this value of R from an initial value of R = 1,800 
is read from Fig. 2, Curve A, as 


u(t — to) /4r2p, = 1.25 & 10 
From Curve B of Fig. 2, the dimensionless distance 
traveled in this time (from an initial value R = 1,800) is 
(X — Xo0)pa/2rpy = 6.10 


hence, the (dimensionless) maximum distance traveled 
by a drop initially at R = 154 is 


Amar. Pa/ 2’ Pw = 8.63 — 6.10 = 2.53 
(8.63 is the range of a droplet initially at R = 1,800). 
In feet, this range is 


2.53 X (6.56 K 107°) X 62.4 


2.53(27) po 
Amar, = = 


Pa 32.2 XK 0.00228 
0.141 ft. 


The time required for the droplet to traverse half this 
distance, for example, may be computed from Fig. 2 
as follows: 

(1) The initial dimensionless position at time zero 
is 6.10. 

(2) Traveling half the ‘“‘range’’ increases the dimen- 
sionless position by 2.53/2 = 1.26—1.e., the dimension- 
less position will be 6.10 + 1.26 = 7.36. 

(3) The dimensionless time (Curve B) is 2.6 X 
10~* corresponding to position 7.36. 

(4) At position 6.10, dimensionless time equals 
1.2 X 107°. 
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TABLE | 
Distance the 
Air Stream 
Distance Traveled Time Required Distance Traveled Moved Parallel 
by Droplet Relative to Traverse Time by Droplet Rela to Wedge in 
to the Air Stream This Distance Required tive to Air Stream This Time Interval 
( Dimensionless ( Dimensionless (Sec Ft Ft 
2.53 0.141 
1.90 3.3 X 10°? 7.46 X 1074 0.106 1.16 
1.26 14 xwH> 3.16 X 10 0.070 0.49 
0.63 0.58 X 107? 1.31 X 1074 0.085 0.20 
TABLE 2 








a=wedge half angle ~ 





Fic. 5 
(5) Therefore, the time interval is 1.4 & 107? (di 
mensionless). 
(6) The actual time is, therefore, 


(1.4 X 1072) (2r)2p, 
t —f= = 


5)? x 62.4 


x (6,56 X 10 
xX 0.368 X 10°° 


to 


3.16 X 
Table | was constructed using the method just de- 


scribed. 


From the known shock angle (w), wedge half angle 


(a), and direction of the relative velocity vector () 
(Fig. 5), the initial distance from the centerline (0) of a 
droplet that strikes the wedge by moving the relative 
distance \} may be computed from trigonometric rela- 


tions as 
b = N{sin w sin (a + B)]/[sin (w — a)] 


The angle (38) of the relative velocity vector may be 
computed from 


; sin @ 
sin p = _ 
Vi1+ (Vi/ V2)? — 2 (Vi/ V2) cos a 
Vs! 
Sin @ 
For the case at hand, the angles are a = 10°, w = 45°, 


B = 45.5°, and 


sin 45° sin 55.5° 
b=) ree " 
sin 35° 


1.01X 


Using these data Table 2 was constructed. 


10~4 sec. 


Average Rate of 
Ice Accretion 
Between This 


Position on Point and Nose 


Initial Wedge at if Liquid Water Average 
Distance from Which Drop Content Is 0.2 Rate of 
Centerline Strikes Gm. per Cu. M Icing 

(d, ft. (Ft. from Nose) (Lbs. Hr./Ft.? In. per Min 
0.0385 0.24 12.0 0.045 
0.071 0.56 10.2 0.036 
0.107 1.27 6.8 0.024 
0.1438 0) 0 


The position on the wedge at which the drop strikes 
differs from the distance traveled by the air stream 
parallel to the wedge surface (Column 5 of Table | 
This difference is the initial displacement of the drop 
relative to the nose of the wedge but parallel to the 
wedge The 
shock wave angle is less than a right angle. 


because the 
This initial 

displacement is equal to (6/sin w (90° — w 

15°, the initial displacement is 


surface. difference occurs 
[sin 
For the case where w = 
equal to b, the initial distance of a drop from the center 
line of the wedge. 

For this particular flight speed and free air temper 
ature, it is found (using the methods of reference 8) that 
icing will not occur because of the aerodynamic heating 
effects. However, at —40°F. and the same flight speed, 
the aerodynamic heating will! be insufficient to prevent 


icirg. 
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Dynamic Stability of a Missile in Rolling 
Flight 


RAY E. BOLZ* 
Case Institute of Technology 


ABSTRACT 


[he paper sets down the equations of motion for a symmetric 
rolling missile with respect to axes attached to the missile. The 
missile may be jet (or rocket) propelled or coasting under acceler- 
ating or decelerating conditions, respectively, wherein the vari- 
able rolling velocity is derived from intentionally or uninten- 
tionally ‘“‘canted”’ fins and/or wings. 

The equations contain a force and moment system that in- 
cludes, in addition to the usual forces and moments, those due to 
magnus effects, misaligned surfaces, canted surfaces, jet misalign- 
ment, and the linear accelerations in the plane normal to the mis- 
sile axis. 

The results present general stability criteria for a rolling missile 
which are summarized in the ‘‘ Discussion of Stability.” 


List OF SYMBOLS 


(1) English Letters 


A = moment of inertia in roll about x axis, pd®, lb. 
sec.? ft. 

| = A pd‘, ft.~! 

B = moment of inertia about y axis, lb. sec.? ft 

B = B/pd‘. ft.— 

a = moment of inertia about z axis, lb. sec.? ft. 

Cp = drag coefficient, dimensionless 

( = wing or tail chord, ft. 

d = diameter of missile, ft. 

G = linear momentum vector 

H = angular momentum vector 

: = unit vectors along y and ¢ axes. respectively 

L, 1, N. = moments about x, y, and 2 axes, respectively 

n = mass of missile, lb. sec.? ft.~! 

” = m/pd’, dimensionless 

‘of = rolling velocity of missile about x axis, rad. per 
sec. 

P; = P/V, rad. per ft. 

1, ? = angular velocities in pitch (about y axis) and in 


yaw (about ¢ axis) 


») area, sq.ft. 

7 = thrust force, lbs. 

T = T/pd?, ft./sec.? 

t = time, sec. 

u,v, W = components of linear velocity of the missile 


along x, y, Z axes, respectively, ft. per sec 
derivatives with respect to time 


V = scalar velocity of missile along flight trajectory 
ft. per sec. 
X, ¥,Z = forces acting along x, y, 2 axes, respectively 


in 
ll 


designations of coordinate axes 
g = distance along trajectory 


Il) Definition of Coefficients 


K = T/V?, 1/ft. 
Kz, Km, = dimensionless coefficients of force and moment, 
F ‘ 
respectively, due to angle of attack or yaw 
KK... Km, = Kz,/d and Km,/d, respectively, 1/ft. 
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kK; , Km, = 
q 7 


A;,, Am, = 
Key Km, - 
Keyg Km yy = 
Kp 
kgaet™ 
Km, K ms\ 
Kz, Ky | - 
| a 
Kn; = 
Rm, 
ZT, Yr/ = 
Nr, Mr\ 
Zr, Y7 = 
Nr,.Mr = 
Ry = 
R; 2 
Fr . 
Mr’ = 
? if 
EK. = 


dimensionless coefficients of force and moment, 
respectively, due to pitching or yawing angu- 
lar velocities 

dimensionless coefficients of force and moment, 
respectively, due to combined roll with pitch 
or yaw 

dimensionless coefficients of force and moment, 
respectively, due to accelerations 

dimensionless coefficient of magnus force and 
magnus moment, respectively 

dimensionless coefficient of axial force (drag) 

dimensionless coefficients of forces and moments 
in 2 and y directions, respectively, due to mis- 
aligned lifting surfaces 

division of coefficients by d, 1/ft 

dimensionless coefficient of moment due to 
canted surfaces 

Km; d, 1/ft. 

forces and moments in = and y directions due to 
jet misalignment, lbs. and Ib. ft., respectively 

Zr/pdand V7/pd, ft./sec.? 

Nr/pd' and Mr/pd*. ft./sec.* 

Ky + ikz, complex coefficient, 1/ft 

Km, + iKm, complex coefficient, 1 /ft 

complex force, Vr + 7Z7, ft./sec.? 

complex moment, Mr + iN7, ft./sec.* 

coefficient of the total force due to jet misalign 
ment, 1 /ft. 

coefficient of the total moment due to jet mis- 
alignment, 1 /ft. 


(11) Greek Letters 


a = 


o> 
ll 


angle of attack (in plane of the velocity vector 
V and the Z axis) 

angle of yaw (in plane of the velocity vector V 
and the y axis) 

angle of cant of surfaces (differential angle of 
offset ) 

angle of roll 

mass density 

complex linear velocity, v + iw 

complex angular velocity, g + ir 

n/v and ¢/v, respectively 

derivatives with respect to z, the distance along 
the trajectory 

component of angular velocity of the missile 
about the x axis due to precessional motion 


(IV) Designation of Derivatives 


Primes (n’, etc.) are derivatives with respect to z, the distance 


along the trajectory. 
Dots (i, etc.) are derivatives with respect to time. 


INTRODUCTION 


HE STUDY OF THE DYNAMIC STABILITY and of the 
motion of high-speed missiles and aircraft is the 
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subject of constant investigation and analysis. Up to 
the advent of the supersonic guided missile and the jet- 
propulsion airplane, the aircraft designer needed refer- 
ence only to the classical stability theory initiated by 
Lanchester and developed by Glauert and others. A 
complete discussion of this may be seen in reference 1. 
Here, the treatment centers about the ability to sepa- 
rate the six equations of motion into two independent 
groups of three equations each, the lateral and longi- 
tudinal groups, and the stability associated with the 
symmetric and nonsymmetric planes is easily obtained. 

The missile development engineer now is faced with 
the dynamic problems concerning rolling missiles and 
the aircraft designer with aircraft having relatively 
large rotating masses in the form of turbines and com- 
pressors. The introduction of a rolling velocity into 
the equations of motion introduces a coupling term 
that prohibits the separation of the equations into 
independent lateral and longitudinal groups, and the 
six equations of motion must be simultaneously treated. 

The problem is by no means new, since the stability 
theory of spinning shells has been a subject of investiga- 
tion by the ballistician for the past several decades and 
ballistic treatments of the subject may be found in 
references 2 and 3. It is the purpose of this paper to 
present the stability criteria for a rolling missile with 
basic symmetry of design but with a limited type of 
asymmetry in the form of misalignment of the lifting 
surfaces and of the propelling jet if employed. The 
analysis includes the effect on stability of the canted 
surfaces that induce the roll either intentionally or 
otherwise. 

The analysis presented here is in no sense a funda- 
mentally new treatment and may be regarded as a 
modification to the work by Fowler, et al.,” and Niel- 
sen and Synge.* In contrast, however, the motion here 
is referred to a coordinate system fixed in the missile 
and rotating with it such that certain asymmetries may 
be most conveniently expressed. The results are not 
restricted to constant-velocity flight or constant rolling 
motion, but only to motion wherein the rolling velocity 








~ 


NY 


Fic. 1. Axis system for analysis. 
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is ‘‘self-induced’’ by canted fins such that it becomes 
a function of the linear velocity itself. 

The results present more general stability criteria 
than previously derived and indicate some interesting 
results concerning the effect of rolling velocity, wing or 
fin cant, accelerations, missile geometry, mass distribu 
tion, thrust and drag, misaligned surfaces and jets, 
magnus moments and forces, and roll-yaw interaction 


moments. 


DISCUSSION OF THE PROBLEMS AND ASSUMPTIONS 


The purpose here is to analyze for the dynamic sta 
bility of a rolling finned and/or winged missile in free 
flight. The missiles will be restricted to essentially sym 
metrical configurations with nonsymmetry due only to 
unintentional small angle of incidence settings (in pro- 
duction assembly) of the wings and or fins or small mis 
alignment of the rocket or jet blast (if any) propelling 
the missile. In the analysis, the following assump 
tions are made: 

(1) Gravitational effects are neglected for sim 
plicity. 

(2) The angular deviations of the missile from its 
trajectory are small such that the sines and angles are 
approximately equal and the cosines may be taken as 
unity. 

(3) The linear acceleration is small enough such 
that the force and moment coefficients may be con 
sidered essentially constant over a reasonable period of 
time in the analysis. For application to flight with a 
large velocity change, a series of stability calculations 
may be made for different values of the coefficients 
selected to cover the whole velocity range in discrete 


intervals. 


EQUATIONS OF MOTION FOR A ROLLING MISSILE 


Consider a missile with cruciform tail and wing sur 
faces, any one of which may be at an angle of incidence 
with respect to the fuselage. A rectangular coordinate 
system is attached to the missile at its center of gravity, 
as shown in Fig. 1. These are the typical stability 
axes and directions employed in the usual aerodynamic 
stability treatment of nonrolling aircraft. The use of 
these axes allows the simple introduction of the effect 
of wing or fin misalignment angle and the effect of jet 
misalignment. 

The six scalar equations of motion may then be 


written from the vector forms, 


>M = dH/dt (1) 
and 
>-F = dG/dt (2) 
and the equations become 
m(u + wq — vr) = DX (3a) 
m(o + ur — wP;) = DY (3b) 
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DYNAMIC STABILITY 


m(w — uq + vP,) = >Z (3c) 
AP, + Cqr — Bar = YL (3d) 
Bg + APir — CPir = YM (3e) 
Cr + BPiq — APig = DON (3f) 


where A, B, Care moments of inertia and the products 
of inertia vanish because of symmetry of mass distribu- 
tion. In these equations, P; is actually the sum of the 
roll rate of the missile P and its own axis and the com- 
ponent of angular velocity w, about the Y axis due to 
the precessional motion resulting. 


P,} = P+ a, (4) 


However, from actual experimental tests of missiles 
with canted fins, w, is, in general, a small magnitude 
compared with that of P (of the order of 0.001P). For 
small values of rolling velocity P, this would not neces- 
sarily follow, although no observations indicate other- 
wise; but here we are in the range of rolling velocities 
where the roll rate and resultant forces and moments 
are extremely small and without important influence. 
Consequently, we neglect w, in the equations of motion 
and replace P; by P without significant effect on the 
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v 
Fic. 2. Sketch of vector 7 


In Eqs. (3) the terms wg, vr, and gr may be neglected 
as second order, since the angular deviation of the 
missile from its flight path is restricted to small angles 
and the angular velocities g and r are disturbance ve- 





stability results. locities. 


From the mass symmetry, B = C, and we may write the equations of motion now simplified and including the 
forces and moments involved. We have, as in the usual stability analysis (see reference 1, page 124, for a detailed 
discussion), neglected the influence of accelerations on the forces and moments involved with the exception of the 
influence of accelerations w and 2. 


(m/pd*)t = (T/pd*) — (Kp/d)V? cos @ = (T/pd*) — (Kp/d)V? (5a) 


(m/pd*) (8 + ur — wP) = (-H,. d)jvV + KarV + K,,Pw + K, — K,,Pw + (Ky/d)V? + (¥Vr/pd*) (5b) 


(m/pd*) (w — uq + vP) = (“Rig d)wV — K.qV — K.,Pv + K..w + K.,Pv + (K,/d)V? + (Zr/pd*) (5c) 


g 


(B/pd*)g + [(A — B)/pd*)Pr = (Kn,/d)WV — KnggV — KnpPv + Kung — Km,V0 + 


6 
KimyPv + (Kym,/d)V? + (Mr/pd*) (5d) 
(Br/pd*) — [(A — B)/pd*]Pq = (—Km_/d)vobV — KingrV — KnyPw — Kn — (Km,/d)wV + Kny Pw + 
a . a 
(K»,/d)V? + (Nr/pd*) (Se) 
Here the coefficients of forces given by the A..'s and the coefficients of moments given by the A,,,’s are dimension 
less coefficients presented mathematically in Appendix (I) and the values of which depend upon the stability deriva- 
tives and geometry of the missile used. The stability derivatives may be evaluated theoretically (see, for exam- 
ple, reference 4) or by experiment. 
In Eq. (5b), a numerical comparison of the terms (m/pd*)d and K..v shows that a ratio Kz .(pd® m) is of the 
order of 10~* or 10~4 even for extremely low density missiles. A similar comparison in Eq. (5c) shows the ratio 
K,, (pd*/B) to be of order 10~*. Accordingly, these terms are unimportant and are, henceforth, neglected. 


We define, for convenience, 


. m m 4 5 Zr Fa } T 
cscieaes ]3 — ]3 4;= ]3” i, = 13 
pd’, pd*, pe pe 


AwmA pd*, B = B/pd*, My = Mr/pd', Mr = Nr/ pd 





and use a bar notation over the coefficients to indicate division by d. Now multiply Eqs. (5c) and (5e) by 7 and 
add, respectively, to Eqs. (Sb) and (5d) to obtain 
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mu = T — KyV? (6a 


ig) — mP(w — iv) = —K, Vv + tw) + K,,V(r — 1g) + V*(Ky + 1K.) + 
(Az, = K.y)P(w — ww) + (Yr + iZr) (6b 


m(o + iw) + mur — 


B(g + th) + (4 — B) (r — iqg)P = Kn Vw — wv) — K,,Vi¢ + ir) — Kn Vv + iw) — 
(Kny — Kmy)P(v + ww) + V2(Kint + 1K m2) + (Mr + iNr) + K, _(@ — 7%) (6c) 
Now we define 
gt+uwr= 9; v+we=yn (7 


where y and ¢ are complex values of the linear and angular yaw velocity, respectively. 


Eqs. (6b) and (6c) then become 


mn — imou + imPn = —K. Vn — 1K,,Vo — (K., — Kiy)Pn + ad + Fy, (Sa 

Bé — i(A — B)Pd + NK m, = —iKy, Va — K,,Vé — (Ka, — Kay) Pa + °K, + Mr - KV (Sb 
where 

Ry= Ky +iky Pr = Pr tidy i 

AK; = Am, + t1Km Mr’ = My + iNeI 


are the complex force coefficients due to misaligned surfaces and jet offset, respectively, and the complex moment 
coefficients due to misaligned surfaces and jet offset, respectively. 
By assumption, we have specified that the angle of yaw and the yawing velocities are small disturbance values, 
and we write « = J cos ¢, where cos @¢ > 1 and sin ¢ > ¢; then 
uw = V — Vsin od = V — Vo¢d (10 


oru = Vandu => V. 


A final step is now taken 
the trajectory—and, since n, ¢, and P are actual velocities, we define » = 


namely, the transformation from the independent variable ¢ to z, the distance along 
Vn., @ = Vo., and P = VP,, where 


n. =7/V,o:. = o/V,and P, = P/V are velocities per unit distance of travel alongs. Then 
n = (dn/dz) (dz/dt) = Vn’ 
or 
a= (Vn. + 9.V’)V; & = (Vo. + 6.V’)V (11) 
where the prime refers to a derivative with respect toz. Eqs. (9) become 
m(V?n.’ + VV'n.) — muVo, + mP.y,V? = —K, Vn, — 1K,,V*¢, + V2K, + Fr - 
UK, — Kiy)V?P 7 
B( Vo,’ + VV") — i(A — B)V?P.d, + i(n.’V? + VV'n2) Km. = —iK, Vn, — = 
K n_ V7. — Kn Vn: — (Kiny — Kny)V?P.m + V*K, + in 
and Eq. (6a) becomes, using Eqs. (10) and (11), 
mVV’' = T — KpV?; ~V'/V = (1/m) [((T/V?) — Kp] (13 
Divide Eqs. (12) by |? and insert Eqs. (13) and (10), 
nin.’ + imP.»:. — imo, + [(T/V?) — Kp]n. = ~—K. 9: — 1K,,¢, — (K,, — Ky) Pa, + K, + (Fr/V? (14a) 
Bo.’ + 1K» m2’ + 1 Km,/m) (T/V?) — Koln. + (B/m) V?) — Kylo. = i(A — B)P.¢, — 


(T 
IR m0: _ K ing: nat Kn n: re (Bins ‘ite se wl te + K; a (M'r 'V?) (14b) 


It has been found, see reference 5, both experimentally and theoretically, that for all of the many geometrical 
missile configurations investigated the ratio of equilibrium rolling velocity P induced by aerodynamic forces to 
linear velocity V, P/V = P.,, is very nearly a constant over a Mach Number range of 1.3 to 4.0 (see, for example, 
Fig. 14, reference 6). Then, for conditions where the linear acceleration is not large compared with the resulting 
angular acceleration in roll, P, is essentially a constant for supersonic flight. 
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Furthermore, if the missile is powered by a ram-jet, we may write the thrust approximately as 7 > A |” and 
/ => KV* toa good approximation | Appendix (2) | if the changes in velocity are not too great. Then, since 


Fr = iM} = K,V?; M,’ = Tip = K.V? 


where / = moment arm, ft.; AK, = 8K, 1/ft.; K. = BK/, dimensionless: and 8 is the complex angle of misalign- 
ment, Eq. (14) becomes completely independent of the flight velocity (for supersonic velocities), except for the depend- 
ence of the force and moment coefficients on the velocity. 

If the missile is rocket-powered, the thrust is essentially independent of the velocity, and we find the equations 
dependent to some extent upon the velocity. Here, then, we restrict ourselves to small enough changes in velocity 
such that [(7'/ 1?) — Kp] is negligible and K7/ ?— constant and W7/V?— constant. However, for a rocket, the 
mass of the missile usually varies considerably with time or distance, and the solution presented here does not 


properly treat thiscase. We shall solve the former case—namely, 


mn.’ + (K — Kp)n. + imP.n. — imo: = —K, 2. — iK.9@: — i(K., — iK.,)Pm. + Ky, + Ky | 
Bo.’ + ine’Km, + (B/m) (K — Kyo: — (A — BYP: + i(Kn,/m) (K — Ry)n = (15) 
—iK »_% sn K ng: = (K mp — K, a) P he - K,, 2: T RK, a K,\ 
SOLUTION OF EQUATIONS OF MOTION If the denominator becomes extremely small under 


5: 3 , ae ; . ; any circumstances, instability results. We investigate 
Eqs. (15) are two linear differential equations with : os oe Cs 

ape ae < under what conditions, if any, this may occur. 
constant coefficients, which we rewrite for simplicity as : 


Since ac — bd = A + 1B, 


n-’ + an. + bod: = 5 (16a) 
Pe ee ee ee (16b) n- = (cs — bt) (A — 1B)/(A? + B (21 
sales But if A° + lene 0, then A- Oand B— 0. 
From an inspection of the coefficients a, 5, c, d there 
a = (1/m)((K — Kp) + RK. + (K,, — K:y)P:|+ appears no feasibility of A and B approaching zero 
iP, simultaneously. 
b = i[(K.,/m) — 1] 
c = (1m) (K — Kp) + (K,,,/B) — 1(A/B)P; + STABILITY CRITERIA 
; i — ; Dynamic stability of a missile implies that, if a small 
d = Uh a B) + [(Kny, — & mu)/ B| P, + disturbance occurs during flight and causes a displace- 
(An,/B) + UKm,/mB) (K — Ko) ment of the orientation of the missile axes from their 
e = 1(Kn,/B) equilibrium position, the effects of the disturbance will 
s = K, + Kk, damp out and the missile will tend to return toward 
t = K+ Kk, equilibrium. Mathematically, this means that the 
: real part of Eq. (19) must always be negative. We in- 
’ , spect Eq. (19), and two criteria for stability then arise. 
nS (I) The real part of (a + c — eb) must be positive. 
Solving Eq. (16a) for ¢,, differentiating, and sub- This means then that 
stituting these into Eqs. (16) yields i, (R ~ RK) R., Kn ‘ Kn.) “en 
n-" + n-'(a + ¢ — eb) + ne(ac — bd) = cs — bt (17) ‘~ hn + hn + B B { — 
Let and since RK. | Kp, K, and Km, are always positive 
n. = Re* (18) coefficients, the relationship indicates stability for al- 


most all cases. The only exception here may be con- 
ceived under the most unlikely condition of zero thrust, 
4 = [-(a + c — eb)/2] + K = 0, and extremely high drag coupled with an ex- 
(1/2)V (a — c + eb)? + 4b(d — ae) (19) tremely small lift coefficient (extremely small surfaces) 

and small damping moment. 


(II) The second criterion for stability is that 


The complementary part of Eq. (17) becomes 


THE PARTICULAR SOLUTION 


jreal (a + c — eb)| > _ 
The particular solution of Eq. (17) is real V (a — c + eb)? + 4b(d — ae) | = | real 


ne = (cs — bt)/(ac — bd) (20) VP+ 10 
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Then the stability criterion becomes 


real (a + c — eb) ?*P + 
(P?/4) > (P? + Q?)/4 


real (a +c — eb)\4 — 
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Since R(a + c — eb) is a number much smaller than 


1, we neglect the fourth-power term and write 


Inserting the values for P, Q, a, c, etc., we get a stability equation, 


A (; : 
4) 1 (K., — Kiy) A Rn. B \2m 
i to 
B BP. 


2 m 


— P > Q?/4 [R (a + c — eb)/P (23 
K ng — Kms K mp King, Kins | 
2B | 2B BP. B a” @ 
= at 7 : : (24 
y(K- Ro), Rig | Km Kms 
7 m ™ m af B B 


In this relation the magnitude of the A’s for a finned missile is of the order of 10~* or smaller, and, as such, the 


products of these terms are neglected. 


stability. 


In some instances the ratio of the moments of inertia 4 /B may be considerably larger than 10 


The equation therefore is not valid for missiles with near-neutral static 


+ and 


for this reason the product terms involving this ratio are retained. 


For the case where the term 4/B is of the same order of magnitude as the A's, Eq. (24) may be further simpli 


fied to 
Kn, 4 K mp 4: Ma ie u P 
—K, . BP, B B B bs 
» > “ (Zo 
BP i (AK a Kp) A. Kn, 
mi m B B | 


For cases where the term (A /B) is of the same order 
of magnitude as the A’s but Ru. BP.” is of the order 
of the square of the other A’s (near-neutral static sta- 
bility), neither equation is valid, and Eq. (23) must be 
used and analyzed in its entirety. 

(III) The third criterion for stability is that Eq. 
(21) is never satisfied. 
efficients, it appears impossible in a practical sense for 


From an inspection of the co- 
this expression to vanish. 


DISCUSSION OF STABILITY 


From Eq. (2+), the following results may be stated as 
applied to a missile in (a) accelerated or decelerated 
supersonic flight where the rolling velocity is obtained 
aerodynamically and varies with the velocity or (b) 
a supersonic or subsonic missile in essentially constant- 
velocity flight. 

(A) Fora statically stable missile with zero or small 
rolling velocity, dynamic stability is assured, since 
Re. is a negative coefficient and | is zero or ex- 
tremely small corresponding to the rolling velocity. 

(B) For a statically unstable missile with zero or 
small rolling velocity, dynamic instability results. 

(C) For rolling rates in the range of '/, to 2 or so 
rad. per ft. travel, dynamic stability is not at all assured 
for a statically stable missile unless relatively large 
numerical values of Ka. are employed. The left side 
of the inequality has an order of magnitude of the cube 
of the coefficients (K’s) compared with the square on 


the right, although the coefficients retained need not 


themselves be of the same order of magnitude. 


(D) For large roll rates, dynamic instability can be 
induced for a statically stable missile, whereas for a 
statically unstable missile, dynamic stability can be ap 
proached and in some cases obtained by large roll rates. 


(E) Large values of thrust have a stabilizing influ- 
ence, while large values of drag have a destabilizing in- 
fluence. 


An increase in missile size (diameter) for given 


(F) 
values of the linear and roll velocities decreases the terms 
| | ) aS | 


a 6 
° a 
d* 


BP d? BP, 


by the square of the diameter with no effect on the re- 
maining terms. Increasing the size of the missile rap- 
idly decreases the stability because of the major im- 
portance of the | term in Eq. (24). It should be 
noted, however, that for a missile with a rolling ve- 
1/d and the 


~ 


locity induced by cant on the fins, P; 
stability is independent of missile size. 


(G) For a statically stable rolling missile with 
small A/B values (as is usual), the coefficients arising 
from the (1) canted surfaces, (2) rolling motion, and (3 
linear cross accelerations contribute toward instability 
tendencies. It is interesting to note that the magnus 
moment opposes this tendency. 
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TABLE I 
FORCE EQUATIONS USED IN STABILITY DISCUSSION 
| 3c P 2 C ¥ 
Forces due to angles z ats =— K, w V( ed.) K, = hac 2 Sy 
| | | of attack and yaw, 3X ‘ - e. d% 2d 
| respectively Y= —K, vV( ed ) Always positive 
| a 
| 
‘= * 
dC P - dC 
Forces due to pitching | Z: cat oS LV _¢ omK, qV(ed>) | K o all £S_ 
and | 9 (ae 2V q | 2q qc\ 443 
> | yawing angular | OV 3 DV 
| ~ | velocities, respectively 
Ys Kg TV(ed ) | Always positive 
| | C | Ye 
| [Forces due to rolling | Z = ; a oa: fs -—K, Pv(od) KS mo ab 
| 3 | Combined with pitching | ) ¥) P | pf (5) 4d 
and yawing angular 3 | 
velocities , respectively | Y= ." Pw(ed ) | Aways positive | 
| ; = _—s 7 
ac 2 ac 
| Forces due to | Z= # es Ky w(ed ) K,0* — + Sc. | 
4 | accelerated motion 4 (Se K 9(Se\ 4a | 
| | 2v) s 2V 
| | Y=Kze V(ed ) | Always positive 
iia on 
| . | s P Note: Derived from Kutta-| Kz, is a dimensionless 
5 | Magnus pitching and side bith he Joukowski lift equation | positive or negative 
4 
| ov forces | Y2- Key Pw( ed) Z=vi where ['~ P | experimental coefficient 
= 2 + 
| | ¢: Bhiwg kes a Sigh | «eel 
Forces due to mis - om m3 z e i, mae 
| 6 | aligned lifting surfaces XC 
| | Y = Ky V (Pd) Ik = bm 
; | ev" 2 | S 
as 
| 7 | Axial force or drag X=-C) 5 S*= —K, V (ed) | Ky = Cp - 
| Always positive 
* The stability derivatives may be evaluated theoretically; see Ref.4 for rectangular fins 
| (for example) and Ref.6 for fuselage. 
— 








(H) For a supersonic missile (coasting or ram-jet- 
powered), the dynamic stability is independent of the 
linear flight velocity magnitude; whereas, in subsonic 
flight this independence does not exist. The moment 
due to accelerated motion in the (y-z) plane results in 
decreased stability, although this could be important 
only for large roll rates. 

Eq. (24) derived from the particular solution results 
in the following conclusions: 

(1) Misalignment of the fins or wings has no effect 
on stability. 

(2) Jet offset for a propelled missile also has no 
effect on stability. 

It is known that asymmetry in mass distribution or 
in lifting surface area can result in instability for a roll- 


ing missile, as shown in the restricted method of refer- 


ence 7. However, such asymmetry cannot be treated 
by the methods of this paper. 


APPENDIX (1) 


The forces and moments acting on a rolling finned 
and/or winged missile are given in Table I and for the 
most part need no further discussion. The stability 
derivatives may be evaluated theoretically (see refer- 
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Fic. 3. Sketch for analysis of canted fins 
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TABLE IL 
MOMENT EQUATIONS USED IN STABILITY DISCUSSION 
| T e | 
| Cn 3 gc 
Moments due to angles | M= a ev Sc: Km wV(ed ) Km : ng —_- | 
| | of attack and yaw, 
respectively o— K.., vV(Od> ) Always positive for 
Moe statically wrstable 
ac 2 aC 2 : 
Moments due to pitch- M = er Ae ey SC=-K,, qV(od*) Km “ae £ = 
9 (ing and yawing angular Oy) 2V 2 q q (55) 4d 
velocities, respectively 
Nom TV (od4 Always positive 
m 
3c 2° 
Moments due fo rolling M=— Kinp Pv(od ) Kon Oe a a 
3 | combined with pitching als : 4d 
and yawing angular N=— K,, Pw(ed) 
| velocities , respectively Always positive 
r OO “" ° 
° 4 dom Sc* 
Moments due to M=KaA, w( d ) K = 
e Mo (Xe\ 
4 | accelerated motion 5 P x SV) 4d 
eK. v(ed ) ia 
ways positive 
on = — Ren aE 
4 
| 5 Magnus pitching and M=K,, Pv(@d ) Kmy 'S @ dimensionless, 
yawing moments mM 4 positive or negative 
N= Km, Pw (ed ) experimental coefficient | 
| ac S 
2 3 m c 
Moments due to mis - M=Km, V (ed) Km, — aa" 
6 aligned lifting surfaces oe - | i 2p 
_ 302 ev 3 dCz Sa 
Moments due to canted M = ya = —Km Vv(ed ) Km, ory ‘= 
7 | surfaces 8 d 
N =— Km Vwed | Positive for positive 
8 roll induced by cant angle | 
*The stability derivatives may be evaluated theoretically, see Ref. 4 for rectangular fins | 








ence 4, for example), and the fuselage effects may be 
included from data such as appear in reference 6. 

The moments due to canted fins may be briefly dis- 
cussed. Consider two canted surfaces as shown in Fig. 
3. The moment resulting from the axial components 


of the lifting forces is a couple—namely, 


ac. V2 
M = — —([8+8+ (¢@ — 88 — aS 
Oa Z 
Oc. aS 
M =| -— ‘)5 Vpd? = —K,, Vupd? 
| (5) go | ” eo 


as appears in Table IT. 
The forces and moments due to jet offset are merely 
the thrust multiplied by the constant misalignment 


angle or merely constants expressed as Zr and J as 
forces and Ny and \/7 as moments. 

APPENDIX (II)—DERIVATION OF THE SIMPLE THRUST 
EQUATION FOR A RAM-JET 


Symbols: 
t = static temperature 
7 = total temperature 
pb = static pressure 
P = total pressure 
C, = specific heat at constant pressure 
M = Mach Number 
VY = velocity 
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DYNAMIC STABILITY 


S = area 
m = mass flow, slugs per sec. 
Derivation: 


From Fig. 4 we write 


Facet = m(Vs — Vo) (1) 


where 
m= paS4 | 4 (2 ) 


Then 


0 , \ to — , 
(=) = T = (1 + : = Me) (4) 


and it can be shown from simple thermodynamics, 
assuming no loss in total pressure during combustion, 
that 


Vo = ¥ 2C,T>[1 a (po Py)" 1 Nh 


- — ar (6) 
Vs = V2C,T3[1 — (po/Po)%~ P7119 

Use of Eqs. (2) to (6) in Eq. (1) yields 
Fone => y¥Mo*Sspo (1 = Vv T> T3) (7) 


or for relatively small changes in flight speed only where 
the ratio 7) 


‘T; does not change greatly we may write 
Eq. (7) in the approximate form, 
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Canted fins or wings 
Relative 
Wind 
Fic. 4. Schematic of ram-jet power plant 


Fret = K V0? 


where A = constant. 
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Two-Dimensional Los 


ses in Turbine Blades’ 


CHARLES A. MacGREGOR 


Propulsion Research Corporation 


SUMMARY 


Tests have been made on three typical turbine blades in a two- 
dimensional cascade tunnel. The were run at varying 
angles of attack, stagger angle, solidity, and downstream survey 
stations at a Mach Number of about 0.4 and a Reynolds Number 
of 6.5 X 105. The static pressure distribution about the blade 
was measured by 27 static pressure taps. The wake behind the 
blade row was surveyed and the loss determined 

Calculations were made of the boundary-layer loss at the vari- 
ous angles of attack, stagger angle, and solidity using the meas- 
The calculated loss was found to 


tests 


ured pressure distribution. 
agree with the measured loss even though the pressure distribu- 
tion varied greatly, It was found that the loss on the suction 
side of the blade was about six times as large as the loss on the 
pressure side. This was the result of the shape of the pressure 
distribution about the blade. 

Discharge measurements made behind the blade row showed 
that the wake spread out as it moved downstream. The loss did 
not increase greatly in the first half-chord length behind the 


blade row, even though the shape parameter did 


INTRODUCTION 


HE RAPID DEVELOPMENT of the turbojet and the 
yew engines, with their great dependency 
upon component efficiencies, has made it necessary to 
design turbines from an aerodynamic point of view. In 
order to provide an effective method of doing this, it 
is necessary to be able to calculate the losses in a turbine. 
The British'~* have made attempts to isolate the losses 
The losses are divided into two gen- 
i.e., profile 


in turbomachines. 


two-dimensional losses 


eral 

losses 
secondary flow losses, tip clearance loss, etc. 
mated that the profile and the three-dimensional losses 
are of about equal magnitude at the design point. This 
paper shows that it is possible to calculate part of the 
losses, the two-dimensional losses, with reasonable ac- 


groups: 


and three-dimensional losses—t.e., annulus loss, 


It is esti- 


curacy. 

A boundary layer grows next to the blade surface 
(the profile loss) because of viscosity effects and regions 
of adverse pressure gradients as the gas flows through 
the blade row. In order to find this loss, it is necessary 
to know the pressure distribution about the blade which 
may be either measured or calculated. Methods for 
calculating the velocity distribution about a blade for 
any solidity or camber have been developed from in- 

Received November 27, 1951. 

* The research on which this paper is based was done during 
1947-1948 while the author was an Aerodynamicst at The Turbo- 


dyne Corporation. The author expresses his thanks to Mrs. C. 


D. Oswalt for performing the computations and reduction of 


data. 
+ Aerodynamicist. 


compressible potential-flow theory. * However, these 
methods do not consider compressibility and are mathe- 
matically complex. An alternate approach to calcu- 
lating the velocity distribution has been the stream- 
filament theories, where the flow equations are inte- 
grated graphically’ or mathematically." * The latter 
method has been compared with experimental values 
and has been found to be in close agreement. Thus, a 
concise, simple method of determining the velocity dis- 
tribution about turbine blades is available. 

The cascade tunnel has been used for a long time in 
the tests of turbine blade profiles. Knowing the inlet 
conditions and surveying the discharge, it is possible to 
find the loss in the fluid as it flows through the blade 
passage. Three blades of greatly different character- 
istics were tested in a two-dimensional cascade tunnel 
in which it was possible to independently vary the 
angle of attack, and solidity. The 
static pressure was measured around the blade by static 
Calculations were then made, using 
distribution about the 
These calculations 


stagger angle, 


pressure orifices. 
the measured static 
blade, for the boundary-layer loss. 
are compared with the measured values obtained from 


pressure 


a discharge wake survey. 










SYMBOLS 
c = chord, see Fig. 1 
H = ratio of displacement thickness to momentum thickness 
K = ratio of momentum thickness to energy thickness 
k = constant in equation for skin friction 
n = constant in equation for skin friction 
POSITIVE NEGATIVE 
ANGLE OF ANGLE OF Se 
ATTACK ATTACK 


¥ |. DISCHARGE ANGLE 


Fic. 1 Diagram of flow through cascade blades 
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SLIDING STRIP 


BLADE GROUP _ 4, AIR 
\. (TO BE TESTED EXHAUST 
2 
/ 
AIR IN — Med 
=. TURNTABLE - . | > TUBE 
TRAVERSE 
7 V4 7 
Ss / \ 
. \© BOUNDARY - LAYER 
~ REMOVAL SLOT 


ADJUSTABLE ENTRANCE 
. SIDE PLATES 


Fic. 2. Cross section of two-dimensional cascade tunnel 





aN 


BLADE NO. | 





a 
S 
Te 
= 
i BLADE NO. 2 
——— 
BLADE NO. 3 
Fic. 3. Blade tested in cascade tunnel. 
P = pressure 
7 = dynamic pressure 
R = Reynolds Number based on chord and free-stream ve- 
locity 
Ra = Reynolds Number based on boundary-layer thickness 
s = distance along blade surface 
S) = position on surface at start of integration 
t = blade pitch, see Fig. 1 
u = local velocity inside boundary layer 
U = local velocity outside boundary layer 
() = free-stream velocity, taken as integrated mean discharge 
velocity 
x = distance normal to discharge flow 
a = angle of attack, see Fig. 1 
38 = stagger angle, see Fig. 1 
y = discharge angle, see Fig. 1 


6* = displacement thickness, JS {1 — (u/U )\dx 
6 = momentum thickness, S (u Uo) [1 — (4/ Ud) |dx 
6) = momentum thickness at transition point 


p = density 


o@ = solidity, chord/pitch 

tT) = skin friction per unit area 

¢ = energy thickness, S u/U )? [1 — (u/U5) |dx 

w = blade loss factor 
Subscripts 

a = axial 

d = downstream 

pr = pressure side 

s = static condition 

su“ = suction side 

t = total condition 

“u = upstream 


CASCADE BLADE TESTS 


Tests have been made on three turbine blades in a 
two-dimensional cascade tunnel shown schematically in 
Fig. 2. The tests were made at varying angles of at- 
tack, stagger angle, solidity, and four survey stations 
downstream of the blades. The tunnel discharged into 
the atmosphere. The blades had a chord of 3 in. and 
an aspect ratio of 2. The blade sections that were 
tested are shown in Fig. 3. These three blades are of 
greatly different shapes and are representative of the 
blade shapes found in a turbine. The pertinent design 
data for these blades are given in Appendix (A). The 
blades were tested in groups of six or seven blades in 
order to eliminate end effects and simulate the flow 
through an infinite grid. The center blade of each 
group was equipped with 27 static pressure taps, which 
were located at the mid-span position. These taps were 
used to obtain the velocity distribution about the blade. 
The tunnel was equipped with boundary-layer removal 
slots, which were used to remove the boundary-layer 
air at the inlet of the blade row. A survey made at the 
discharge of the blade row showed that the flow was 
uniform and two-dimensional. 

A traverse was made with a claw-type yawmeter, 
used by the null method, at the mid-span height of the 
discharge of the center passage. In addition, measure- 
ments were made of the temperature and total and static 
pressure upstream of the blade row. The blade loss 
factor, which is defined as follows: 


energy lost 


blade loss factor = w = > (1) 
energy available 


which, for incompressible flow, is 


| f t Pu = Pr 
oss factor = ow = 
P Pe 


tu 


to 


was found by plotting the survey points and integrating 
over the blade passage, assuming a uniform discharge 
static pressure equal to atmospheric pressure. 

The tests were all run at a Mach Number of about 
0.4 and a Reynolds Number of about 6.5 & 10° based 
upon the blade chord length. 











106 JOURNAL OF THE AERONAUTICAL SCIENCES 


METHOD OF CALCULATION 


The method used to calculate the boundary-layer 
momentum thickness about the blade is that of refer- 
ence 10, where, in order to integrate the von Karman 
boundary-layer momentum equation, 


dé HT + 2 6dq T) a 
= = 
ds 2 g ds 2g 


This equation was evaluated graphically from the 
measured pressure distribution in two parts: one for 
the laminar region and one for the turbulent region. 
The point of laminar separation was calculated by the 
method of reference 11 and was assumed to be the 
transition point. The distance between the velocity 
peak and the calculated laminar separation point was 
only about 5 per cent of the chord length. The values 
of n and k were determined, as recommended in refer- 
ence 10, from the initial and final values of R, for the 
turbulent region. A value for H of 1.4 was used. 

The boundary-layer momentum thickness was calcu- 
lated from Eq. (5); however, this momentum loss can 
be expressed in terms of a total pressure defect. The 
total pressure for incompressible flow at some point in 


the boundary layer is 
P, = £, + (1/2)o8 (6a) 
and the total pressure in the free stream is 
FP, = P, aad (1 2)pUo? (6b) 
For a uniform discharge static pressure, the loss in 
total pressure at some point in the boundary layer at 
the trailing edge is then 
(AP,) = (1/2)p(U0? — u?) (7) 
Integrating this loss across one blade passage at the 
trailing edge and weighting for flow, 


me ny p 
I (AP,)u, dt = J a (Uo? — u*)uq dt (Sa) 
0 0 “= 


and letting 


- Ne... 
(AP)u, dt (Uo? — u*)u, dt 
> 0 0 2 


Af, = _ = - 
Hal Ual 





(Sb) 


where the bar values represent integrated mean value. 
It follows, then, from the definition of the loss factor 
that 


AP, 1 f ( u? \ 
o= ——_ = -—.); at 
(1 2)pl 07 tf cos Y 0 l o- l 0 as 4) 
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the skin-friction coefficient is expressed as a function of 
the boundary-layer Reynolds Number, 


and the boundary-layer shape; therefore, H is assu:ned 
constant over the chord. The resulting integrated 


boundary-layer momentum equation is 


0-0 ee LO "OO Ty 


Expressing in terms of the distance normal to the 
discharge flow instead of parallel to the trailing edges, 


"s/Cosy 4 u 
/ (: ee ) dx 
0 Us Uo 


wo =* 4 


t cos ¥ 


| *x/COS 4 u- ( Uu ) 
ee dx| 
kao, 
l + i 
x cos 4 Uu Uu 
| (. = - dx | 
0 Uo Uo 


From the definition of 0, ¢, and A and for the values 
at the trailing edge equal to the sum of the values on the 


then 


(9b 


pressure and suction sides, then 


(*) ( o \(“ + ') (10 
ee Ne pr+su \COS ¥ K . 


This then determines the total pressure loss in terms of 
the sum of the momentum thickness to chord, solidity, 
discharge angle, and boundary-layer shape parameter 
K at the trailing edge. 


RESULTS 


A comparison of measured and calculated momentum 
thickness for the three blades at various blade settings 
is shown in Fig. 4. The measured values have been 
corrected for the traverse not being normal to the dis- 
charge flow. The comparison is made using the mo- 
mentum thickness, since it could be obtained directly 
from the discharge measurements and calculated di- 
rectly using the measured pressure distribution. In 
general, the agreement is good between the calculated 
and measured momentum thickness and is felt to be 
within the experimental error. The accuracy of the 
measured momentum thickness is limited because of the 
number ef points that were taken in the traverse in 
the wake. There were usually from six to ten points 
taken within the wake region. This close agreement 
also indicates that there was no pressure drag loss 
arising at the trailing edge of the blade or any appre- 
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ciable loss when the pressure and suction side boundary 
layers join. 

Fig. 5 shows a comparison of the measured and cal- 
culated blade loss factor. The calculated blade loss 
factor was determined from Eq. (10), which involves 
the thickness, the discharge 


calculated momentum 


angle, and the ratio of the momentum thickness to the 


IN TURBINE BLADES $07 


energy thickness (A). The value of A used in these 
calculations was 1.280, which was the average value 
the 0.1 
downstream of the blades. The agreement between the 
measured and calculated loss factors is not so good as 


determined from measurements made chord 


the comparison of the measured and calculated momen- 
tum thickness. This is probably the result of using 
this constant value for A as representative of the aver- 
age of the pressure and suction surfaces at the trailing 
edges. 

The blade setting shown in Figs. 4 and 5 does not in- 
clude any where the flow has separated from the sur- 
face, since the calculations would not be applicable. 
The velocity distribution about the blade varied greatly 
with stagger angle, solidity, and angle of attack. At 
zero angle of attack, the velocity on the pressure side 
reached a peak at about 10 per cent of the chord from 
the leading edge and then accelerated to the free-stream 
velocity. Only at large negative angles of attack did 
the velocity at the peak reach the free-stream velocity. 
On the suction side, the flow accelerated to a velocity 
peak at about 50 per cent of the chord and then diffused 
to the free-stream velocity at the trailing edge. In 
general, the velocity was greater than the free-stream 
velocity over most of the suction surface. The meas- 
ured velocity distribution at three angles of attack is 
shown in Fig. 6. Because of the difference in velocity 
distribution over the two surfaces, the calculated loss of 
the suction surface was about six times as large as the 
loss of the pressure surface. The effect of an increase 
in stagger angle and solidity was to decrease the velocity 
peak on the suction surface and, hence, to decrease the 
loss. 


The measured momentum 
from the integral of the measurements made in the wake 
behind the blade. Because of turbulent mixing, as 
the wake moves downstream, the wake becomes wider 
and the velocity distribution in the wake changes. 
The variation of /7 with downstream distance is shown 
in Fig. 7 which changes from its trailing-edge value to 
The value shown in Fig. 7 is 


thickness was obtained 


unity far downstream. 
for the sum of the pressure and suction sides, and it 
changes rapidly in the first half-chord length behind 
The variation of momentum thickness with 
The loss 


the blades. 
downstream distance is also shown in Fig. 7. 
increases slightly as the wake spreads out and is mixed 


up. 
CONCLUSIONS 


An analysis of turbine blade cascade test data indi- 
cates the following conclusions: 

(1) The two-dimensional losses of the flow through 
a turbine blade cascade have been found to agree with 
the calculated boundary-layer loss. In order to evalu- 
ate the loss, it is only necessary to know the pressure 
distribution about the blade. 
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: TABLE | 
1.5 Cascade Blade Design Data 
Blade No. 1 Blade No.2 Blade No.3 
° Typical section of: Rotor trip Rotor root Nozzle 
= 1.07 Entrance angle, deg. 22.5 15.8 ) 
5 Turning angle, deg 60 ; 85 52 
Surface length, in. 6.85 7.03 7.48 
o Chord, in. 3.09 3.00 3.42 
5r- Trailing-edge thickness, 
< in. 0.045 0.063 0.064 
3 4 
> 0 SUCTION SIDE 
= PRESSURE SIDE (3) The wake spreads out as it moves downstream 
S from the blade. However, the additional loss is small 
> 5 in the first half-chord behind the blade. 
1.0 (a) O° ANGLE OF ATTACK APPENDIX (A) 
The three blade sections that were tested correspond 
1.5 | to typical turbine blade sections. Blades | and 2 cor 
respond to rotor blade section, and Blade 3 to a nozzle 
° blade section. The blades were of die cast aluminum. 
2 1.0; To obtain a good surface finish, the blades were coated 
> with a primer and then hand polished. Some of the 
é blade design information is given in Table 1. 
._ 
a 
« 
REFERENCES 
> Py SUCTION SIDE 
3 PRESSURE SIDE ! Howell, A. R., Fluid Dynamics of Axial Compressors, Pro 
° ceedings Inst. Mech. Engr., Vol. 153, p. 441, 1945 
a 2 Howell, A. R., The Present Basis of Axial Flow Compressor 
* Sr Design, British Aeronautical Research Council, R. and M. No 
2095, 1942. 
3’ Ainley, D. G., Performance of Axial-Flow Turbines, Pro 
1.0 (b) +5° ANGLE OF ATTACK ceedings Inst. Mech. Engr., Vol. 158, p. 230, 1948 
(Continued on page 432 
1.5r- 
| 
3 
2 10F- | 
=" “7 
- || 
° — a 
— 5 +— 
- = 
= 7 
« 
> SUCTION SIDE 
5 PRESSURE SIDE 
° 
be 
w 
> or 
1.0 (c)+10° ANGLE OF ATTACK 





Fic. 6. Velocity distribution about Blade 1 at various angles of 
attack. £8, 13.8°; and a, 1.25 


(2) The loss on the suction side of the blade was 
about six times as large as the loss on the pressure side. 
This is because of the shape of the velocity distribution 
around the blade—i.e., the diffusion of the fluid down- 


stream of the ‘‘throat’’ on the suction side. 
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CHORD LENGTHS DOWNSTREAM 


Fic. 7. Variation of H and 6/C with downstream distance 
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Ripple-Type Buckling of Sandwich Columns’ 


A. CEMAL ERINGEN? 
Illinois Institute of Technology 


SUMMARY 


The problem of local instability of thick sandwich columns 
having orthotropic cores is treated under two different types of 
loading: (a) axial compressive loads acting at the edges of face 
plates and (b) bending moments acting at the ends of the sand 
wich column. End points of the face plates are assumed to be 
pinned. The present method is, however, applicable to general 
end conditions and for arbitrary axial and transverse loading 
It is shown that the deflected shape is neither symmetric nor 
antisymmetric, but asymmetric. Numerical calculations and 
curves are made for various types of cores employed in industry. 
The results are found to be in good agreement with the experi- 


mental data 


NOTATION 


FE, Ex, Ef = Young’s moduli of faces 1 and 2, or common 
modulus if they are identical, lbs. per sq.in 

Ex, Ey = Young’s moduli of core in the x and y direction, 
respectively, lbs. per sq.in 

G = shear modulus of core, lbs. per sq.in 

hh, Ie = moments of inertias of faces 1 and 2 with respect 
to their median planes, respectively, in.* 

i = length in x direction, in. 

m = number of half wave lengths in x direction 

0 = subscript to indicate boundary 

5 = subscript to indicate Fourier sine transform 

t = thickness in y direction, in 

ti, to, ty = face thicknesses for faces 1 and 2, or common 
thickness if the faces are identical, in. 

u, 7 = components of the displacement vector, in. 

a = Cartesian coordinates 

m Vxy } = Poisson’s ratios 


( contraction in y direction) 
Vay ; : = - 
? extension in x direction 


v Vyx \ 


Scr = wrinkling stress, lbs. per sq.in. 
Ox, O75 = normal stresses, Ibs. per sq.in. 
T = shear stress, lbs. per sq.in. 


Subscripts after comma indicate differentiation (¢7,, = Oox/Oy) 


INTRODUCTION 


yp THE PAST 10 YEARS, many investigators 
treated the problem of buckling and bending of 
sandwich type of beams, plates, and shells.'~'* In 
all of these works various simplifying assumptions are 


Received October 31, 1951. 

* The paper is condensed from the section on application of 
the author’s thesis submitted in partial fulfillment of the require- 
ments for the degree of Doctor of Philosophy in Applied Me- 
chanics at the Polytechnic Institute of Brooklyn, June, 1948. 
The work contained in this thesis was part of an investigation 
carried out at the Polytechnic Institute of Brooklyn under con- 
tract with the Office of Naval Research of the U.S. Navy (Con- 
tract No. N6onr-263-Task Order II). The author is indebted 
to Dr. N. J. Hoff, his thesis adviser, for the criticism and advice 
received. 
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made so as to make the problems tractable. These as- 
sumptions are mostly concerned with the stress distribu- 
tion in the core. Otherwise, a boundary value problem 
of elasticity for a two- or three-dimensional medium 
under general boundary conditions must be solved. 
Assumptions made in literature either eliminated or 
oversimplified the transverse stress components in the 
core, which in turn either eliminated or oversimplified 
the local instability. Thus, little attention is devoted 
to this ripple-type instability. As a result, different 
authorities seem to hold opposing views regarding the 
phenomenon of ripple-type instability of wrinkling. 
In particular, European investigators found that 
ripple-type instability would always occur in an asym- 
metric pattern, while other authors obtained wrin- 
kling in a symmetric pattern.' The buckling load 
is known to be largely dependent upon the deflected 
shape pattern in sandwich structures. Therefore, 
present work was undertaken to investigate the wrin- 
kling for completely arbitrary deflected shape patterns 
of faces. Hence, symmetric, antisymmetric, or skew 
patterns would be special cases of the problem treated. 
This analysis required the solution of a two-dimensional 
mixed boundary value problem of elasticity for a rec- 
tangular orthothropic medium to obtain the correct 
stress and deformation distributions in the core in terms 
of arbitrary distribution of the mixed boundary 
values. 

The present problem became tractable as soon as the 
foregoing boundary value problem was solved? (refer- 
Thus, the 


> 


ence 13 or Navy report footnote, page | 
problem of ripple-type buckling of a thick sandwich 
beam having orthotropic core is solved under axial 
compressive loads and under end moments. In this 
way it was possible to prove that the wrinkling should 
occur in an asymmetric pattern rather than the pat- 


terns found by earlier investigators. Experimental 


evidence is seen to be in favor of this result." * Buck- 
ling stresses calculated from the present theory were 
always lower than the ones given by reference 1. Dif- 
ference between the two theories, while small for lower 
face and core modulae of elasticity, is considerably 
larger for higher ones. Experimental results published 
in reference 6 represent excellent agreement with the 
present theory (see Fig. 14). 

t Because of the lengthy mathematical analysis involved in 
this boundary value problem, it is intended to present this part 
as a separate paper. Consequently, only the necessary formu- 
las will be given here. 
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(1) Stress DISTRIBUTION 


Consider a rectangular orthotropic elastic medium 
of length L, width ¢, and thickness unity, with the axes 
of orthotropy parallel to the x and y directions (Fig. 1) 
The differential equations of the boundary value prob- 
lem without body forces are 


A ({O°u Oru 077 
( - + (1+ A) = () 


y \Ox? oy" Ox Oy 7 
Ow Asov' O°u 

i + (1+ A) = 0} (1) 
Ox" u | oy" ) Ox OV 


M (7) = v (= 
1 — uwv\G — =f 


Eqs. (1) are solved'* under the following mixed bound- 


ary conditions: at y = 0, 
u(x, 0) = u,(x, VO); v(x, O) = v,(x, O) 
ity= 
gtx, f) = u(x, t); wx, t) = v(x, #) 
tx = 0 (2 
o,(0, v) = a,,(0, y); 7(O, vy) = v,(0, y) 
atx=L 


v(L, vy) = w(L, y) | 





N E, 


1952 


SCIENCES—JU 


where subscript o denotes the given functions at the 
boundary. 
The result of the analysis in the special case of 


u,(x, 0) = u(x, t) = v,(0, y) = v,(L, y) = 
uy (0, y) = uw (Ly) =O (3 
E,/2G >ut (u vy)’ (4 
are 


oe 2 j mr 
= - <| — U,(m, y) + 
eater i] ™ Um, 3 


- mr : 
yV, y(m, VY) | Vos(m, 0) + | — ] Um, vy) + 


; ) . mmx 
vl, (mt, y) | vs(m, t) > sin 
f L 


gE 2 > j 7 
Ty = ; > (Mm, Vv) — 
= Mv L m 1 ; 


mr. ; 
Mo Ug(m, vy) | vos(m, 0) +1 V2, - ie 


nur U ) . mm 
m (m, Vv) | v.(m, t) > sin 
Zz . { ib 


Lim, vy), Vim, v) are known 


pages 52, 


where L,(m, vy), V,(m, y), 


functions (see reference 13, 55, and 56 or 


Navy report, pages 42-45). 
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(2) DIFFERENTIAL EQUATIONS 


Consider a sandwich column composed of a thick 
orthotropic core, for which a boundary value problem 
was previously solved,'* and the faces, as shown in 
Fig. 3. Face thicknesses ¢; and f. are assumed to be 
small as compared to core thickness ¢ so that boundary 
functions u,(x, 0), v(x, 0), u,(x, ft), v(x, t) may be taken 
to coincide with the deformations of median lines of the 
faces. 

The equilibrium condition of faces that are subjected 
to axial forces P;, Ps, and reaction stresses o, and rf 
occurring at the contacting surfaces of the faces and 
the core due to presence of the core leads to the foliow- 


ing differential equations: 


d*y,(x, 0) . d*v,(x, 0) 


(EI); : + P, Saree ee ee. ay 
dx} dx? 
t, dr (x, 0) 
$ dx 
(6) 
- QUA, B) dv,(x, t) 
(ET) + Ps —— + @,(x, t) — 
dx4 dx? 


to dr (x, t) 


2 dx 


The notations and the sign convention are shown in 
Figs. 3 and 4. 

Since it is assumed that the faces are thin, the vari- 
ations of the shear and normal stresses across their 
thickness are negligible, and the effect of shear stress 
on the curvature can be neglected. 

In the problem of buckling it is known that the de- 
formations u,(x, 0), u,(x, t) are higher order quantities 
than v,(x, 0), v(x, ¢). In addition, the end support con- 
dition requires that v,(0, y) = v,(L, vy) = 0, since the 











i 


friction between ends of the core and the load applying 
device prohibits any deformation in y direction. Fur- 
ther, u, -(0, vy) and u, -(L, y) can also be taken equal 
to zero, since, as a result of ;t; > E,t, 

a,(0, y) = ox(L, y) =O, v(0, vy) = o(L, y) = 0 


and Hooke’s Law reducing the boundary condition on 


stresses in Eqs. (2) to 
Uy, (0, y) = uo, AL, y) = 0 


Consequently, Eqs. (5) can be used to compute the 
boundary stresses ¢, and 7 needed in Eqs. (6). 


Substitution of o,(x, 0), 7(x, 0), o,(x, t), and r(x, t) as given by Egs. (5) into Eqs. (6) leads to 


d4y,(x, 0) d*v,(x, 0) Poe E 
(EI) +P + = 
, dx4 dx? i, >» 1 — py 


[Vo, y(m, O)v,.(m, 0) + Vi, y(m, O)z,.(m, t)| + 


Gh ff a mr 0) +4U 9 ; mn) Mw 
: a m, VO) Vos(m, V) ,(m, O)v,.(m, t) > sin = 0) 
alii L = EF? L 
, | (7) 
_ d4y,(x, t) i d7v,(x, t) : a he , : 
(ET). + P,» = + > [V,. y(m, t)dos(m, 0) + Viz y(m, t) v%s(m, t)| + 
dx4 dx? Leimaz1 \l — wv . 
Gte Vr. Ne. o+lu ) + mr (mr ) _ mrx 9 
1 U,. (m, t)vos(m, 0) (m, t) Vo3(m, t) sin = 
a i ia L iL L 
Eqs. (7) are solved for simple-support conditions at x = 0 and x = L—namely, 
v,(0, 0) = 2,(L, 0) = », --(0, 0) = 2, --(L, 0) = 0 (face 1)} 
‘ ; S 
v(0,t) = v(L,t) = Vi, (0,0) = %, (ZL, t) = 0 (face 2) 
Eqs. (9) express the fact that the deflection and the bending moments at the supports are zero. 
(3) SOLUTION 
The definition of finite Fourier sine transform is 
*1 
v,(m, y) = y) v(x, y) sin (max/L) dx eo = x2 = L) i) 
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I 


Application of the finite Fourier sine transform to Eqs. (7) and the use of Eqs. (8) lead to two homogeneous linea 


equations in v,,(m, 0) and z,,(m, ¢). 


Wiu(m)v,.(m, 0) + Wio(m)z,.(m, t) = Ol (10 
Wo1(m)v,.(m, 0) + Wao(m)v,.(m, t) = OF 
where 
Wilm) = E\hi(mr/L)* — Py(mr/L)? + (Gt,/2) (mr/L) [U,. ,(m, 0) + (mr/L))] + 
E,/(1 — uv) |e, ,(m, 0) 
Wie(m) = (Gt/2) (mr/L) U, ,(m, 0) + [E,/(1 — wv) | V2. ,(m, 0) 
(1] 
Woa(m) = (Gte/2) (mr/L)U, ,(m, t) + [E,/(1 — uv)] Vo, y(m, t 
Wo(m) = Eelo(mar/L)* — Po(mmr/L)? + (Gte/2) (mr/L) [U1 ,(m, t) + 
(mr/L)| + [E,/(1 — wv) | Vi y(m, t 
To find the nonzero solution of these equations, the determinant of the coefficients is set equal to zero. Hence, 
P,P, — BP, — B,\P; — C= 0 (12 


Here B,, B,, and C are expressions that are functions of physical and geometrical quantities related to the faces 
and the core. 


In the case where P; and P, are both compressive loads, ¢,, is the critical value of stress o defined by 


P, = oti, P, = ote (13 
When the faces have equal thickness and are made of the same material, /,/; = /2,/,. Thus, Eq. (12) can be sim- 
plified to 
ao. = (1/2)E,(mat,/L)? + (1/2)Gu,(m, 0) + }[E,/(1 — wv) ]? (L/mat,)? [22(m, 0) — v27(m, 0)] + 
(G 2)2u2(m, 0) | ' (14 


provided F,/2G > uw + (u/v) ’*, where 


l l ‘ 
u(m, O) = 1 + ¢ + dvo(m, 0) [3 — ') sinh 2;/ sinh ct (sinh 2,¢ — 7 sinh 22f) 
é J 


l 


u»(m, 0) = dv.(m, 0) (cosh 2,¢ — cosh 2f)(sinh 2;¢ — 7 sinh 22f) 


v(m, 0) = v(m, 0) (sinh 2f cosh gf — 7 sinh 2f cosh 2;/)(sinh 2,¢ — 7 sinh 20f) 


; os iy. iy os : , 
vo(m, O) = e(sinh z,¢ — 7 sinh 22f) [ — cosh 2,f cosh got + - (; + ') sinh 2;/ sinh st 92 
‘ é ] (i) 


21 t i 


c = (1/2) (1 + A)A*[(ai? — A2)—! + (ee? — A2)—4], = w)l+[1 — (Q/w)*) °} 


ta 


d = (1 + A) (mm/L)2,(2,7 — A*)—',, w? = [(E£,/2G) — pw]|(mr/L)? 
e = (1/2) (2:2 — 2) (2:2 — A?) [(1 + A)? (ma/L)d22,]-!, 24 = (u/v)(mr/L)? 


7 = (2/22) (2.2 — A) (232 — A2)—!, A? = (p/A)(ma/L)? 


A similar procedure can be used to obtain o,, for the case of orthotropy represented by /,/2G < uw + (u/v) 
In this case, however, new expressions are needed for a, and 7 to be used in Eqs. (6), since stresses and strains in the 
core calculated on the basis of this inequality are different from those given by Eqs. (5). Here, only the resulting 
o-, is given [see Navy report, footnote, page 1, reference (13) ]. 


Ser = (1/12) E,(mat,/L)? + (1/2)Guy'(m, 0) + S[E,/(1 — wv)? (L/mat,)? [0,'2(m, 0) — v2'2(m, 0)] + 
J tl J J L J 

vY /é 9 ,9 U's: > 

(G/2)* u2"*(m, 0); (16 


provided E£,/2G <u + (u v),° where 











linear 


(10 


(11 


(15 


e sim- 


vy) a 
n the 
iting 


(16) 








RIPPLE-TYPE BUCKLING OF SANDWICH COLUMNS $13 


uy'(m, O) = 1 + c’ + d’v2'(m, 0) (sinh? z,’t + sin? 22’t) (sinh 2,'t cos 2.'t + 7’ cosh 2,’t sin 22’t)~! 

us’(m, 0) = h’ve'(m, 0) (sinh 2;’¢ sin 22’t) (sinh 2,’f cos 22’t + 7’ cosh 2't sin 22’) 

v%1'(m, 0) = (1/2) ve’(m, 0) (sinh? 2;'t — 7’ sin? 22’t) (sinh ’\f cos 22’t + 7’ cosh 2,'t sin 2.’t)~! 

ve'(m, 0) = e’(sinh 2,’t cos 22’t + j’ cosh 2,’t sin 2.’t) (sinh? 2,’t — 7”? sin® 29’t) 

c’ = (1 + A) A2(w? — A?) (Q? — 2w?A? + A4)-! , 
d’ = —(1 +A) (mr/L) 2;/ (Q? — X?) (Q4 — Qw2rA? + V4)! 


e’ = 232,'(Q24 — 2r2w? + A4) [(1 + A)*A2(maA/L) (Q? + X?®) 
h’ = (1 + A) (mr/L)Q? [22’(Q? + A?) J] 
j! = (21'/2.’) (Q? — A?) (Q? + A?)-! 


a’ = [(1/2)(Q? + w’)]”, Zo’ = [(1/2) (Q? — w?) 


(4) BUCKLING OF A SANDWICH BEAM SUBJECTED TO PURE BENDING 


A sandwich beam subjected to bending moments .V/ at both ends can be looked upon as a beam having one face 


submitted to compression load P; = \//(t + ¢,), while the other is submitted to tension load P, = — P,. 
By use of the relations £,/; = F2J2 and tf; = ft, and of Eq. (12), a, is found to be 
Or = —[E,/(1 — pwv)| (L/mat,)o\(m, 0) + E,/(. — pv) ]?(L/mat,)*v2?(m, 0) + (G/2)*[u172(m, 0) — 


= : —_ = mrt; ) 
us?(m, O)| + [(1/12) Ey(mat/L)?]? + (1/12) E,G I? ui(m, 0) ¢ IS 


(5) THe DEFLECTED SHAPE FOR BUCKLING UNDER AXIAL COMPRESSION 
In the following analysis it will be proved that the ripple-type of buckling occurs according to an asymmetric 
pattern :* 
(a) Possibility of Symmetric Buckling 


When the buckled shapes of the faces are symmetric with respect to the median line of the sandwich parallel to 
the x axis, the buckling is called symmetric. Therefore, for symmetric buckling, 
v,(x, 0) = —v,(x, ¢) (19) 
Substitution of Eq. (19) into Eqs. (10) gives 
Wiul(m) = Wy(m), Wo(m) = Wo(m) (20) 
Combining Eqs. (10) and (20), it is found that for symmetric buckling to occur the following equation must be 
satisfied : 
J sinh 2,¢ (cosh zo + 1) — 7 sinh Zot (cosh gf + 1) ol ’ 
é - Seg : = (21) 
| 1 — cosh 2¢ cosh got + (1/2) [7 + (1/7)] sinh af sinh zef f 


It can be seen that e cannot be zero, since 2; ¥ 2: and \* = 2,” is in contradiction with the inequality E,/2G > 
u + (u/v) ‘*on the basis of which the solution was derived. Then Eq. (21) reduces to 


(2:7 — A?) sinh 2,t/2, (cosh af + 1) = (2? — A?) sinh 2ef/z. (cosh got + 1) (22) 

since 7 ~ 0, 2; and 2: are finite. Let 
F(z) = (2? — A?) sinh gz ¢t/z (cosh 2f + 1) (23) 
It can be seen that F(z) is a monotonically increasing function. Suppose the contrary is true, then dF/dz = 0 
at some points. Hence, sinh 2f/zst = (A*/s* — 1)/(A*/s* + 1). This, however, is impossible since sinh zt > at 


and (7/27 — 1)/(A*/z? + 1) < 1. Furthermore, z¢ > 0. Therefore, as was to be proved, dF/dz > 0. Hence, 
F(z) > F(z) or Eq. (22) cannot hold. 
Hence, it may be concluded that symmetric buckling is not possible. 


* Proof is given here for orthotropic cores satisfying E,/2G > w + (u/v)'/?. A similar method may be used for cores satisfying 
E,/2G < p+ (p/v)'”? 
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(b) Possibility of Antisymmetric Buckling 


An antisymmetrically buckled shape is obtained when 








V(x, 0) = v,(x, t) (24 

Following the same procedure as that used in the previous discussion, it is found that 
sinh 2,f(cosh gf — 1) — 7 sinh ef (cosh gt — 1) - 
= 0 (20 


1 — cosh 2,¢ cosh Zo + (1/2) [7 + (1/7)] sinh 2;¢ sinh 2of 7 


This equation is satisfied only for = 0. Otherwise, with ¢ ¥ 0, it can be seen that the following must be true: 


F(z,) = F(z.) (26 
where 
F(z) = (2? — i?) sinh 2t/z (cosh zt — 1) (27 
F(z) is a continuous monotonic function in 0 < zs < If F(z) is not monotonic, then dF/dz = 0 at some points 
in the interval. This leads to 
(2S 


sinh 2t/zt = (1 — A?/27)/(1 + A?/s?) 
Eq. (28) cannot be satisfied except for ¢ = 0, for which Eq. (25) must be examined. A limiting process indicates 
This, however, is obvious since, in this case, the sandwich beam becomes a 
Hence, it may be concluded that anti 


that Eq. (24) is satisfied for ¢ = 0. 
beam consisting of two faces without core which fails by Euler buckling. 
symmetric buckling is not possible for a sandwich beam. Therefore, the buckled shape pattern must be an asymmetric 


one. 


The assumed core and face materials and their proper 


(6) NUMERICAL CALCULATIONS AND DISCUSSION 
A comparison between 


ties are given in Figs. 5 to 13. 


Numerical examples are presented in order to study the present theory and that of reference | is shown in 


the effects of variations in the elastic constants of the Fig. 14. 

core: ,, £,, G, v, 4, Young's modulus of the faces, the Some hypothetical materials are chosen with assumed 
elastic properties as given in Figs. 5 to 13. It is be- 
lieved that materials similar to these may be produced 
by using certain synthetic resins in combination with 


core thickness, and the face thickness. 
The results of these calculations are compared with 
previous theoretical and experimental results. 
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Fic. 7. Buckling of sandwich made of core No. 1.3 and 24ST 
Alclad faces. E,/2G > p+ (p/v)!“?. 
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Fic. 8. Buckling of sandwich made of core No. 1.3 and 75ST 
faces. E,/2G> pw t+ (p/v)!72. 


nylon fibers arranged in different patterns [see reference 
14 and Fig. 2 (a and b)]. These materials may be 
classified into two groups according to the present theory 


as: 


(a) Elastic materials with elastic properties satis- 
fying the inequality 


E,/2G > w+ (p/v)'” 
(b) Elastic materials satisfying the inequality 


E,/2G <u t+ (p/v)' 
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Fic. 9. Buckling of sandwich made of core No. 1.3 and 75ST 
faces. E,/2G> wt (p/v)!? 
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num-alloy faces. E,/2G < p+ (p/v)!/?. 


Both cases are studied here. The general conclusion 
is that the effect of variations in E, is found to be the 
most important (see Figs. 5 and 13). 

Young’s modulus of the faces is also an important 
factor, but, as can easily be seen, it appears only in the 
first term of the formula for o,, which represents Euler 
buckling of each face [see, for instance, Eq. (14)]. 
Hence, it becomes more important as the wave length 
decreases. 

The effect of variations in the core thickness is found 
to be small (see Fig. 9). 

The parameter (mrt,/L) is found to be an important 
factor. Plots of o,, vs. (mmt,;/L) always show a mini- 
mum (see Figs. 5 to S and 10 to 13). 
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Fic. 12. Buckling of sandwich made of core No. 3.2 and Papreg 


faces. E,/2G = wp + (p/v)!/?. 


A comparison between the present theory and test 
results published in reference 1 is shown in Fig. 14. 
Good agreement is indicated between the present theory 
and the tests. The discrepancy does not exceed 6 per 
cent. The deviation between the two theories is, in 
general, within 30 per cent of the values given by the 
present theory for the materials investigated. The 
present theory always gives lower buckling stress. 

The theory of reference 2 agrees well with the experi- 
ment for bending of cantilever beams, while a larger dis- 
crepancy is noticed in case of buckling. 

The method of reference 3 assumes either symmetric 
or antisymmetric buckling, both of which are proved 
impossible in the strict sense. However, in experiments 
it is often difficult to detect the differences in buckled 
shapes of faces. 

Reference 4 considers the Euler type buckling under 
initial stresses. 

The method of reference 5 may be used for sandwich 
plates and beams when the normal stresses across the 
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Fic. 18. Comparison of o.- for pure bending and compression, 
Core No. 1.1 and 75ST faces. E,/2G > w+ (p/v)!”2, 
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While the as- 
sumptions made lead to satisfactory results in the study 


thickness of plate are not important. 


of overall bending and buckling, they preclude con- 
sideration of the wrinkling phenomenon. 
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RIPPLE-TYPE BUCKLING 

In Fig. 5 of reference 7 a comparison is shown be- 
tween test and theory which takes the initial waviness 
into account. However, it is noticed that all experi- 
mental points are condensed around the yield stress of 
the aluminum alloy 75ST faces of the specimens. There- 
fore, material failure rather than instability may have 
occurred. However, initial waviness may be important 
for extremely thin sandwiches having strong cores for 
which Euler type of buckling occurs. 

More recent papers are concerned primarily with the 
Euler type of buckling of plates*~" for which only con- 
sideration of average stresses across the core seems to be 
sufficient to express the overall buckling phenomenon 
reasonably well for engineering approximation. This, of 
course, eliminates the type of lengthy analysis given in 
reference 13. However, ripple-type of wrinkling is 
excluded because of this oversimplification. Therefore, 
a correct buckling stress for plates having rather soft 
thick cores cannot be given. While an analysis similar 
to that used in reference 13 would be tedious, if not im- 
possible, the method of reference 12 avoids this difficulty 
with reasonable engineering accuracy. 

Finally, it may be of interest to mention that the 
method of the present paper may be applied to inves- 
tigate ripple-type buckling of sandwich cylinders."! 
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A High-Speed Schlieren Technique for 


Investigation of Aerodynamic Transients’ 


WALTER S. BRADFIELD? ano WALTER Y. FISH? 
Unwersity of Minnesota 


SUMMARY 


Development ot a simple high-speed schlieren technique for 
investigation of aerodynamic transients is described, and sample 
photographs are shown. A spark light source was developed 
which was capable of delivering bursts of up to 250 photographs 
at a frequency as high as 16,000 pl otographs per sec. A simple 
inexpensive camera for recording the data photographically was 
devised and constructed in the laboratory. The photographic 
units were adapted to the schlieren equipment that was on hand, 
and data were recorded at a Mach Number of 3.5 ia an inter- 


mittent wind tunnel. 


INTRODUCTION 


k bu ROLE OF TRANSIENT PHENOMENA in supersonic 
flow has grown in importance during the past few 
years. The advent of 
engine, flight at supersonic speeds, and the widespread 


the aerothermodynamic jet 


application of the supersonic wind tunnel to research 
problems have stimulated the modification and further 
development of basic optical techniques to this end. 
The technique described herein was developed in con- 
nection with a study of wind-tunnel starting transients, 
which is sponsored by the Air Materiel Command under 
contract AF 33 (038) 7078.° 

Specifically, it was desired to study nonstationary 
shock-wave propagation and reflection and the effect of 
these transients on the initiation of flow in a ‘‘choking”’ 
supersonic diffuser. The whole starting process was 
expected to take place in less than 0.020 sec. It was 
considered necessary to obtain at least 30 to 50 schlieren 
photographs during this interval in order to adequately 
define the flow. With the instrumentation completed, 
150 photographs were easily obtained during the speci- 
fied period (corresponding to 7,500 frames per sec.), and 
frequencies as high as 15,000 frames per sec. were re- 
corded. The exposure time is roughly 2 microsec 

SPARK LIGHT SOURCE 


The light source developed to photograph the phenom- 
ena discussed above grew out of the basic single spark 
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unit using the Libessart gap as the actual spark head, 
An undesirable feature of these single-shot units proved 
to be the link that led to the present design—i.e., when 
the supply potential was raised beyond a critical point, 
the gap broke down and fired sporadically without a 
trigger. 

The design requirements as indicated by the above 
discussion can be summarized as follows: 

(1) A spark repetition rate upward of 7,500 per sec. 

(2) Sufficient intensity for use with 
apparatus. 

(3) Duration of approximately 20 millisec. per burst 

(4) Individual short 
assure sharp exposure of rapidly moving phenomena. 


schlieren 


spark duration enough to 
(5) Ability to control beginning and end of spark 
train or burst. 
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Fic. 2. Schematic diagram of the repetitive spark light source 
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INVESTIGATION OF AERODYNAMIC 


These requirements were met using the system to be 
described ; a simplified diagram (Fig. 1) is included for 
purposes of discussion. . 

Basically, the repetitive spark light source consists 
of a power supply of high storage capacity, a switch, 
istor-condenser components of the charge-discharge 


res 
close-spaced gap where the electrical 


circuit, and a 
discharge energy 1s converted into light. With the 
input voltage adjusted at some point higher than the 
critical (breakdown) gap voltage, the switch is closed and 
C charges up toward the source voltage; more, however, 
when the critical voltage is reached, the gap becomes 
ionized (breaks down) and the condenser is rapidly 
discharged through the gap until the voltage (/) can 
no longer sustain an are and ionization ceases. The 
condenser then recharges at a rate determined by R-C 
and the sequence is repeated. The arc, or spark dura- 
tion is determined by the value of C and the resistance 
of the ionized gap. Thus, the repetition rate is estab- 
lished by the R-C constants and input voltage level, 
while the duration depends on C and R,. 

The actual circuit used is shown schematically in 
Fig. 2. The power source supplies 0 to 15,000 volts at 
5300 milliamp., and its storage facilities are augmented 
by the 2 mfd. condenser at the spark input. The addi- 
tional vacuum tubes and components are necessary 
to control the starting time and duration of the burst. 
To describe a typical operating sequence (with SW-3 
closed and the thyratron, 2D21, cut off), the hydrogen 
thyratron switch tube 4C35 is held nonconducting by 
the negative battery bias. The 6J6 relay tube is con- 
ducting and the relay, Rly 1, is closed. Consequently, 
the supply voltage would be applied to C and the gap, 
except for the 4C35 switch. At the moment that it is 
desired to start the spark pulses, SW-3 is opened and 
the network R,; and R, becomes a voltage divider that 
applies a sudden positive bias to the grid of the 2D21, 
through C-2, as a pulse. Firing of the 2D21 trigger 
tube causes the 4C35 to ionize, thus allowing full 
application of the source voltage to the gap and con- 
denser C through charging resistor R. This action 
initiates the repetitive sparks, which continue until 
the 6J6 (relay tube) is cut off by action of the positive 
cathode voltage and at a time later determined by C,, 
the 840-ohm cathode resistor, and normal relay inertia. 
When the relay opens, the high-voltage circuit is broken 
and the sparks cease. The condenser is shunted across 
SW-1, and the 4C35 absorbs the transient and prevents 
arcing in the circuit. 

It has been determined experimentally that the spark 
duration is in the neighborhood of 2 microsec.; this is 
sufficiently short to ensure good exposures. 

In early experiments with the apparatus, it was 
found difficult to obtain bursts of a high repetition rate 
for any appreciable length of time (with respect to the 
desired 20 millisec.) without the gap becoming com- 
pletely ionized. This difficulty was overcome by piping 


in air, under pressure, to the enclosed gap area or 
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Fig. 3. It is assumed that the air pressure drives 
the ionized particles from the enclosure, thus allow- 
ing the gap to deionize and extinguish between 
sparks. 

The mechanical construction of the spark head and 
associated components is such as to facilitate assembly 
of all high-voltage parts with minimum lead lengths and 
also to allow its installation on the optical apparatus 
(Fig. 4). 


CAMERA 


It was necessary to devise a camera that would per- 
mit recording photographs at the frequencies anticipated 
and which would be adaptable to the schlieren optical 
system planned for use with the operation. Since the 
time duration of the process under investigation was 
relatively short, it was found possible to utilize fairly 
short strips of film for recording purposes. This sug- 
gested immediately the possibility of using a simple 
drum camera. The camera constructed consists of 
an aluminum drum 20 in. in diameter which is rotated 
at a constant speed by an electric motor. A drive from 
the motor to the camera is effected by a pulley system 
in which driving pulleys of variable diameter permit 
continuous variation of the drum speed between 2,000 
and 4,000 r.p.m. (Fig. 4). The drum is enclosed by a 
plywood box (Fig. 5), and, in order to facilitate loading 
the film during the daylight hours, a light, tight, felt 
loading bag was attached to the top of the box. In 
loading, the film is fastened to the outer perimeter 
of the drum by means of double-sided Scotch tape. 
Standard 35-mm. film is cut to length for this pur- 
pose. 

The objective lens assembly was made up of avail- 
able optical components as shown in Fig. 6. It con- 
sists of one achromat element and one meniscus element 
arranged as shown. The focal length and the image to 
object ratio are variable, although the assembly was de- 
signed with a 25-ft. object distance in mind. Since 
the lens system projects an image as small as 4 by 6 
mm., the camera will accommodate spark frequencies 
to 16,000 per sec. 

A simple mechanical shutter has served to protect 
the film from fogging, and for this application the tim- 
ing of the cycle is known with sufficient accuracy to 
permit the manual operation of the shutter. The 
shutter remains wide open for the duration of the aero- 
dynamic phenomenon. The shutter has been left 
open for periods as long as 3 min. (with the drum run- 
ning) without objectionable fogging of the film. 


SCHLIEREN OPTICAL SYSTEM 


A standard double-pass schlieren optical system was 
adapted for use with the spark light source and drum 
camera. <A diagram of the optical setup is shown in 
Pig. 7. 
0.060 in. in diameter, and it emits lightasa conical beam 


The spark light source is of approximately 
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Fic. 7. Schlieren optical diagram 














Fic. 8. Drum camera schlieren series showing flow and P- 
shock passing over 10.88° subcritical diffuser; 7,000 frames per 
sec. 


due to the shape of its negative electrode. The beam is 
picked up by the condensing lens, which produces an 
image of the light source on the aperture of the schlieren 
optical head. A right-angle prism then directs the beam 
of light through the tunnel test-section window be- 
hind which is situated a spherical mirror whose con- 
jugate focus is located at the knife edge of the schlieren 
system. Through another set of prisms in the ob- 
jective head, light is passed into the objective lens 
system of the drum camera. 


APPLICATION AND RESULTS 


A schlieren system of this particular type is extremely 
useful in its application to nonstationary supersonic 


(Continued on page 432 





On t 


John V 
Fulbrig 
Februa' 


where 


Thi 
forma 
respec 
of the 
soluti 
ner-H 
monic 
form « 
tion* | 

An 
in as 
poten 


in th 
been ' 
treatt 
the h 
the o 
gener 
made 
know 
is qui 

We 
is bo 
y+ 
the I: 


=O 


P. 
per 


an 
en 
im 








Readers’ Gorm 





RIEF REPORTS of investigations in the aeronautical sciences and discussions of papers pub- 


lished in the JOURNAL are presented in this special department. 
pleted as quickly as possible after receipt of the material. 
itself responsible for the opinions expressed by the correspondents. 


Publication will be com- 
The Editorial Committee does not hold 
Contributions should not exceed 


800 words in length. 


On the General Solution for Unsteady Motion of a Rectangular Wing in Supersonic Flow 


John W. Miles* 


fulbright Lecturer, Auckland University College, Auckland, New Zealand 


February 28, 1952 
W: CONSIDER A WING TIP projecting on the quadrant z = 0, 


x > 0, y > 0 in a supersonic flow of Mach Number  di- 


| "sg mB \(x-—t)+y 
dé 
wf 0 B-\(x-§&)—y 


x ea 


a(x, y, t) = 


2 


(x -¢ 


where a is the sonic velocity and 


(2) 


This solution was obtained by carrying out a Fourier trans- 
formation with respect to ¢ and a Laplace transformation with 
respect to x and comparing the result with the Laplace transform 
of the corresponding problem in steady flow, for which a general 
solution is available (either by Evvard’s method! or by the Wie- 
ner-Hopf approach?). This general procedure, at least for har- 
monic time dependence, is due to Magnaradze,* but the general 
form of the final result was suggested to us by Stewartson’s solu- 
tion’ for w of the form f(x)e"™ 

An alternative approach to the general problem is to expand w 
in a series of terms like y"w,(x, t), for which the corresponding 
potentials may be derived from solutions of form 


Q\rt nee 1 ; 
¥, = ar ~ cos|{ a +-—]0| F,(x,4,r) (3) 
oy 2 


in the cylindrical polar coordinates (x, r, 6). This method has 
been applied by Rott® in the special case nm = 0, following Lamb’s 
treatment of the half-plane diffraction problem.® Rott discards 
the higher order solutions (x > 1) of form (3) in consequence of 
the order of the singularity exhibited at r = 0, but for the more 
general case ¢,, as obtained by integration of Eq. (3), can be 
The F, (x, r, t) may be determined from a 
‘x, but this 


made to vanish as r /?, 
knowledge of the potential on the Mach cone r = 8 

is quite involved for n # 0. 
We remark that the domain of integration for (£, 7) in Eq. (1) 
is bounded by the leading edge and the pseudo Mach lines 7 = 
£)/8| cos 6 and y [(x — £)/B8] cos 6, except where 


the latter is off the wing, in which case it is to be replaced by its 


y+ [(x — ¢) 


* On leave from University of California at Los Angeles 


wilt, n, £ — (M/Ba) (x — &)| dn 
— By — n)?]' ; 


rected along the x axis. If the prescribed, normal velocity of the 
wing is w (x, y, t), positive down, a general solution for the velocity 


potential on the upper surface of the wing is 


1 "3 *«/2 ee 8B-'\(x—£)cos 0 + y 
ner / dy / de dé x 
il 0 0 0 8B~'(x—£)cos @-—y 


welt, n, t — (M + sin @cos y) (x — §£)/B*a] dn ' 
(1) 
[(x — &)? cos? @ — By — n)?] 
reflection (in the side edge) » = [(x — £)/8!cos@ — y. At the 


lower limit 6 = 0, these lines delimit the “equivalent area’’ pre- 
scribed by Evvard,' and the solution appears as a generalization 
of Evvard’s result for a rectangular wing. Unfortunately, this 
interpretation has no direct extension to more general plan forms, 
where the lower limit of integration for n depends on x and é inde- 
pendently [rather than simply on (x — £)], in consequence of 
which the Laplace transform may not be effected by convolution. 

Finally, we remark that the result [Eq. (1)] can be extended 
to any wing having a supersonic leading edge of arbitrary shape 
and straight subsonic side edges (not trailing), provided that there 
is no interference between the latter. To effect this extension we 
transform (via a modified Lorentz transformation) the potential 
equation to the normalized wave equation 


de’! + byy + bez = Gr'x’ (4) 


in which form it is invariant after a (true) Lorentz transformation 
with respect to (x’, y). Such a transformation can be used 
(cf. the oblique transformation of Jones,’ as modified by Hayes* 
to form a group) to bring the port edge to y = 0, while the super- 
sonic leading edge can be filled in to x’ = 0 on the proviso that 
0 over the filled-in portion. We hope to discuss such ex- 


tensions in more detail in a subsequent paper 


we 


ws 
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the stability of a system without the necessity of finding all th 5 
roots of the characteristic equation (see, for example, referent aft 
1). Briefly, it states that, for any function f(p), B 
Z=N+P 7 ee 
W 
+ The R. T. Jones approximation is used here to represent the Wagner p 
Function—i.e., | (s) = 1 0.1656 ~%"4 0.335¢ ~°% see Referer 5, J 
t Strictly speaking, the criterion used here is the more basic Routh c a | 





terion, since it involves encirclements of the origin rather than of the 


Fic. 1. Nyquist plot for bending-torsion flutter point 
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where Z and P are the number of zeros (roots) and the number of 
poles located in the positive real half of the complex plane, respec- 
tively, and where N is the number of encirclements of f(p) about 
the origin as p maps the boundaries of the positive real half of the 
complex plane. For the system given by Eqs. (5), the character- 
istic equation /(p) will be of the form 

ip) = Ko + Ki(1/p) + Ko(1/p?) + Ka(1/p*) + 
K(1/p*) (8) 
It is to be noted, though, that some of the multipliers K,, will 


contain the polynomial approximation to the Theodoresen C(k) 
function 

; 0.5p? + 0.2808p + 0.01365 

( (p) = one _ (9) 

p? + 0.3455p + 0.01365 

which enters from the aerodynamic coeflicients of Eqs. (6). How- 
ever, since the roots of the denominator of Eq. (9) are both neg 
itive real numbers, the characteristic equation f(p) will contain 
no poles in the positive real half of the complex plane. * 

The basic system will therefore be stable if there are no net en 
circlements of f(p) about the origin as p maps the positive real 
half of the complex plane. 

To apply this Nyquist criterion, a value of velocity v is first 
set into Eqs. (5). The determinant is then evaluated for various 
values of p = 1k ranging from p = + 1 to p = 10%, especially 
in the vicinity of w near a natural frequency. The counterclock 
wise circuit of p about the positive half of the real plane is com 
pleted by continuing in the conventional manner with due regard 
for passing from p = i0* to p = iO” and from p = i to 
It is to be noted that for this problem, f(p) — 


b= +71 


K(1/p*) as p — 0, thus resulting in two counterclockwise en- 
circlements of the origin in passing from p = 10* to p = 10 It 
should be further noted that the value of the determinant of Eqs 


for p = ik will be nearly identical to that of Eqs. (2), the only 


difference being that Eqs. (5) use the R. T. Jones approximation 
for the C(k) function, while Eqs. (2) utilize the exact expression 
in terms of Bessel functions.? Hence, in plotting the Nyquist 
diagram, the more familiar sinusoidal form [Eqs. (2)] may be 
used directly, the less familiar Laplace operator form being used 
only to investigate the behavior when k approaches zero. 

The preceding analysis was applied to a three-dimensional two 
degree of freedom bending-torsion case. The flutter speed, as 
worked out by the standard A.M.C. method, resulted in vp = 
855 m.p.h. Three Nyquist plots, two for the stable 400 and 800 
m.p.h. speeds and one for the unstable 900 m.p.h. speed, are 
shown in Fig. 1. 

The Nyquist method of solution is convenient when there are 
many degrees of freedom in the flutter analysis. The method can 
also be used in synthesizing servo systems where aeroelastic 
flutter characteristics must be included, since these character- 
istics can be introduced as the airplane block of a more compli- 
cated automatic control system block diagram. However, it 
should be mentioned that, although stability or instability is defi- 
nitely determined, the degree of stability is found only in a lim- 
ited sense. Also, a range of different forward speeds v must be 
investigated to ensure that the system is stable throughout the 


flying range of the iirplane. 
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Comments ‘‘On the Deflection of Swept 
Cantilevered Surfaces’’ 


David B. Hall 
Structures Department, The Glenn L. Martin Company, Baltimore 
February 26, 19592 


I THEIR PAPER, “‘On the Deflection of Cantilevered Surfaces,”’ 

in the December, 1951, JoURNAL OF THE AERONAUTICAL 
ScIENCES, Professor Martin and Mr. Gursahaney have reported 
“effec 


tive elastic axis’’ and have expressed the belief that this will also 


finding experimentally that a swept-tapered plate has an 
be the case for a typical airplane wing. The following analysis 
indicates that the latter is unlikely. 

Essentially, a swept cantilever of high aspect ratio is an ordi 
nary cantilever with a skewed support which produces localized 
stress disturbances. The stresses in the vicinity of the support 
may cause the beam, under bending or twisting moments, to ro- 
tate about other than the bending or twisting axes. One of the 
requirements for the existence of an “effective elastic axis’’ is 
that in the normal portion of the structure there can be a family of 
lines such that 
leave that line straight; some one of them should remain hori 


a twisting moment about any one of them will 


zontal also 


Consider a symmetrical beam (Fig. 1) having section properties 


as follows: 


The deflection of a point on a line y = mx will be 


Me dx M, ” dx 
w= dx + mx 
C s* C2 x" 


Mp lie My x? 
_ + ml T 
aq (1 —n)(2 —n) C2 l n 


(linear terms 


To fulfill the definition of an ‘effective elastic axis,"’ we must have 
w = (linear terms only) along a line y = mx, and we must have 
Mz = mMy. Making these substitutions, we find 
Q/c = GJ/EI = 2 n (1 
For a trapezoidal plate nm = 1, and GJ should be equal to E/.* 
Actual values are approximately as follows: 
EI = Ebd? 


GJ & Gbd3/3 = 0.375 Ebd?/3 = 0.125 Ebd 


12(1 — uw?) = 0.094 Ebd* 


These values are not exactly equal but may be near enough to 


account for the results found experimentally. However, for 


* Eq. (1) applies even through the equation for w assumes a special form 


forn = 1 
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AXIS OF SYMMETRY 
x= —- - 
LINE _y = mx —g 
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TYPICAL AXES 











GJ = 2EI. 


typical wing structure, in which n is approximately 4, the cur- 


vature due to twist reverses, and we find that Eq. (1) demands GJ 


and El of opposite signs or, more realistically, that under an in- 
clined torque the only straight line on the wing is inclined the op- 
posite way. Thus, the existence of an ‘‘effective elastic axis’’ ap- 
pears to be a unique property of a special kind of beam of which 
the tapered plate is an example. Another interesting example 
would be a prismatic beam so constructed that GJ = 2EI. Such 
a beam could be twisted about a double infinity of axes (Fig. 2) 
and, when cantilevered from a skewed support, could have a 
single infinity of effective elastic axes. 


On an Extension of Slender-Wing Theory 


Mac C. Adams and W. R. Sears 


Graduate School of Aeronautical Engineering, Cornell University, 
Ithaca, 


February 24, 1959 


— FAMILIAR SLENDER-WING THEORY of Munk and Jones! 
may be interpreted as the result of expanding a solution of 
the three-dimensional small-perturbation equations in powers of 
a width parameter and retaining only the leading term. In this 
note, we shall show how an expansion of this type can be ob 
tained for thin wings in supersonic steady flow. Our procedure 
makes use of Laplace-transform methods, previously applied by 
Ward,? Stewartson,* and others. It gives us higher order terms, 
which appreciably extend the scope of the slender-body theory 
Consider first the symmetrical (nonlifting) slender wing 
sketched in Fig. 1. Let x, y, z be orthogonal Cartesian coordi 
nates, as indicated, and let corresponding dimensionless coordi 





u | M>t 





J 


y= sh,(x/t) 
y =sh,(x/2) 
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nates be defined by — = x//, n = y/s, and ¢ = 2/s, the Constants 
Zand s being defined in Fig. 1. If the thic bia distribution jg 
given by z = +/(x, y) the slope of the surface is ¢ = Ot/dx = 
o(f,). Further, let the plan form of the wing be defined by the 
left and right edge curves » = /,() and n = h,(£), respectively 
The differential equation satisfied by the velocity potential , 


is 
Bk pez = Onn T OE 


where & denotes the width parameter s// 
In terms of the Laplace transform of ¢, which we shall call 

(Pp, 0, £), Eq. (1) reads 
(Bkp)*d = dann + ee (9 


In Eq. (2) the proper first-order boundary conditions of yanish 
ing velocity components at the leading Mach cone have bee 
introduced. An elementary solution of Eq. (2) which is useful 
for the construction of nonlifting wings is 


do = f(p, m)Kol Bkpr)dn 


where r? = ¢? + (n — »,)? for any value of n; and Ko is the modified 
Bessel function of the second kind. In the &, n, ¢-space, this soly 
tion represents a line source of variable strength, lying along th 
line 7» = m in the &, n-plane. Its strength is proportional to th 
inverse Laplace transform of f(p, m), say f(£, m), and a continuous 
distribution of such line sources in the é, 7-plane can be used t 
satisfy the boundary condition at the wing surface. 

Similarly, a distribution of line doublets, which is obtained by 
differentiation with respect to ¢, can be used to construct th 
flow past a lifting wing without thickness; in this case the ek 
mentary solution of Eq. (2) is 


do” = Bk pcr 'o( p, m)Kil Bkpr dn 

The relation between the distribution functions f or g and the 
surface slope or surface potential, respectively, is found by cal 
culating the behaviors of these elementary solutions at the line 
singularities; the results, in terms of Laplace transforms, are 


For nonlifting wing: 


l Fos 
f( p,m) = ornP.m, +0) = - l o(p,n 
7 


For lifting wing without thickness: 
gl p, m) = (1 t )p( p, ", TU) 


In principle, therefore, the superposition of line singularities in 
this way suffices to construct the flow past supersonic thin wings, 
either lifting or nonlifting. The velocity potential is obtained by 
integrating the effects of the elementary solutions across the span 
and inverting the Laplace transforms; the properties of the Bessel 
functions are such that only elements in the forecone of any point 
contribute to the result. In fact, the results of the inversion are 
just the familiar surface integrals over source or doublet disiribu- 
tions of linearized supersonic wing theory 

To obtain the desired expansion in the width parameter §s 
we now restrict ourselves to regions near the wing and expand 
the Bessel function Ko for small values of its argument—viz 


K,o( Bkpr) ~ —[1 + (BRpr/2)?| In (BRpr/2) 
y¥ + (1 — y) (Bkpr/2)? + O|( Bkpr) In ( Bkpr 3 
where y denotes the Euler constant. The elementary potential 


> 


function dé can now be obtained by multiplying Eq. (3) by 
f(p, m)dm and interpreting the operators term-by-term with the 
help of the convolution theorem. In the subsequent superposi 
tion by integration over 7m, there is some difficulty at the leading 
edge, but this can be overcome; the result is 
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Bk U 
P eo - J som(%) ~ ot Ie S"(t) In (& — t) X 
"2 oe — sd o? 
e+ — f eaeres -( as Jue 
* hy 
o(é, m) | (1 — y)r? — r? In (*) dn — 
of fi a he 


r2a(t,m)dm (4) 


pel) © [y + In (é — ¢)] dt 
de O& JO hi 


S(~) denotes the cross-sectional area of the airfoil at & 





where 
i.e, 
"2 
S‘(é) = a(t, n)dn 
Similarly, 
ln oF 
Wino) =— Ss as,m) [am — | eer’ 
? Jo #0 r2 2r oF 
1 ' (=) ; Bk? 93 
gs, m) | ¥ + ~ in ear | a 
» #0 Z 2 2r oF 
ly — In (E + ¢)] dt g(t, mdm (5) 
0 Jo #0 


In Eq. (5) we denote the integration across the span by the sym- 
bol “@ # 0” to remind the reader that contributions of trailing 
vortices must be included when they are present. Such vortices 
extend rearward behind all trailing edges, just as in Jones’s work. 
(4) have their counterparts in 
the last 


The first three terms in Eq. 
Ward’s results? for slender pointed (nonplanar) bodies; 
two terms are the additional ones supplied by the present theory 
for the special case of a thin wing. Similarly, the first term of 
Eq. (5) is Jones’s result, and the other two are the new ones 
that are expected to extend his theory to broader wings. 

For the direct calculation of flow about a specified lifting wing 
without thickness, Eq. (5), as usual, becomes an integral equa 
i.e., the derivative ¢-(£, n, 0)/sU must be put 
To solve 


tion for g(£, m) 
equal to the slope of the wing surface, say a(&, 7). 
for g(£, m), we find that a method of successive approximations is 
We begin by neglecting the second and third terms of the 
li.e., the terms corresponding to the second and 


useful. 
right-hand side 
third right-hand terms in Eq. (5)], and thus our first approxima 
tion to g is that of Jones’s theory, which we designate g". This 
value is subsequently substituted in the second and third terms, 
yielding a new integral equation for the second approximation, 
g), still of the same form as Jones’s integral equation. This can 
be solved immediately for g‘?. The process might be repeated 
as often as desired, but a consideration of orders of magnitude dis- 
closes that further iterations are beyond the accuracy of Eq. (5) 
itself. An explicit formula for g in terms of a(t, 7) can be 
written down by means of this procedure, but, in order to con- 
serve space, it will not be reproduced here. (The process is not 
quite so straightforward for wings for which the Jones solution 
yields discontinuous pressures on the wing; in such cases a more 
careful treatment is required in the region of the discontinuity. 
This will be given in a more complete report.) 

An interesting test case is the well-known flat delta wing at in- 
cidence. For this case, ‘‘exact’’ results of linearized theory are 
available,‘ and the results of our present theory are found to be 
identical with these through terms of the orders carried. In 
Fig. 2, for example, the lift-curve slope, in the form C;,/ak, is 
plotted against the width parameter 8/2 for ‘‘exact,’’ Jones, and 
the present theories. Moreover, for this special case, because of 
the type of dependence of the results on Bk, the curves of Fig. 2 
actually display the relative accuracies of the three approxima- 
tions in predicting local pressure coefficient and all the engineer- 
ing values that are proportional to this coefficient. The new 
method is seen to extend considerably the slender-wing theory. 
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dicted by three approximations. 


In a longer paper we intend to present the theory in more detail 
and to show its application to several typical airfoil problems, in- 
cluding the design of wings for minimum drag. ‘‘Exact”’ linear- 
ized theory has not proved tractable for such (nonconical) prob 


lems, while Jones’s theory is limited to very slender plan forms 
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Diffraction of Blasts by Axisymmetric Bodies* 


H. F. Ludlof and M. B. Friedman 
Daniel Guggenheim School of Aeronautics, New York University 


February 10, 1952 


A’ ANALYTIC METHOD developed by C. S. Gardner and H. F 

Ludloff¢t and later extended in the thesis-of L. Ting,f has 
now been applied to the diffraction of blasts advancing over 
axisymmetric bodies of arbitrary profile f(z). Besides serving a 
practical purpose, the computations were undertaken to explore 
a spatial shock configuration with regard to its Mach effect and 
its ‘‘slipstream.’”’ Suppose that the shock front is of arbitrary 
strength and originally plane and that the blast struck body is 
relatively slender. Then, during the passage, the blast will con- 
sist of a plane incident shock, a spherical weak reflected shock, 
and a curved diffracted Mach stem. The air stream that has 
traversed the Mach stem will be rotational and, for sufficiently 
strong blasts, supersonic with regard to the body. The problem 
is treated as a time-dependent perturbation, using the method of 
Lorentz transformation and ‘‘reflexion.’’}’ { But in the present 
problem, the solutions have to be assumed in the form of “‘re- 
tarded potentials.” Fixing cylindrical Lorentz coordinates (x, 
r, t) in the undisturbed air downstream of the shock, the problem 


may be stated in terms of the disturbance pressure p\'’, which is 


of 0( €*): 
Differential Equation: 


“Ap —c (1) 


2p? = 0 


* Investigation supported by Navy funds 
T See review article by H. F. Ludloff in Advances of Applied Mechanics 


Vol. IIT, 1952 
t See article by L. Ting and H. F. Ludloff, JourNat or THE AERONAUTICAL 
Sciences, Vol. 18, p. 143, 1951. 
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Boundary Conditions: (Surface) at r = 0,x <0,t>0: 
rp) = —Ao[f’ (u) + f(u)-f"(u)] (2) 


were u = a(t + Xo-x). 
(Shock Front) at x = 0,r >0,t> 0: 


D.. 1°? (p) = 0 (3) 


where D,, ,°?) is a differential operator of the second order in x 
and ¢. The parameters Ao, a, and X are simple expressions in 
terms of the shock velocity U, relative to the airfoil and the ve- 
locity U and Mach Number & of the shock relative to the un- 
disturbed flow behind. 

It is necessary to note that at the intersection point of shock and 
body surface, there is a singularity in the derivative, even for a 
conical body, characterized by 

lim (rp? = —Ao[f’ %at) + f(at)-f"(at)] 
) 


r=0, x > ( 
however, 


lim (rp-) = —Bo[f’ *at) + f(at)-f"’(at)] (4) 
x=0,r— 0° 
In accordance with the “reflexion principle,’’ Eq. (3) can be 
replaced by a condition on the “extended axis’’ of the body: 
Atr =0,x>0,?>0: 
3 
rp) = : e Ai[f’ %ui) + fl(ui)-f"(u)] (5) 
i=1 
where u; = a(t — A;x). Also, Bo, Ai, and \; are simple functions 
of Up, U, and M. 
It is now essential that, because of the properties of the wave 
equation, the solution be assumed in the form of a ‘‘retarded po- 
tential’’—integral 


l 0 nit, t — (R/c)] 
p(x, r,t) = re ' dé ail R + 
re elt, t — (R/c) 
dé a It (6) 
0 R f 
where R = [(x — £)? + r?]'/? is the real distance between source 


point (, 0) and field point (x, 7). It turns out that the ‘‘retarded 
temporary source strength’s” g; and go are equal to expressions 
(2) and (5), respectively, after x and ¢ have been replaced by £ and 
t — [R(&)/c]. Note that now only one integration is involved. 

After substituting Eqs. (2) and (5), Eq. (6) represents the solu- 
tion for a general profile form f(z). 

For the particular case of a blast passing over the conical nose 
of a missile, the solution in physical (conical) coordinates — = 
x’/ctand n = r’/ct reduce to 

: 1 t¢— M 
pY(E, n) = K-sinh™! . _ 4 
V1 — M? 7 


3 1 t—-M—a«)) 
> Ai sinh J “it (7) 
M2 n { 


i=0 V1 = 


where the w; are simple functions of £ and 7. 
The density field is given by 


pep) = (1/c?) [pM (E, n) — 2- pM M, n)] (8) 


Also, K and Q are simple expressions in Uo, U’, and MV. 

In Fig. 1 the surfaces of constant pressure and density are 
plotted for a blast passing over a cone of haJf vertex angle « = 
10°-'. The shock strength is characterized by (U» — U)/c=1.21, 
so that this plot can be compared with the two-dimensional plot 
in article by Ting and Ludloff.{ The following features emerge: 
The rotationality of the flow is again restricted to the domain 
TOC, since the reflected shock, in general, is of O(e*); at 7, its 
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strength is of even higher order. It is seen that in the present 
case, the splitting off of corresponding isopycnics and isobars is 
far less strong, and the density gradient along any isobar is Weak 
There are no longer singularities at 0 and D; in neither point qo 
several isopycnicsend. Thus instead of a contact discontinuity 
line there is, here, a domain (70C), over which the required B 
tropy change is spread. The isobars are conical surfaces jn the 
Karman-Moore region and follow this general trend, to some 
extent, also inside the time-dependent domain; apparently the 
strong singularities located on the axis of the conical body haye 


a dominant influence on the picture. 
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Effect of Wind-Tunnel Nozzle on Steady-Flow 
Nonuniformities 


H. S. Ribner and M. Tucker 
Lewis Flight Propulsion Laboratory, N.A.C.A 
March 3, 1952 


T “HERE IS AN AREA of common ground between the extensior 
of Prandtl’s concepts! reported under the above title by Pro 
fessor Corrsin? and another recent investigation.* Suppose the 





velocity in the low-speed section upstream of the nozzle, although 
purely axial, is not quite uniform across the tunnel. Then, ac 
cording to an energy consideration by Prandtl,! the nonuniform 
ity is considerably smoothed out downstream of the contraction, 
the more so the greater the contraction ratio. This smocthing | 
effect on longitudinal perturbations (but not on lateral perturba- 
tions) is still observed for the unsteady turbulent fluctuations in 
the wind-tunnel stream. Corrsin has considered the former case 
of steady perturbed longitudinal filaments, and has extended 
Prandtl’s analysis to include (a) the effect of density (~ tem- 
perature~!) perturbations and (b) the compressibility of the 
flow. We have concerned ourselves with the latter case of the 
turbulence carried by the wind-tunnel stream, utilizing the spec- 
trum concept, and have included (b) but not (a). 

In the spectral or Fourier integral representation, a turbulent 
field may be regarded as a superposition of plane, transverse, 
sinusoidal waves with an apparently random distribution of 
orientations and wave lengths. The mutual interference of the 
waves during their transport through the wind-tunnel contrac- 
tion may be neglected under suitable circumstances; then the 
contraction effect may be computed for individual plane waves 
and the overall effect obtained by superposition. 

Of this aggregation of plane waves, we now limit attention to 
those representing purely longitudinal disturbances: The axial 
velocity varies sinusoidally across the tunnel. These are then 
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of the steady type considered by Corrsin, since the convection of 
these waves by the main stream introduces no unsteadiness. It 
was found that, when the general result was specialized to such 
waves, there was obtained for any stream line 

( U; { l ) 


Us = Us)uy 


where subscripts 1 and 2 refer to sections upstream and down 
1 of the contraction, respectively; U’ is the main-stream ve- 


strean 
The same relation 


locity; and u is the perturbation velocity. 
holds for a rather wider class of waves.® 
Corrsin’s Eq. (8) reduces to the above equation for our case of 
no temperature perturbation (his @,; = 0). The methods of der- 
ivation were totally different: He employed considerations of en- 


ergy, whereas we employed considerations on the transport of 


vorticity. 
When the complete aggregation of plane waves is superposed 
applied between 





to represent turbulence the simple relation (1) 
root mean square (r.m.s.) values] is, of course, appreciably modi- 
A rough approximation to our equation® in the range 1 < 
< is 


fied 
U2/l 


us (r.m.s.) = (U/U2)"*8u, (r-m.s. ) (2) 
when the initial turbulence is isotropic (no preferred directions) 
and the wind-tunnel contraction is axially symmetric. This re- 
sult is limited to circumstances‘ where the decay of the turbulence 


during its transport through the nozzle is small in comparison. 
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Comments on ‘‘Supersonic Flow of a Two- 
Dimensional Jet’’ 


E. B. Klunker and Keith C. Harder 
Langley Aeronautical Laboratory, N.A.C.A., Langley Field, Va. 
January 16, 19592 


I Pal’s! PAPER, the mathematical formulation of the prob- 
lem of the flow of a two-dimensional supersonic jet into sub- 
sonic stream appears to be in error, since the boundary conditions 
given by Eqs. (6) and (7) of his paper are valid only for positive 
values of x. For negative values of x, the proper boundary 


conditions are 


0¢:/OyY = O¢g2/Oyv = 0 


on the boundary y = 1. No condition can be prescribed for 
0¢2/0x on this boundary for x negative, since the external flow is 
subsonic. In addition, the boundary condition given by Eq. 
(Sa) must also be satisfied atx = — ©. 

The terms of Eqs. (12) and (13) are solutions for the flow 
infinitely far downstream of the exit which satisfy the boundary 
conditions given by Eqs. (5), (6), and (7). Thus, any superposi- 
tion of these solutions will satisfy the same boundary conditions. 
As noted, however, these are not the correct boundary condi- 
tions for the flow upstream of the jet exit. Consequently, the 
complete solution cannot be determined through a superposition 
of only solutions of the type given by Eqs. (12) and (13). 

A formal solution satisfying the proper boundary conditions 
may be formulated by considering the disturbances in the sub- 


FORUM {27 


sonic field to arise from a distribution of sources on the interface 
forO<x< +o. The strength of the source distribution is deter? 
mined by equating the pressures and flow directions at the inter- 
face. 

Consider a coordinate system such that one interface lies on 
Then, in 


the axis y = 0 and the jet exit is on the axis x = 0. 


terms of source distributions, 


1 *© v(t, 0) 
u(x, 0) = * dt (1’) 
Bor _x-e 
and by the Allied integral of Fourier’s integral? 
By ” ult, 0) 
v%(x,0) = — * dé (2") 
7 x—é 


—« 


The supersonic solution may be expressed 


—- &(1 + y))] (3’) 


u(x, ¥) =f Ix + B(1 + y)) +f [x 
B iflx + Bi (1 + y)) — f[x — Bs (1 + y))} (4’) 
The boundary conditions are satisfied 


where f(z) = 0 fors < fh. 


if 
v (x, O) + v(x + 28,0) = 
Be Peo Men? 1 “© Tw (€, 0) v2 (, 0) 
- * a — dé (5’) 
B: Pio Mio? x 0 x—E& x+28,-—¢é 
forx > 0. v(x, 0) = 0 forx < 0. 


The solution of Eq. (5’) together with 


1 ”  t (€,0) (x — &) 
U2 (x, y) = dé 
Br Jo (x — &)? + Be? y? 
py {* v2 (é, 0) 
v2 (xX, y) = es oh, = 
r Jy (x — £)? + Byty? 


determines the subsonic flow 

The functions f(x + 6,) and f(x — 8,) may be related to the 
values of uw. and v on the interface, and the supersonic flow in- 
side the jet may be determined from Eqs. (3’) and (4’). 
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Author’s Reply 


S. |. Pai 
University of Maryland 
March 17, 1952 


—sSsangad SPEAKING, as Messrs. Klunker and Harder pointed 
out, my solution for the case of subsonic uniform stream 
holds true only far downstream of the nozzle. As far as the sub- 
sonic stream is concerned, the present problem is similar to the 
subsonic uniform flow over a wall, which is straight for — » < 
x <0 and which is wavy for 0 <x +o. The only difference in 
these two cases is that the wavy forin in my problem is not given 
but should be determined by the condition on the boundary of 
the supersonic jet. In my analysis I calculate this wavy form 
downstream approximately by assuming a solution of the form 
of Eq. (9) of my paper,! which is implicitly for the flow condition 
far downstream from the nozzle and which should satisfy only the 
boundary conditions given in my paper. Since I am interested 
mainly in whether the supersonic jet in a uniform stream has 
periodic structure or not, my analysis is sufficiently accurate for 
the purpose, and the main conclusions of my simplified analysis 
that the supersonic jet in uniform subsonic stream has almost 
periodic structure and that the approximate wave length in- 
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creases with the Mach Number of the surrounding stream are _ profile, which determines the point of separation. Both are cal- 
correct. However, if one wants to calculate the flow field up- culated by simple quadrature. 

stream or near the exit of the nozzle, more refined theory should From Prandtl’s boundary-layer equations one obtains by a 
be used. well-known procedure the equations for momentum and energy 
My analysis for a supersonic jet in a supersonic uniform stream __ in the following form for the case of axial symmetry: ’ 


holds true for all x, because in a supersonic stream the disturbance (d/dx) (U?R8) 8 dU a 
has no influence upstream. TR 4 7 ds = oF 
REFERENCE (equation of momentum ); 
1 Pai, S. I., On Supersonic Flow of a Two-Dimensional Jet in Uniform (d/dx) ( U%RAS) : d+t 
Stream, Journal of the Aeronautical Science, Vol. 19, No. 1, pp. 61-65, Janu U3R =2 pU3 2 


ary, 1952 
(equation of energy). 

These equations resolve into those of the two-dimensional cas 
by omitting the radius R(x) of the cross section of the body of 
revolution. U(x) is the potential flow, #(x) is the momentum 
loss thickness, 5* is the displacement thickness, 6** is the energy 
loss thickness, H = 6*/3, H = 6**/3, ro is the shearing stress 
at the wall, and (d + ¢) is the work done by the shearing stress 

Subtracting Eq. (2) from Eq. (1), one gets 


An Approximate Method for the Calculation of 
the Laminar and Turbulent Boundary Layer 
by Simple Quadrature for Two-Dimensional 
and Axially Symmetric Flow 


P dA w= odU i d+t Aa” 

= — |) é — 3 
E. Truckenbrodt ° dx U dx pU3 pU? 
Dr.-Ing., Institute. of Fluid Mechanics, Technical University, Both the shearing stress, r)/pU*, and the work of the shearing 


Braunschweig, German = : : : e 
9, y stress, (d + t)/pU%, can be given as functions for the dimension 


January 15, 1952 less quantities Udi/y and H. Furthermore, the velocity profiles 





of the boundary layer are assumed to be one-parametric, and, in 


dimensional and axially symmetric case, valid for laminar flow, as equation, Eq. (2), and not from the equation of momentum, 
well as for turbulent flow. The method consists of calculating Eq. (1), as in most other methods. This brings some advantages 


‘ _— NUMEROUS APPROXIMATE METHODS! * for the calculation consequence of this, there exists a general relation between H 
t of the laminar and turbulent boundary layer have been and A. With regard to this, Eqs. (2) and (3) are two equa 
4 treated from a common point of view. In this way it has been tions for the two unknown quantities 3(x) and A(x). The mo 
~ . . . . . . 

fl possible to establish a uniform approximate method for the two- mentum loss thickness is calculated from the energy loss 
- 

je 

a 


the momentum thickness and a shape parameter of the velocity for carrying out the integration. 


With some simplifications, one gets, for the momentum loss thickness, 


r i+ wn Pahl U 3 + 2n R l1+un x’ 1/(1 + n) 
le+@) "fl, Ga) G4] 
3 2 x /l Ue 1 l 4 


1 


—_- SEYMTTE 


l (U/U.e){R/l) 
where the constant of integration is 
GC. = [((Ui/Ue)*R:/l) (3/1) |' +” 5 


U,, Ri, and 3; are the values in the point x = x). 

Here, cy means the coefficient of skin friction of the flat plate. There is m = 1 for laminar flow and » = 1/6 for turbulent flow. Ir 
the two-dimensional case, R has to be omitted in Eqs. (4) and (5). Eq. (4) can also be used for rough walls by inserting the cy value 
of the rough plate. 

For the convenience of calculation, it was useful to introduce, instead of 7, the parameter 


L= JdAH -1)A 6 
This shape parameter of the velocity profile is calculated from Eq. (3). With some simplifications, one obtains 
g U 1 é U e 
L=~*L,+in— + b — In — | de’ 
€ Ui g ve U; 


where ¢ is obtained from 


rom [o a (2) 4 os ‘( U ) + "(Ay +n P “| . 
2 ie ee i ; 


Here, b and c have the following values: 
Laminar: b = 0; c= 8.0 


_ Ud 
Turbulent: b = 0.07lg — 0.23; c = 4.0 
v 


The calculation of — can be done very simply, because the expression in brackets in Eq. (8) also occurs in Eq. (4 
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The whole calculation is done in the following way: . At first, 
» one gets the momentum loss thickness from Eq. (4); then the 
T shape parameter L from Eq. (7) in connection with Eq. (8). 





Separation occurs in the laminar case for —L > 0.018 and in the 
024 -— tt Truckenbradt [10] —— turbulent case for —L > 0.13 to 0.18. 
Howarth [6, . — ; ,; 
(a) The applicability of the method has been shown in the laminar 

017 —- case for the Howarth* flow U(x) = U.a[l — (x//)], Fig. 1, and 
in the turbulent case‘ for the airfoil NACA 65(216)-222 at an 
7. | , 7 angle of incidence of a = 10.1°, Fig. 2. In both cases, there is a 

002 006 006 008 QO le Qt l comparison with other methods and with measurements. More 
details of our calculations will be published shortly." 


=e — ; ht ete 
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Fic. 1. The laminar boundary layer of the Howarth flow 
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U (x) = Uy 1 — (x/))|. 0 = exact solution, Howarth’; und turbulent Reibungschlichten bei ebener und 
Strémung, to be published in the Ingenieur Archiv 


++10° 
40 


= approximate solution, Truckenbrodt.” 
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Suggested Way of Comparing Overall Efficiency 
of Helicopters 


W. Z. Stepniewski 
Director of Aerodynamics, Piasecki Helicopter Corporation 
February 19, 19592 


So STUDENTS OF THE HELICOPTER still consider the figure 

of merit of the rotor as a proper measure of the efficiency 
with which energy brought into the cycle is transformed into the 
useful work of creating thrust. This approach of evaluating all 
rotary aircraft on the basis of the aerodynamic merits of the rotor 
alone can be justified to some extent within selected groups of 
mechanically driven helicopters equipped with engines of similar 
specific fuel consumption. However, in such cases as jet heli- 
copters, the variations of propulsive efficiency of jets and the 
thermodynamic efficiencies of cycles used may be so great that 
the whole approach based on the figure of merit evaluation be- 
comes of doubtful value. In order to establish a common de- 
nominator for the evaluation of the overall energy efficiency (75. ) 
of all types of helicopters, it is suggested that it should be ex- 
pressed as the ratio of the ideal power (7) required to produce a 
given thrust to the potential power (7) represented by the heat 
energy of the fuel actually used in a unit of time, 


Noe = Pia IP not (1) 


Denoting hourly fuel consumption of the helicopter by F.C. 
(Ibs. per hour) and assuming that heat value of presently used 
fuels is 18,500 B.t.u. per lb., the potential heat power can be ex- 
pressed (in horsepower) as P,., = 7.27 F.C. In hovering, the 
ideal power is simply expressed in terms of gross weight (W), disc 
loading (w), and air density (p), and the overall efficiency be- 
comes 

toe, = WV w/p/5,640(F.C.), (2) 
where the disc loading w (lbs. per sq.ft.) is based on the total 
swept area of the rotor (s), and (F.C. ), is the hourly fuel consump- 
tion in hovering (Ibs. per hour). 

In forward flight, the overall efficiency of different helicopters 
should either be evaluated at some definite speeds (say 100, 120 
m.p.h.), or still better at the flying speed, Vu »:., corresponding 
to the longest range. Substituting in Eq. (1) for P,,, its expres- 
sion in horizontal flight and simplifying, the following expression 
for overall efficiency at optimum flying speed (in m.p.h.) is ob- 
tained, 

Nosop. ™ W?/11,80027R2o( F.C. opt. Vopr. (3) 

Having the values of ,, (say from statistics or general consider- 
ations) for various types of helicopters, together with their use- 
ful load/gross weight ratios, it is possible to compare them from 
the point of view of their pay load carrying capacity in hovering 
and in forward flight. In the case of hovering this can be done 
by calculating the hourly fuel consumption from Eq. (2) and 
comparing it with the amount of fuel available at a given pay 
load. Similarly, in the case of range, (F.C. )opt. is computed 
from Eq. (3); then from the known flying speed, Vopt. and 
amount of fuel at a given pay load, the corresponding range is ob- 
tained. Through this procedure curves showing pay load vs. 
time in hovering or payload vs. range can be established. 

In order to g've some feeling as to the order of magnitude of 
the moe, and noe, Values, these efficiencies are briefly discu sed 
for the example of mechanically driven helicopters. 

Hourly fuel consumption in hovering of this type of helicopter 


can be expressed as 


s.f.c mW I 
(F.C.), = kW — + rR%pz,, V3 (4 
=: ( V2 8 en ) / 


p 


where new symbols are: 7; = the transmission efficiency (jn. 
cluding cooling and antitorque losses), k = the ratio between ideal 
and actual induced power in hovering (see page 94, reference |), 
o = rotor solidity, €;, = equivalent profile drag coefficient, V, = 
tip speed (ft. per sec.), and s.f.c. = specific fuel consumption 
of the engine (lbs. per b.hp. per hour). It is easy to show that 
(see page 104, reference 1) maximum power loading will be ob- 
tained when the second term in the brackets of Eq. (4) equals 
one-half of the first one or, in other words, when the disc load- 
ing is 

w = (1/2)pVPEW(ot,,/k)? (5 


Substituting Eq. (4) into Eq. (2) and assuming that the dis 
loading w is as given by Eq. (5), the following simple expression 
for the overall efficiency in hovering is obtained :* 


Noe, = 7 7.27(|k + (1/2)|s.f.c. (6 

The following numerical values may be considered as representa. 

tive for a single-rotored helicopter: , = 0.84, k = 1.1, and sfc. 

= ().5 lb. per b.hp. per hour; hence, 7, ¢, = 9.145. It may be of 
' } 


interest to note that for a ram-jet helicopter this efficiency will 
probably be in the neighborhood of Nor, = 0.015 or less, while for g 
well-designed pressure jet it may reach or exceed 7, = 0.08 

The engine power, P?, required in horizontal flight of mechani- 
cally driven helicopter when expressed in ft.lbs. per sec. is 


1 e 2. kW? 1 
P= pfV3 + Py» + : 
me \2 2rR2p J 


where the flying speed V is in ft. per sec. and the new symbols are 
f is equivalent flat plate area (including drag component due to 
the rotor profile drag), P,, is the rotor profile power (in ft.lbs 
per sec.), and & is, in this case, the ratio of the ideal to the actual 
induced power in forward flight. For a constant specific fuel 
consumption, the best range will be obtained when P/V becomes 
aminimum. This condition is fulfilled when 


Substituting Eq. (8) into Eq. (7), changing units of power to 
hp., and multiplying the thus-transformed Eq. (7) by s.f.c., the 
hourly fuel consumption is obtained, and the overall efficiency at 


optimum flying speed becomes 


14.54] k ane te f 
= wo ," S22. ({ 
ee V(W/P yyw) 


where (W’/P,,,) may be called the profile power loading (Ibs. per 
hp.) and Vopr. isin m.p.h. 

For a general discussion it would be more desirable to express 
Vopt. in terms of design parameters. This can be done by solv 
ing Eq. (8) for V. Unfortunately, because this equation is of the 
fourth degree, no simple exact solution can be obtained. How- 
ever, ai approximate value of Vopt. can be obtained by neglecting 
the Pp, term. Then Vopt. in m_p.h. would be 


> ~ 4 ° / 
Vopt. = 1.47V kww,/p? (10 


where wy = W’/f (Ibs. per sq.ft.) may be called the equivalent 
flat plate area loading. Substituting Eq. (10) into Eq. (9), the 
following approximate formula is obtained: 


825 ‘ kw ,p? 


some, = 11/14.54 | b+ sfc. (Il 
ges er (W/Pr) Nw 


For a single-rotored helicopter with the following characteris 


* It should be noticed that quite often for structural and other reasons 
the disc loading of practical helicopters deviates from its optimum value 
will be lower than given by Eq. (6 


as given by Eq. (5); hence, Noe, 
h 
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tics: w = 5 Ibs. per sq.ft., W/P pr = 40 Ibs. per hp., ws = 800 Ibs. 
per sq.ft. (very clean machine) and assuming that, in the con- 
sidered forward flight, k = 1.1, 7 = 0.86, and s.f.c. = 0.45 Ib 


per hp. per hour, the optimum overall efficiency in horizontal 
se 2 * 
fight would be ee, = 0.043. 
REFERENCE 
Stepniewski, W Z., Introduction to Helicopter Aerodynamics, Rotorcraft 


Publishing Committee, Morton, Pa., 1951 


* The above outline represents only the private opinion of the writer 


On the Stability of Plane Parallel Laminar 
Flows to Two- and Three-Dimensional 
Disturbances* 


Martin Lessen 
The Pennsylvania State College, State College, Pa. 


February 27, 1952 


( FLOW EQUATIONS of a viscous compressible fluid may be 
given as follows: 


Continuity: 


a + (puke = 0 
Force: | 
pit; = —pi + (mii, &)e + (CHU k, i) — (2/8) (Cutts, «&,)i 


p (1) 


Energy: 


pep! = (IKT:) x — [Pp + (2/3)uur, eur, a + 

MC UK, Met + UK, i,k) | 
Stale: 
o=pp,7T) w=u(T); K = K(T) 


where the variables are defined as 


x; = rectangular Cartesian space coordinates 
{ = time coordinate 

u; = velocity component in 7 direction 

p = density of fluid 

p = pressure 

TY = temperature 

un = coefficient of viscosity (absolute) 

K = conductivity (thermal) 


If a plane parallel steady flow with a small time variable three- 


dimensional disturbance superimposed on it is assumed, such that 


uy = U(xe) + u(x) exp [i(ax, + Bx; — act)| 
U2. = (xX) exp [i(ax1 + Bxy — act)] 

u, = w(xe) exp [i(ax; + Bx; — act)] 

p = pA(xs) + p(x2) exp [i(ax, + Bx; — act)] 

pb = pt p(x2) exp [i(ax, + Bx3s — act)] 

T = T(x2) + T(x.) exp [i(ax + Bx3 — act)] 
Me = p(X.) + p(x) exp [i( ax, + Bx; — act) | 


K (x2) + K (x) exp [i(ax,; + Bx; — act) | 


ll 


the equations for the disturbance may be obtained by elimi- 
nation of the steady-state terms from Eq. (1) as satisfying them- 
selves and, in addition, dropping the terms nonlinear in the dis- 
turbance quantities in the usual manner. It should be noted 
at this point that, generally, absolutely parallel flows cannot 
exist. 


* This work was carried out under O.N.R. Contract NONR- 


656(01). 
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Continuity: 
ia(U — c)p + pliau + 0’ + iBw) + p'v = 0 (2a) 
Force: 
_— F ; 4 
plia(U — chu + U'v] = —iap +a] u” — 3 a? + B? Ju + 
la ? ap ? $ 2 , 
3 v — 3 v| + p'(u’ + tav) + U"p + U'p’ § (2b) 
; : , 4 
blia( UU — cw] = —p’ +2 2” + (a* + 8*)v + 
la , ip : P 4 : 2 ’ 
uo + wita2e v,— tau — pw + tal "u (2c) 
3 3 3 j 3 
. s * . ” 4 
dlia(U — chw) = —iBp + pl] w” —[ a? + 3 B? Jw — 
ap iB ee ; 
u v’ | + g'(ipv +w’) (2d) 
3 3 
Energy: 
icp[ia(U — c)T + Tv] = K [T" — (a? + B?)T) + 
R'T' + KT" + K'T’ — pliau + v' + ifw) + 2iU’ X 
(u’ + iav) + w(U’)? (2e) 


The following transformation of dependent variable will now 
be tried similar to that of Squire? after multiplying Eq. (2b) by 
a, Eq. (2d) by 8. and adding: 


au + Bw = au 


a® + Bp = q? 
(a /a)p = »p 
(a/ajp =h 


The continuity and force equations transform directly to two- 
dimensional form with the exception of the viscosity-disturbance 
terms. The energy equation also transforms to two-dimensional 
form except for the terms 24U’(u’ + iav) and a viscosity dis- 
turbance term. It should be noted at this point that, even for 
a classical incompressible fluid, the aforementioned energy equa- 
tion terms do not satisfy the Squire transformation. 

In the region where U is not approximately equal to c, the 
solution of the disturbance equations is given primarily by the 
inviscid terms. These terms obey the Squire transformation. 
In the “inner friction’ 
the temperature (and therefore density) disturbance has an unim- 
portant effect on the eigenvalue problem,' the viscosity disturb- 
Therefore, for 


’ 


region where U =~ ¢, if it is assumed that 


ance and the energy equation can be neglected. 
this region also, the Squire transformation is applicable. 

Thus, if instability occurs for a three-dimensional disturbance, 
it will occur for a two-dimensional disturbance at a lower Reyn- 
olds Number for the same steady-state flow, sincea@ > a. It is 
therefore only necessary to examine the stability of parallel flows 
to two-dimensional disturbances. 

The equations for a two-dimensional disturbance follow. 


Continuity: 
ia(U — c)p + pliau + v') + pv = 0 
Force: 


plia(U — chu + Uv] 
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Energy: 
Acy[tia(U — c)T + Tv] = — pliau + v’) 

If it is assumed that the temperature (and density) disturb- 
ances have a higher order effect on the eigenvalue problem 
throughout the flow régime, the energy equation can be dropped 
entirely. This is rigorous for the case of an incompressible fluid 
with temperature-variable viscosity and for the laminar mixing 
of dissimilar fluids. 
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A High-Speed Schlieren Technique for Investigation of Aerodynamic 
Transients 
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flow problems of any description. In particular, it 
lends itself most admirably to the study of short-period 
transients. Other cameras are available*® which permit 
the investigation of lower frequency effects over a 
longer period of time. Fig. 8 shows a series of photo- 
graphs obtained during the starting process of a super- 
critical supersonic diffuser showing the effect of a nor- 
mal shock wave propagating through a steady super- 
sonic flow in the.shock configuration in that flow. The 
exposure time of each frame is in the order of 2 to 4 


microsec., and the frequency of exposure is about 
7,000 thousand frames per sec. 
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